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Abstract 

I first give an overview of the thesis and Matrix Product States (MPS) representation 
of quantum spin systems on a line with an improvement on the notation. 

The rest of this thesis is divided into two parts. The first part is devoted to 
eigenvalues of quantum many-body systems (QMBS). I introduce Isotropic Entan- 
glement (IE), which draws from various tools in random matrix theory and free 
probability theory to accurately approximate the eigenvalue distribution of QMBS 
on a line with generic interactions. We then list some open related problems. Next, I 
discuss the eigenvalue distribution of one particle hopping random Schrodinger oper- 
ator in one dimension from free probability theory in context of the Anderson model. 

The second part is devoted to ground states and gap of QMBS. I first give the 
necessary background on frustration free Hamiltonians, real and imaginary time evo- 
lution of quantum spin systems on a line within MPS representation and a numerical 
implementation. I then prove the degeneracy and unfrustration condition for quan- 
tum spin chains with generic local interactions, including corrections to our earlier 
assertions. Following this, I summarize my efforts in proving lower bounds for the 
entanglement of the ground states, which includes partial new results, with the hope 
that they inspire future work resulting in solving the conjecture given therein. Next 
I discuss two interesting measure zero examples where the Hamiltonians are carefully 
constructed to give unique ground states with high entanglement. One of the exam- 
ples (i.e., d = 4) has not appeared elsewhere. In particular, we calculate the Schmidt 
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numbers exactly, entanglement entropies and introduce a novel technique for calcu- 
lating the gap which may be of independent interest. The last chapter elaborates 
on one of the measure zero examples (i.e.. d = 3) which is the first example of a 
frustration free translation- invariant spin-1 chain that has a unique highly entangled 
ground state and exhibits signatures of a critical behavior. 

Thesis Supervisor: Peter W. Shor 
Title: Professor 
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Chapter 1 



Many-Body Physics and an 
Overview 

Physics is concerned with specification and evolution in time of the state of particles 
given the laws of interactions they are subjected to. This thesis is devoted to better 
understanding of quantum aspects of Many-Body Systems. 

1.1 Phase Space 

The phase space of a classical system of N particles is a 6A^ dimensional space where 
each particle contributes 3 spatial coordinates and 3 momenta. A point in this space 
fully specifies the state of the system at any given time; the motion of this point 
in time specifies the time evolution of the system. Even if the laws of interaction 
are precisely formulated, the analytical solution of the equations of motion, given 
the exact initial state, can be daunting The complexity is due to interactions. 

""^Uncertainties in the initial conditions can give rise to further compHcations caused by chaotic 
behavior for positive Lyapunov exponents 
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For example general solution of the 3— body problem has been an open problem 
for roughly 350 years[lj. Despite the interaction being 2— body, the correlation in 
time, of the distances between pairs of masses makes the analytical solution hard to 
obtain. In such a case, one can resort to computational methods to simulate and 
make predictions with controllable accuracy by keeping track of 6A^ real parameters. 

The state of a quantum system of interacting particles is defined by a ray in the 
Hilbert space, the dimension of which is a multiplicative function of the dimension 
of individual Hilbert spaces (see [21 |3l section 2.1] for nice expositions of quantum 
theory) . Mathematically, 



where we denote the Hilbert space of the i particle by Tii and, following Dirac's 
notation [21 section 20], the pure state of the system by the vector ip). We restore | 
only when ambiguity of the label of the vector with linear operators preceding it may 
arise. To fully specify the system one needs to specify ci^ complex numbers (assuming 
each particle has d degrees of freedom), which makes the study of quantum many- 
body systems (QMBS) computationally intractable. This, compared to simulation of 
classical many-body systems, is an additional obstacle we face. Hence, the complexity 
in studying QMBS is in interactions as well as state specification. 

It is worth mentioning that there are quantum systems that are particularly 
simple, without having as simple a classical analogue. An example of which is a 
quantum bit (quhit), whose state consists of only two points in the phase space. 

In this thesis we confine ourselves to finite dimensional Hilbert spaces where a 
general pure state of an A^— body problem each of which having d states is 



N 



^)e(g)-H. 




i=l 
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^) = ^^^-^^ Ii)®12)®---®in). (1.1.2) 

where repeated indices are summed over. 

In classical physics one can pick a single particle among interacting particles, 
specify its state at some time and predict its evolution subject to the fields of the 
other particles. If after some time the fields impressed on this particle by the remain- 
ing particles diminish, the particle becomes free and uncorrelated from the remaining 
A^ — 1 particles. In contrast, in quantum mechanics, interaction can lead to quantum 
correlations (entanglement) that persist even if the interaction strength diminishes, 
by say taking the particles far apart. For example two electrons can interact for 
some time and end up in the entangled state 1^1,2) = 00) + 11) |39l see discussions 
on EPR pairs]. This peculiar feature of quantum mechanics, as of yet unexplained 
by classical physics, is a radical departure from the latter. Quantum computation 
and quantum information science make use of entanglement as a resource to do tasks 
that are classically difficult or impossible to do in reasonable time such as |H |5] . 

Whether, for a given real linear operator, any eigenvalues and eigenvectors exist, 
and if so, how to find them is in general very difficult to answer [2, p. 32]. 

In this thesis I mainly focus on QMBS systems with Hamiltonians (except in 
chapter 3) 

N-l 

if = V-.-i. (1.1.3) 



See Figure 1.1.1 for the sparsity pattern of H. 



Eq. 1.1.2 is the most general equation for a pure state of N particles irrespective 



^One could further simplify the notation by denoting the state by ) = -0) = t^'i"'*™ ) h )i2 '■')«« j 
but I do not think it is worth the trade offs. 
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Figure 1.1.1: Sparsity pattern of H. The nonzero elements for k odd are shown in 
blue and k even in red. 



of the type of interactions or configuration (see Fig. 1.3.1). However, the state of a 
system for any given problem has inherent constraints such as an underlying lattice 
that reduces the spatial degrees of freedom, a range of interaction, or underlying 
symmetries that the system has to obey. It is therefore, reasonable to suspect that 
the physical properties of any given problem could be well approximated by a re- 
duced number of parameter say ~ poly (A^) if one understood the effective degrees 
of freedom well enough. 

Non-commutativity of the interaction in QMBS is responsible for the richness 
of possibilities such as various phases of matter and quantum computing. Simulta- 
neously, it is accountable for formidability of finding eigenvalues and corresponding 
eigenvectors. The exact computation of eigenvalues alone, on a line, has shown to be 
QMA-Hard [31j. Energy eigenvalue distributions are needed for calculating the par- 
tition function and calculating other transport properties (Part I). The eigenvectors 
specify the corresponding states of the system. Bulk of matter usually finds itself in 
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its lowest energy [2H1 P- 48]; hence, the urge to comprehend the low lying states in 
condensed matter research (Part 11). 

1.2 Part I: Eigenvalues 

Consider the general problem of predicting the eigenvalue distribution of sums of 
matrices from the known distribution of summands. In general this is impossible to 
do without further information about the eigenvectors. However, any progress in this 
direction is extremely desirable as many problems are modeled by non-commuting 
matrices. For example, the Schrodinger equation has a kinetic term and a potential 
term; often the former is diagonalizable in the Fourier and the latter in position 
space. However, the sum does not have an obvious global basis nor a distribution 
that can trivially be inferred from the known pieces. 

For the sake of concreteness suppose we are interested in the distribution of the 
random matrix M = Mi + M2 where the distributions of Mi and M2 are known. 
There are two special cases worth considering: 

1. The summands commute. In this case, one can find a simultaneous set of 
eigenvectors that diagonalizes each of the summands. In the language of linear 
algebra of diagonal matrices, the eigenvalues add. When the eigenvalues are 
random, this connects us to the familiar classical probability theory where the 
distribution of the sum is obtained by a convolution of the distribution of the 
summands. 

2. The summands are in generic positions. In this case, the eigenvectors of the 
summands are in generic positions. It is a fascinating fact that this case also has 
a known analytical answer given by Free Probability Theory (FPT) |84j |48] . 
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The eigenvalue distribution of M is given by the "free convolution" of the 
distributions of Mi and M2. 

There are many interesting questions that one can ask. How "free" are general non- 
commuting matrices? What is the relationship between commutation relation and 
freeness of two matrices? To what extent does the Fourier matrix act generic? Can 
a large class of non-commuting matrices be analyzed using a convex combination of 
the two extreme cases discussed above (see Isotropic Entanglement)? 



Suppose the local terms in Eq. 1.1.3 are generic (i.e., random matrices), can we 
utilize the existing tools of random matrix theory to capture the eigenvalue distribu- 
tion of H given the distribution of if;^/+i's? Despite, the local terms being generic, 
H is non-generic. The number of random parameters grow polynomially with N 
whereas H is dimensional. Fraction of sparsity of if is < {N — 1) d~^^~'^^ (Figure 



1.1.1). 



Since the exact evaluation of the density of H is very difficult [31], one can use 
two known approximations. As far as parameter counting is concerned the quantum 



problem falls nicely in between the two extreme case (Figure 1.2.1 and Isotropic 
Entanglement). 

Chapter 2 gives a detailed description of Isotropic Entanglement which gives the 



distribution of Eq. 1.1.3 and elaborates on the discussion of this section. 

Chapter 3 is on the distribution of one particle hopping random Schrodinger 
equation and Anderson model. 



1.3 Part II: Eigenvectors 

Is it possible to capture the essential physics of the system accurately enough with 
an efficient simulation with a much smaller x ~ Polj{N), spanning only a small part 
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The Unitary Group U(d^) 

Isotropic Quantum Classical 




Number of Random Parameters 



Figure 1.2.1: Parameter counting 



of the full Hilbert space of the system? In our case, the qudits are arranged on a 
ID lattice and only have nearest-neighbor interactions. We could thus expect that 
a reduced space might suffice for our needs. This concept is common for the vari- 
ous approaches proposed for efficient (tractable on a classical computer) numerical 
investigation of quantum many body systems such as the Density Matrix Renor- 
mahzation Group [35], Matrix Product States |37], Tensor Product States [36] and 
projected entangled pair states (PEPSs) |100j . For gapped one dimensional systems 
MPS ansatz is proved to suffice |112j . Consider a general Hamiltonian for ID open 



chain with generic local interactions as given by Eq. 1.1.3, where each is a 

(P X matrix of rank r and the total number of particles is A^. 

It is interesting to ask under what circumstances can there be a degeneracy of 
ground states? Moreover, when is the ground state of the whole system {H in Eq. 
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1.1.3) also the local ground state of all Hi i^is ; i.e., the system is "unfrustrated" ? 
We answered these questions for ID spin systems with generic interactions |133] 
(Chapter 6). We found that in the regime where r < (P/4 the system is unfrustrated 
with many ground states; moreover for r < there can be product states among the 
ground states. For sufficiently large N the system is frustrated for r > (P/4. 

The next natural question is: how entangled are the ground states in the regime 
d < r < (P/A7 The entanglement can be quantified by the Schmidt rank (see MPS 
below). We call a state highly entangled if its Schmidt rank is exponentially large in 
the number of particles. In this regime, it is straightforward to show that among the 
many ground states, there are no product states and that there exist states with high 
amount of entanglement with probability one. The latter can be shown for example 
using results of algebraic geometry [TT] [TU] . I have been trying to show that all the 
ground states, in this regime, have Schmidt ranks that are exponentially large with 
probability one. Despite some partial results, the goal has not been fulfilled. 

Using a genericity argument one can show that, to prove results in the generic 
case, it is sufficient to find an example of local terms whose ground states all have 
large Schmidt ranks. We have not yet succeeded in finding such examples in the 
regime of interest d < r < d'^/A. There are, however, interesting examples for which 
there is a unique ground state with exponentially large entanglement in the frustrated 
regime, i.e., probability zero case (Chapters 8 and 9). 

Below I give a background on Matrix Product States on with an improvement on 
the notation. 

Chapter 4 discusses the unfrustration condition and a numerical code I developed 
to study spin systems on a line with local interactions and without translational 
invariance. I find the ground states using imaginary time evolution. I then provide 
the proofs and corrections of our previous work regarding the unfrustration condition 
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and degeneracy of quantum spin chains with generic local interactions. 

Chapter 5 summarizes various attempts I made in proving a lower bound on the 
Schmidt rank of the ground states of generic spin chains. It includes unpublished 
results and two ways one can potentially prove the conjecture given there. 

Chapter 6 discusses two interesting measure zero examples {d = 3 and d = 4). 
I include the combinatorial techniques for calculating the entanglement entropies. 
The d = 4 example has not been published elsewhere. 

Chapter 7 elaborates on the d = 3 example of Chapter 6 and has a novel technique 
for calculating the gap that may be of independent interest. 

1.3.1 SVD and Matrix Product States (MPS) on a Line 

In this section, in order to avoid confusion, we restore the summation symbols. As 
stated earlier the state of a composite system is a vector in the tensor product of the 
Hilbert spaces of the constituents. Suppose we have the pure state of a composite 
system, how can we express it in terms of the pure states of the constituents? This 
is answered by singular value decomposition (SVD). This application of SVD in 
quantum information theory is called Schmidt Decomposition p] . |39j . 

Theorem 1. (Schmidt Decomposition) Let Hi and if 2 be Hilbert spaces of dimen- 
sions n and m respectively. Assume m < n. For any vector 1^12) G iii (g) H2, there 
exists orthonormal sets {r]i),--- ,rin)} C Hi and {ui),--- ,ujm)} C H2 such that 
^) ~ Siii -^i Vi) ® ^i) , where Aj are non-negative and , as a set, uniquely deter- 
mined by ip) . 

The number of nonzero A's are called the Schmidt rank denoted by x- 

Comment: % < m < n; the Schmidt rank is no greater than the minimum of the 
dimensions of the two Hilbert spaces. 
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Figure 1.3.1: Some general interaction quantum spins whose pure state is repre- 
sented by Eq. |1.1.2 



Comment: Schmidt Decomposition can be thought of as an expansion of a vector 
in bases of the subsystems. The need for having more than one expansion coefficient 
(i.e., X > 1) indicates that the state is not separable in any basis (i.e., subsystems 
are entangled). A quantifier for entanglement is therefore x- 

DMRG, natural representation of which is MPS, beautifully utilizes Schmidt 
Decomposition to capture low energy properties of QMBS on a line or trees |33l |3lj [7] . 
Though any state can be expresses as a MPS, the time evolution is only naturally 
implemented for the line or a tree. Some attractive features are 1. MPS gives a 
local description of QMBS on a line 2. MPS allows for a systematic truncation of 
the Hilbert space to capture the low energy properties with controllable accuracy. 
Below I give a derivation of MPS similar to [96j and improve the notation by making 



it more compact (Eq. 1.3.8). 



Suppose we have a chain with sites and take 1 < n < A^. First, perform a 
Schmidt decomposition of the chain between sites n — 1 and n as 



Xn-l 



^) = E ^£"1^ 1 ® <^«.-i)n,-,iV, (1.3.1) 
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where the states on the left and on the right of the division form an orthonormal 
bases for the respective subsystems of the state of ip). The Schmidt rank Xn-i is the 
minimum number of terms required in this decomposition. Recall that Xn-i is at 
most equal to the minimum dimension of the Hilbert spaces of the two subsystems. 

Next, the Schmidt decomposition for a split between n and n + 1 gives 



^) = E-^S K:)l,-n ® 0.„)n+l,-,N. (1.3.2) 



These two decompositions describe the same state Eq. 1.1.2 allowing us to combine 



them by expressing the basis of the subsystem n, ■ ■ ■ , iV as 



d Xn 

0a„_,)n,...,iV = E E ^ailajal ® ^« Jn+l,- ,Ar, (1-3.3) 

i„ = l On = l 

where we inserted Aq"^ for convenience. This gives us the tensor r'"' which carries an 
index i„ corresponding to the physical states i„) of the n spin, and indices 
and a„ corresponding to the two consecutive divisions of the system. Since 4>a„-i) 
and 6a„) are orthonormal states, the vectors A and tensors F obey the following 



normalization conditions. From Eq. 1.3.2 we have 



Xn 



a„=l 



while 1.3.3 implies 



d Xn 

= E E r£;„^HrH-,„„AN = 5^„_,,^,_^, (1.3.5) 

in = l an = l 

and 
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^a„_i/l,..., «-l ^a„_, /«,..., 



n-1 




n 



n+1 




n 



^a„ )«+!,..., 



N 




Figure 1.3.2: Decompositions that define F'"! in the derivation of MPS. 



d Xn-l 

, If) \ = \^ A [n-l]pN,in* X [n-l]p[n]-'n ^ X 

i„=l a„_i = l 



;i.3.6) 



One can do what we just did for every site \ <n < N and get the MPS represen- 
tation of the spin chain, denoting open boundary conditions and periodic boundary 
condition by OBC and PBC respectively, 



^2112. ..ijv 



J- l,ai J- ai,Q2 ■ ■ ■ ajv-1,1 
oi,...,ojv_i=l 

X 

E-pn,[l] ■p«2,[2] _ _ _ -pjjv,[Af] PRP 



;i.3.7) 



"0:"lv;"iV-l=l 



Comments: It is customary to absorb the A's into F's and omit them as we did 
in the foregoing equation. The upper hmit x = niaxi<„<Ar {Xn}- 



I beheve the MPS notation given by Eq. 1.3.7 can be improved 
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TV d 



^) = ^ <^ (8) E r fe) ? 0)®"^ OBC (1.3.8) 

p=i 1^=1 J 

^) = Ti^V I (g) 5^ r (tp) a'^ i 0)®^ PBC 

p=l ip=l J 

where 0"*^ defined by a*'' |0) = ip) are the generahzed Pauh operators, each F (ip) = 
pbl.^p —s. Tal'yi,ap, for a given ip, is a x x % matrix and V denotes an ordering with 
respect to p of the tensor products. The subscript x the trace reminds us that the 
trace is on x x % part of F's and not the physical indices ip. Lastly, we can simplify 
notation by assuming repeated indices are summed over to get 



ij) = Ti^V !^(^^T (ip) a'^^ Of^ PBC (1.3.9) 

Note that we now need at most ~ Ndx^ parameters to specify any state. There 
are orthogonality conditions on F's that further reduce the number of independent 
parameters needed. 



1.3.2 MPS on a Tree 

The MPS representation relies on splitting the system into two subsystems and 
making use of Schmidt Decomposition which applies for loop-less configuration of 
spins such as chains and trees. The derivation above can be generalized (see [HBJ) to 
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a T,a , a ,,a a.a^, 



n-1 V 



[n] 




'a , '^ot 

n—l ' 



Figure 1.3.3: MPS representation 

spins on a fc— child tree. 



Xp 



n^.bondsE^S n^^sitesErXil..,.^ J (1.3.10) 

y ap=l J \ i„ / 

where , ■ ' ' ^ (^n^ are indices corresponding to the k bonds rii, - ■ ■ , coming out 
of site n. Each index appears in F tensors and one A vector. The normahzation 
conditions for a MPS description of a state on a tree are analogous to Eqs. 1.3.4[1.3.5 
and ll.3.6[ We have 



Op 

E E ^^C^:.,>^--->^ r-N!,„„^= (1.3.12) 
and the other variations with n2, - ■ ■ ,nk- 



32 



Part I 
Eigenvalues 
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Chapter 2 



Isotropic Entanglement 



In this first part of tlie tliesis I focus on the density of states or QMBS. This chapter 
treats the eigenvalue distribution of spin chains, though some of the theorems apply 
in higher dimensions. We treat generic local interactions, where by local I mean 
every interaction term acts nontrivially on L consecutive spins. This chapter also 
appears in [H [9] . 



2.1 Elusive Spectra of Hamiltonians 

Random matrix techniques have proved to be fruitful in addressing various problems 
in quantum information theory (QIT) [T0l[32l[TT] and theoretical physics [121 [13]. 
condensed matter physics (CMP), quantum spins with random exchange have been 
extensively studied with surprising phase transition behaviors [T7t [TH [T3] . 

Moreover, constraint satisfaction problems are at the root of complexity theory 
[T5| ITU] . The quantum analogue of the satisfiability problem (QSAT) encodes the 
constraints by a Hamiltonians acting on quantum spins (or qudits) [20]. Interest- 
ing phase transition behaviors have been seen for random QSAT with important 
implications for the Hardness of typical instances[2Il [221 [23]- 
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While the apphcation of the aggregate of eigenvalues may be different in CMP, we 
point out a few examples where they are studied [211 |26l |25l |27] . By comparison, in 
QIT and CMP the ground state and the first few excited states have been well studied 
to date EH ESI ESI EZl E21 ES] • To the best of our knowledge the aggregate has 
not been as well studied perhaps in part because of the complexity [31] and perhaps 
in part because the low energy states have captivated so much interest to date. 

Energy eigenvalue distributions or the density of states (DOS) are needed for 
calculating the partition function [281 p. 14]. The DOS plays an important role in 
the theory of solids, where it is used to calculate various physical properties such 
the internal energy, the density of particles, specific heat capacity, and thermal con- 
ductivity. Quantum Many-Body Systems (QMBS) spectra have been elusive for two 
reasons: 1. The terms that represent the interactions are generally non-commuting. 
This is pronounced for systems with random interactions (e.g., quantum spin glasses 
[111113]). 2. Standard numerical diagonalization is limited by memory and computer 
speed. The exact calculation of the spectrum of interacting QMBS has been shown 
to be difficult [3I|. 

An accurate description of tails of distributions are desirable for CMP. Isotropic 
Entanglement (IE) provides a direct method for obtaining eigenvalue distributions 
of quantum spin systems with generic local interactions and does remarkably well in 
capturing the tails. Indeed interaction is the very source of entanglement generation 
[Sni Section 2.4.1] [ID] which makes QMBS a resource for quantum computation [H] 
but their study a formidable task on a classical computer. 

Suppose we are interested in the eigenvalue distribution of a sum of Hermitian 
matrices M = YliLi ^i- general, Mj cannot be simultaneously diagonalized, con- 
sequently the spectrum of the sum is not the sum of the spectra. Summation mixes 
the entries in a very complicated manner that depends on eigenvectors. Nevertheless, 
it seems possible that a one-parameter approximation might suffice. 

Though we are not restricted to one dimensional chains, for sake of concreteness, 
we investigate interacting d-dimensional quantum spins (qudits) on a line with 
generic interactions. The Hamiltonian is 
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Figure 2.1.1: Odd and even summands can separately be locally diagonalized, but 
not the sum. The overlap of the two subsets at every site generally requires a global 
diagonalization. 



iV-l 

H =J2h'-^ ^ Hi^-,i+L-i®hN-i-(L-i), (2.1.1) 
1=1 

where the local terms Hi ... i^l_i are finite x random matrices. We take the 
case of nearest neighbors interactions, L = 2, unless otherwise specified. 

The eigenvalue distribution of any commuting subset of H such as the terms with 
/ odd (the "odds" ) or / even (the "evens" ) can be obtained using local diagonalization. 
However, the difficulty in approximating the full spectrum of H = H^^^^^ + H^^^^ is 
in summing the odds and the evens because of their overlap at every site. 

The intuition behind IE is that terms with an overlap, such as and if;+i^;+2, 

introduce randomness and mixing through sharing of a site. Namely, the process of 
entanglement generation introduces an isotropicity between the eigenvectors of evens 
and odds that can be harnessed to capture the spectrum. 

Random Matrix Theory (RMT) often takes advantage of eigenvectors with Haar 
measure, the uniform measure on the orthogonal /unitary group. However, the eigen- 
vectors of QMBS have a more special structure (see Eq. 2.3.10). 

Therefore we created a hybrid theory, where we used a finite version of Free 
Probability Theory (FPT) and Classical Probability Theory to capture the eigen- 



value distribution of Eq. 2.1.1 Though such problems can be QMA-complete, our 
examples show that IE provides an accurate picture well beyond what one expects 



from the first four moments alone. The Slider (bottom of Figure 2.2.1) displays the 
proposed mixture p. 
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2.2 The Method of Isotropic Entanglement 



2.2.1 Overview 

We propose a method to compute the "density of states" (DOS) or "eigenvalue 
density" of quantum spin systems with generic local interactions. More generally 
one wishes to compute the DOS of the sum of non-commuting random matrices 
from their, individually known, DOS's. 



We begin with an example in Figure 2.2.2 , where we compare exact diagonaliza- 
tion against two approximations: 

• Dashed grey curve: classical approximation. Notice that it overshoots to the 
right. 

• Solid grey curve: isotropic approximation. Notice that it overshoots to the left. 

• Solid black curve: isotropic entanglement (IE). 



Dots: exact diagonalization of the quantum problem given in Eq. 2.1.1 



The classical approximation ignores eigenvector structure by summing random eigen- 
values uniformly from non-commuting matrices. The dashed curve is the convolution 
of the probability densities of the eigenvalues of each matrix. 

The isotropic approximation assumes that the eigenvectors are in "general posi- 
tion"; that is, we add the two matrices with correct eigenvalue densities but choose 
the eigenvectors from Haar measure. As the matrix size goes to infinity, the resulting 
distribution is the free convolution of the individual distributions |18]. 

The exact diagonalization given by red dots, the dashed and solid grey curves 
have exactly the same first three moments, but differing fourth moments. 

Isotropic Entanglement (IE) is a linear combination of the two approximations 
that is obtained by matching the fourth moments. We show that 1) the fit is better 
than what might be expected by four moments alone, 2) the combination is always 
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+ 



+ 



add random variables = 
Classical convolution of densities 



add random variables = 
Classical convolution of densities 



Several options for adding 
random variables 



Section lll-B,C 



Isotropic 
p = 



Section lll-D 



Section lll-A 



Classical convolution = 
dv^ * 




Quantum convolution = dv'' = dv. ffl„dv„ 

A q b 

Isotropically Entangled 



Classical 
p = l 



Figure 2.2.1: The method of Isotropic Entanglement: Quantum spectra as a convex 
combination of isotropic and classical distributions. The Slider (bottom) indicates 
the p that matches the quantum kurtosis as a function of classical {p = 1) and 



isotropic (p = 0) kurtoses. To simplify we drop the tensor products (Eq. 2.3.7) 
in the local terms (ellipses on top). Note that isotropic and quantum convolution 
depend on multivariate densities for the eigenvalues. 
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N=5, d=2, r=4, trials=500000 
0.05 1 1 1 1 




Eigenvalue 

Figure 2.2.2: The exact diagonalization in dots and IE compared to the two approx- 
imations. The title parameters are explained in the section on numerical results. 

convex for the problems of interest, given by < p < 1 and 3) this convex combi- 
nation is universal depending on the parameter counts of the problem but not the 
eigenvalue densities of the local terms. 

Parameter counts: exponential, polynomial and zero. Because of the locality of 
generic interactions, the complete set of eigenstates has parameter count equal to a 
polynomial in the number of spins, though the dimensionality is exponential. The 
classical and isotropic approximations have zero and exponentially many random 
parameters respectively. This suggests that the problem of interest somehow lies in 
between the two approximations. 

Our work supports a very general principle that one can obtain an accurate 
representation of inherently exponential problems by approximating them with less 
complexity. This realization is at the heart of other recent developments in QMBS 
research such as Matrix Product States [331 |3l], ^ind Density Matrix Renormalization 
Group where the state (usually the ground state of ID chains) can be adequately 
represented by a Matrix Product State (MPS) ansatz whose parameters grow linearly 
with the number of quantum particles. Future work includes explicit treatment of 
fermionic systems and numerical exploration of higher dimensional systems. 
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Inputs 




— > 


Output 





Figure 2.2.3: Inputs and outputs of the IE theory. See section [273] for the definition 
of dv^^. 



2.2.2 Inputs and Outputs of the Theory 



In general we consider Hamiltonians H = -f^odd ~'~ -^even' where the local terms 
that add up to H^^^ (or -f^gyen) fo^'^ ^ commuting subset. All the physically relevant 
quantities such as the lattice structure, A^, dimension of the spin d and the rank r 
are encoded in the eigenvalue densities. The output of the theory is a < p < 1 by 
which the IE distribution is obtained and du^^ serves as an approximation to the 
spectral measure. The inputs can succinctly be expressed in terms of the dimension 
of the quantum spins, and the nature of the lattice (Figure 2.2.3[ ) . 

2.2.3 More Than Four Moments of Accuracy? 



Alternatives to IE worth considering are 1) Pearson and 2) Gram-Charlier mo- 
ment fits. 



We illustrate in Figure [2 . 2 . 4| how the IE fit is better than expected when matching 
four moments. We used the first four moments to approximate the density using the 
Pearson fit as implemented in MATLAB and also the well-known Gram-Charlier 
fit [50]. In [l9j it was demonstrated that the statistical mechanics methods for 
obtaining the DOS, when applied to a finite dimensional vector space, lead to a 
Gaussian distribution in the lowest order. Further, they discovered that successive 
approximations lead naturally to the Gram-Charlier series [50]. Comparing these 
against the accuracy of IE leads us to view IE as more than a moment matching 
methodology. 



The departure theorem (Section 2.4.2) shows that in any of the higher moments 



(> 4) there are many terms in the quantum case that match IE exactly. Further, we 
conjecture that the effect of the remaining terms are generally less significant. 



41 



4-moment matching methods vs. IE: binomial local terms 



4-momen1 matching methods vs. IE: Wishart local terms 





• Exact Diagonalization 






— I.E. Theory: p = 0.046 






— Pearson 








— Gram-Charlier 





















• Exact Diagonalization 

— I.E. Theory: p = 0.42817 
-■-■Pearson 

— Gram-Charlier 




Eigenvalue 



Eigenvalue 



Figure 2.2.4: IE vs. Pearson and Gram-Charlier 

2.3 Spectra Sums in Terms of Probability Theory 

The density of eigenvalues may be thought of as a histogram. Formally for an 
m X m matrix M the eigenvalue distribution is [l5l p. 4]|l6l p. 101] 



^ m 

dvM{x) = - Ai (M)). 



(2.3.1) 



For a random matrix, there is the expected eigenvalue distribution [17], |1H1 P- 



362] 



-E 



X^5(^-A. (M)) 



i=l 



(2.3.2) 



which is typically a smooth curve p. 101] [46, p. 115]. 

The eigenvalue distributions above are measures on one variable. We will also 
need the multivariate measure on the m eigenvalues of M: 

dvM {x) = The symmetrized joint density of the eigenvalues. 

Given the densities for M and M', the question arises: What kind of summation 
of densities might represent the density for M + M'7 This question is unanswerable 
without further information. 

One might try to answer this using various assumptions based on probability 
theory. The first assumption is the familiar "classical" probability theory where 
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the distribution of the sum is obtained by convolution of the density of summands. 
Another assumption is the modern "free" probabihty theory; we introduce a finite 
version to obtain the "isotropic" theory. Our target problem of interest, the "quan- 
tum" problem, we will demonstrate, practically falls nicely in between the two. The 
"Slider" quantifies to what extent the quantum problem falls in between (Figure 



2.2.1 bottom). 



2.3.1 Classical 

Consider random diagonal matrices A and B of size m, the only randomness is in 
a uniform choice among the m! possible orders. Then there is no difference between 
the density of eigenvalue sums and the familiar convolution of densities of random 
variables, 

du'^ = dvA * dvB- (2.3.3) 

Comment: From this point of view, the diagonal elements of A, say, are identi- 
cally distributed random variables that need not be independent. Consider Wishart 
matrices |52] , where there are dependencies among the eigenvalues. To be precise let 
a G be a symmetric random variable, i.e.. Pa has the same distribution as a for 

(a, \ 

all permutation matrices P. We write, A = 



= diag(a). Similarly 



for B. 



Comment: The classical convolution appears in Figure 2.2.1 in two different ways. 
Firstly, in the definition of A (or B) , the eigenvalues of the odd (or even) terms are 
added classically. Secondly, A and B are added classically to form one end of the 
Slider. 

2.3.2 Free and Isotropic 

Free probability [I8t is recommended] provides a new natural mathematical "sum" 
of random variables. This sum is computed "free convolution" denoted 

du^ = duA ffl duB- (2.3.4) 
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Here we assume the random matrices A and B, representing the eigenvalues, have 
densities dvA and dvB- In the large m limit, we can compute the DOS oiA + Q^BQ, 



where Q is a /3— Haar distributed matrix (see Table 2.1). 

Comment: In this paper we will not explore the free approach strictly other than 
observing that it is the infinite limit of the isotropic approach (i.e., t — )■ oo in Eq. 



2.3.5). This infinite limit is independent of the choice of /3 (see Table 2.1). 





Real R 


Complex C 


Quaternions EI 


"Ghosts" 




1 


2 


4 


general (3 


Notation 


Q 


U 


S 




Haar matrices 


orthogonal 


unitary 


symplectic 


/3— orthogonal 



Table 2.1: Various /3— Haar matrices. 

We define an isotropic convolution. The isotropic sum depends on a copying 
parameter t and (3 (Table 2.1). The new Hamiltonian is the isotropic Hamiltonian 



ISO 



H,so = {A' ® It) + Q-^^ {It ® B') Qis, (2.3.5) 

where is a /3— Haar distributed matrix, A = A' ^It and B = % ® B' . For the 
copying parameter t = d, Hiso has the same dimension as H in Eq. 2.1.1 however, 
t > d allows us to formally treat problems of growing size. We can recover the free 
convolution by taking the limit: limj_>.oo du^^°^^'*'> = du^ . The effect oi is to spin 
the eigenvectors oiIt®B to point isotropically with respect to the eigenvectors of 
A. We denote the isotropic eigenvalue distribution by 



dyiso(M = dUA SisoiM dUB 

omitting t and /3 when it is clear from the context. 



(2.3.6) 



Comment: In Eq. 2.3.5, the and B init® B, can appear in any order. We 
chose this presentation in anticipation of the quantum problem. 

Comment: In this paper we primarily consider t to match the dimension of H. 
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2.3.3 Quantum 



Let du'^ denote the eigenvalue distribution for the Hamiltonian in Eq. 2.1.1 This 
is the distribution that we will approximate by du^^ . In connection to Figure 2.1.1 
the Hamiltonian can be written as 



^ = ^odd + ^even = Yl I ® ^M+i ® ^ + Yl I ® ^M+i ® I- (2-3.7) 



i=l,3,5,-' 



1=2,4,6,- 



We proceed to define a "quantum convolution" on the distributions dOj^ and dOs, 
which is /3-dependent 



du'i^^'^ = do A fflg dvB- 



(2.3.8) 



In general, without any connection to a Hamiltonian, let dvA and dvs be sym- 
metric measures on W^'^ . We define dz/'^'^^) to be the eigenvalue distribution of 



H = A + Q^'BQg 



where Qq 



gf ^ with 



(2.3.9) 



I, and Qf^ =ld 



N oM 



Q 



{A) 

q 



Id ® [(8)S'"' Qf^\ ® Id N even 

(2.3.10) 

and each q['^ is a /3— Haar measure orthogonal matrix of size d^ and 1^ is an 
identity matrix of size d. 

Comment: A, B and Qq are c?^ x d^ . 

Comment: In our examples given in this paper, we assume the local terms are 
independent and identically distributed (iid) random matrices, each of which has 
eigenvectors distributed with /3— Haar measure. 
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The tensor product in (2.3.10) succinctly summarizes the departure of the quan- 



tum case from a generic matrix as well as from the classical case. First of all the 
number of parameters in Qq grows linearly with whereas in Q it grows exponen- 
tially with A^. Second, the quantum case possesses isotropicity that makes it different 
from the classical, whose eigenvectors are a point on the orthogonal group (i.e., the 
identity matrix). 

Comment: General /3's can be treated formally [53]. In particular, for quantum 
mechanical problems /3 is taken to be 1 or 2 corresponding to real and complex entries 
in the local terms. /3 = 4 corresponds to quaternions. 

Definition. The Hadamard product of two matrices Mi and M2 of the same size, 
denoted by Mi o M2, is the product of the corresponding elements. 



Lemma 1. The elements of Qq defined in Eq. 2.3.10 are (dependent) random vari- 
ables with mean zero and variance . 



Proof. Here expectations are taken with respect to the random matrix Qq which is 
built from local Haar measure matrices by Eq. 2.3.10 The fact that E ( Qg'^'' 



Eg 



Orfiv follows from the Haar distribution of local terms. Thus E {Qq) = 



by independence of Qq^^ and Q\^' . Further, each element in Qq involves a dot product 



between columns of Q^j^'' and Qq^' ■ In every given column of Qq'^' any nonzero entry is 



(A) 



a distinct product of entries of local Q's (see Eq,2.3.10). For example the expectation 
value of the 1, 1 entry is E (^qi^\j^qi^\j^^^ = E (^qi^^j^^^^ E (qi^^qj^^^- Because 
of the Haar measure of the local terms, this expectation is zero unless i = j. We 
then have that 
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® Id) fid ® 



9 



(g)fifrf-v,.) ri, 



d2 



A^odd 
even 



(2.3.11) 



where Ji = i x i matrix of all ones. We use facts such as {Ji/if' = {Ji/i), {Ji/i) ® 
(Ji/i) = (Jp/i^) and the variance of the elements of an i x z /3— Haar matrix is 

□ 



2.4 Theory of Isotropic Entanglement 

2.4.1 Isotropic Entanglement as the Combination of Classi- 
cal and Isotropic 

We create a "Slider" based on the fourth moment. The moment of a random 
variable defined in terms of its density is = J x^dvM- For the eigenvalues of an 
m X m random matrix, this is — ETrM'^. In general, the moments of the classical 
sum and the free sum are different, but the first three moments, mi, m2, and 
are theoretically equal [W, P- 191]. Further, to anticipate our target problem, the 
first three moments of the quantum eigenvalues are also equal to that of the iso and 
the classical (The Departure and the Three Moments Matching theorems in Section 



2.4.2). These moments are usually encoded as the mean, variance, and skewness. 

We propose to use the fourth moment (or the excess kurtosis) to choose a correct 
p from a sliding hybrid sum: 



rfz/" ^ rfz/'^ = pdu" + (1 - p)diy''° (2.4.1) 

Therefore, we find p that expresses the kurtosis of the quantum sum (7^) in terms 
of the kurtoses of the classical (72) and isotropic (72^*°) sums: 
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il = in2 + {i-p)ir^ p=^:^—^o- (2-4.2) 

12 12 

Recall that the kurtosis 72 = where is the variance. Hence kurtosis is the 
correct statistical quantity that encodes the fourth moments: 



ml = ^ETr {A + U^BuY , mf = ^ETr {A + Q^BQ)\ ml = ^ETr {A + Q^BQ^Y 

(2.4.3) 

where 11 is a random uniformly distributed permutation matrix, Q is a /3— Haar 
measure orthogonal matrix of size , and Qg is given by Eq. 2.3.10 Unless stated 



otherwise, in the following the expectation values are taken with respect to random 
eigenvalues A and B and eigenvectors. The expectation values over the eigenvectors 
are taken with respect to random permutation U, [5—Haar Q or Qq matrices for 
classical, isotropic or quantum cases respectively. 

2.4.2 The Departure and The Matching Three Moments 
Theorems 

In general we have the i^^ moments: 



mr = -W£i{A + Q^BQY 
ml = -ETi (A + Q^BQgY , and 
m'r = —ETi (A + U^BUy . 

where m = d^ . If we expand the moments above we find some terms can be put 
in the form ETr (^A^^ B"^^ Q and the remaining terms can be put in the form 
ETr |. . . Q^B-^Q,A-^QJB-^Q, . . .} . The former terms we denote non- departing; 
the remaining terms we denote departing. 

For example, when i = 4, 
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rriA 



m 
1 

m 



ElTr 



+ iA^Q^BQ + iA^Q^B^Q + AAQ^B^Q + 2 (AQ'^BQ)^ + B^ (^.4.4) 
A^ + AA^Q^BQ^ + AA^Q^B^Q^ + AAQIB^'Q^ + 2 (AQjBQq)^ + B' 



-^EjXr A^ + AA^lfBIi + AA^lfB'^Ii + AAlfB^Ii + 2{AIi^BIif + B^ | 



where the only departing terms and the corresponding classical term are shown as 
underlined and bold faced. 

Theorem. (The Departure Theorem) The moments of the quantum, isotropic 

and classical sums differ only in the departing terms: ETr |. . . B-^Q,A-^Ql B-^Q, . . .} . 



Proof. Below the repeated indices are summed over. If A and B are any diagonal 
matrices, and is Q or Qq or 11 of size mxm then E = 1/m , by symmetry and 
by Lemmajljrespectively. Since the ETr i^AQ^BQ,^ = E (^qfjUibj), where expectation 
is taken over randomly ordered eigenvalues and eigenvectors; the expected value is 



(^) E {a,bj) for any i or j. Hence, ^ETr {AQlBQ,) = E (aibj) = E (a^) E (6 



which is equal to the classical value. The first equality is implied by permutation in- 
variance of entries in A and B and the second equality follows from the independence 
of A and B. 

□ 



Therefore, the three cases differ o nly in the terms ^ETr [AQ'^BQ] , |ETr {AQ^BQg)' 
and ^ETr {AU'^BuY in Eq. 



2.4.4 



Theorem. (The Matching Three Moments Theorem) The first three moments 
of the quantum, iso and classical sums are equal. 
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number of 
sites 


number of odds 
or evens (A^ odd) 


number of odds 
or evens {N even) 


size of 
H 


Number of 
copies 


Dimension of 
quidits 




K — 2 


k AA - ^ k - ^-'^ 
odd 2 ' even 2 


m = 


t 


d 



size of local terms 


jth j^oment 


Z*'^ cumulant 


mean 


variance 


skcwness 


kurtosis 


n = d^ 


nil 


1^1 






7i 


72 



Table 2.2: Notation 



Proof. The first three moments are 

mj'^ = ^ETr {A + B) 
m^*) = i^ETr {A + QlBQ.f = i^ETr {A^ + 2AQlBQ, + B^) 
mi'^ = ^ETr {A + Q^BQ,)' = ^ETr {A^ + SA^Q^BQ, + 3AQlB^Q, + B^) , 

(2.4.5) 

where Q, is Q and Qq for the iso and the quantum sums respectively and we used 
the familiar trace property Tr(MiM2) = Tr(M2Mi). The equality of the first three 
moments of the iso and quantum with the classical follows from The Departure 
Theorem. 

□ 



Furthermore, in the expansion of any of the moments > 4 all the non-departing 
terms are exactly captured by IE. These terms are equal to the corresponding terms 
in the classical and the isotropic and therefore equal to any linear combination of 
them. The departing terms in higher moments (i.e.,> 4) that are approximated by 
IE, we conjecture are of little relevance. For example, the fifth moment has only two 
terms (shown in bold) in its expansion that are departing: 



ms = ;^ETr (^A^ + bA'^QlBQ, + SA^QjB^Q, + SA^QlB^Q, + 5A (AQ^BQ.)' 
5 (AQ^BQ.)^ Q^BQ. + SAQlB^Q. + fi^' 



(2.4.6) 



By the Departure Theorem the numerator in Eq. 2.4.2 becomes. 
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72 - l2 



Q _ ISO 



E<Tt 



{aq^bq,Y-{aq^bqY]] 



m 



(2.4.7) 



and the denominator in Eq. 2.4.2 becomes 



72 - 72 



[AYl^ Blif - [AQ^ BQY 



m 



where as before, Q is a /3— Haar measure orthogonal matrix of size m 

T 



Qq^] Qq^^ given by Eq. 



evaluation of p reduces to t 
12X81 



2.3.10 



tie eva. 



(2.4.8) 



and denote the fourth cumulants. Therefore, 



nation of the right hand sides of Eqs. 2.4.7 and 



Below we do not want to restrict ourselves to only chains with odd number of 
sites and we need to take into account the multiplicity of the eigenvalues as a result 
of taking the tensor product with identity. It is convenient to denote the size of the 



matrices involved hj m = = tn^, where n = d"^ and k 



N-l 



and t is the number 



of copies (Section 2.3.2 and Table 2.2). 



2.4.3 Distribution of A and B 

The goal of this section is to express the moments of the entries of A and B (e.g., 
and m'^-^) in terms of the moments of the local terms (e.g for odd local terms 
m^'^^ , rui^'^) . Note that A and B are independent. The odd summands that make 
up A all commute and therefore can be locally diagonalized to give the diagonal 
matrix A (similarly for B), 
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iV-2 

A = I (g) (g) I 

i=l,3,-- 
N-1 

B = ^ I (g) Aj (g) I, 

i=2,4,-- 



(2.4.9) 



where Aj are of size (P and are the diagonal matrices of the local eigenvalues. 

The diagonal matrices A and B are formed by a direct sum of the local eigenvalues 
of odds and evens respectively. For open boundary conditions (OBC) each entry has 



a multiplicity given by Table 2.3 



OBC 


N odd 


N even 


A 


d 


1 


B 


d 


Ci2 



Table 2.3: The multiplicity of terms in A and B for OBC. For closed boundary 
conditions there is no repetition. 



Comment: We emphasize that A and B are independent of the eigenvector struc- 
tures. In particular, A and B are the same among the three cases of isotropic, 
quantum and classical. 

We calculate the moments of A and B. Let us treat the second moment of A [B 
is done the same way). By the permutation invariance of entries in A 



m^^E{al) = E (aSJ) + ■ ■ ■ + Alf ) 

= E [A; (A^) + A; (A; - 1) A^^U^^)] (2.4.10) 
= fcm°^d + A; (A; - 1) m??^ 

where expectation is taken over randomly chosen local eigenvalues, m2^^ = E (A?) 
and m^^*^ = E (AjAj) for some uniformly chosen i and j with i ^ j. The permutation 
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invariance assumption implies E (a^) = E (a^) for alH = 1 ■ ■ ■ m. 

Comment: The key to this argument giving is that the indices are not sensitive 
to the copying that results from the tensor product with I^; at the boundaries. 

Next we calculate the correlation between two diagonal terms, namely rrii^ = 
E (ajOj) for i ^ j. We need to incorporate the multiplicity, denoted by t, due to the 
tensor product with an identity matrix at the end of the chain, 



m{m — \) 



E 



n \ " n „ 



ii,--- ,ifc=l 



A; (A; - 1) E (A)^ + k {prob (A^) E (A^) + prob (A1A2) E (A1A2)} 
k{k-l) m^dd + I (t^fc-i _ 1) ^odd ^ ^^^k-i _ ^oddj 



where, prob (A^ 



and prob (A1A2) 



m—1 



Similarly for B. 



2.4.4 Evaluation and Universality of p 



75-72' 



Recall the definition of p; from Eqs. |2.4.2[ |2.4.7| and 2.4.8 we have, 

ETr {Ali^BIif - ETr {AQlBQ^f 



1 — p 



ETr {ART Buy - ETr (AQ^BQ) 



2 • 



The classical case 



(2.4.12) 



IeIY {AlfBuY = l^E±a^b^ = E (a^) E (6?) 

i=l 



A B 
1712 "^2 • 



(2.4.13) 



Comment: Strictly speaking after the first equality we must have used 6,rj instead 
of bi but we simplified the notation as they are the same in an expectation sense. 



The general form for the denominator of Eq. 2.4.12 is 



-ETr 



m 



{AU^BuY-{AQ^BQy 



m 



E {afbf - ttittkbjbp (qjiqjkqpkqpi)} ■ (2.4.14) 
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moments 


expectation values 


count 


ml 


^ ) — m(m/3+2) 


2 

m 


m2mu 


E (|gi,igi,2| ) - „(^;3+2) 


2m^ (m — 1) 






(m — 1)^ 




E {113(124) — n(n/3+2)(n-l) 





Table 2.4: The expectation values and counts of colliding terms in Q when it is 
/3— Haar distributed. In this section we use the first row; we include the last three 
rows for the calculations in the appendix. 



It's worth noting that the arguments leading to Eq. 2.4.18 hold even if one fixes A 
and B and takes expectation values over 11 and a permutation invariant Q whose 
entries have the same expectation value. The right hand side of Eq. 2.4.18 is a 



homogeneous polynomial of order two in the entries of A and B; consequently it 
necessarily has the form 

1 



m 



-ETr {AU^BU) ^ - [AQ^ BQ) ^ =ci{B,Q)m^ + ca (5, Q) m^^ 

but Eq. 2.4.14 must be zero for A = I, for which = m^i = 1. This implies that 
Ci = — C2, allowing us to factor out (m^ — m^^). Similarly, the homogeneity and 
permutation invariance of B implies, 



-ETr 



m 



{au^bu)^-{aq^bqY = { 



[mo — m 



The right hand side should be zero for B = I, whereby we can factor out ( 



mo — m 



^ ) 



-ETr 



m 



{AU^BUY - [AQ'^BQY = [mi - m^^ (mf - mfj / (Q) , (2.4.15) 



where m^ 



Efa? 



m 



B 



E (b'j), and m^^ = E (0^,0^) 



mfi 



E (6j , 6j ) . Moreover 



/ (Q) is a homogeneous function of order four in the entries of Q. To evaluate / (Q), 
it suffices to let A and B be projectors of rank one where A would have only one 
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nonzero entry on the i position on its diagonal and B only one nonzero entry 
th 

on the j position on its diagonal. Further take those nonzero entries to be ones, 
giving m^i 



— mfi = and = mf = 1/m, 



-ETr 



m 



But the left hand side is 



-ETr 



m 



{AU^BUY - {AQ^BQY 



-J{Q) 



(2.4.16) 



m 



n ^ — ^ ^ — ^ 



m 
1 

m I m 



- E (4) 



where, we used the homogeneity of Q. Consequently, by equating this to / (Q) /m^, 
we get the desired quantity/ (Q) = |l — mE (g^^) } . 
Our final result Eq. |2.4.15| now reads 



-ETr 



{AU^BUY -{AQ^BQY = ( 



rrir, — m 



rrin — m; 



^,){l-mE (gj)}. 

(2.4.17) 

The same calculation where each of the terms is obtained separately yields the same 
result (Appendix). In this paper p is formed by taking Q to have a /3— Haar measure. 



Expectation values of the entries of Q are listed in the Table 2.4 



We wish to express everything in terms of the local terms; using Eqs. 2.4.10 and 
2.4.11 as well as tn^ = m, 



^2 - ^ti 



tk {n — 1) n^^^ 

m — 1 
tk {n — 1) n^~^ 

m — 1 



m. 



odd _ ™,odd 
'"'1,1 



[m 



even 
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giving 



-E 



m 



Tr {AU^BUf - Tr {AQ^BQ)' 



m, 



)dd _ _odd^ /_even _ _even\ 
: '"'1,1 i y'h '"'1,1 ) 



X 



km {n — r 
nim — 1) 



{1-mE (g^} (2.4.18) 



We now proceed to the quantum case where we need to evaluate 
1 



-E 



m 



{AlfBYiY-Ti{AQlBQqy 



In this case, we cannot directly use the techniques that we used to get Eq. |2.4.18 
because Qq is not permutation invariant despite local eigenvectors being so. Before 
proceeding further we like to prove a useful lemma (Lemma [2]). Let us simplify 
the notation and denote the local terms that are drawn randomly from a known 
distribution by Hi i^i = H^'-'> whose eigenvalues are A; as discussed above. 

Recall that A represents the sum of all the odds and Q~^BQg the sum of all the 
evens, 



H, 



odd 



E 

Z=l,3,5/ 



I ® H'-^^ ® I, and H, 



even 



J2 l^H^^^^l 



'=2,4,6/ 



Hence, the expansion of ^E Tr (^AQ^ BQgj 



amounts to picking an odd term, an 
even term, then another odd term and another even term, multiplying them together 



and taking the expectation value of the trace of the product (Figure 2.4.1 ). Therefore, 
each term in the expansion can have four, three or two different local terms, whose 
expectation values along with their counts are needed. These expectations are taken 
with respect to the local terms (dense rf^ x random matrices). 

The expectation values depend on the type of random matrix distribution from 
which the local terms are drawn. The counting however, depends on the configuration 
of the lattice only. We show the counting of the number of terms, taking care of the 
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A 
Q'BQ 

A 
Q'BQ 



Tr{AQ^BQ,y 



can be visualized 



Figure 2.4.1: The terms in the expansion of 
as picking an element from each row from top to 
has k of the local terms corresponding to a chain with odd number of terms. Among 
fc^ terms roughly fc^ of them differ among the classical, isotropic and quantum cases 



Dottom and multiplying. Each row 



(See Eqs. 2.4.19 and 2.4.20). An example of such a choice is shown by diamonds. 



boundary terms for an open chain, along with the type of expectation values by 
which they need to be weighted: 

For odd (fc odd terms and k even terms) 



Four i/W's ■.k'^{k- if d^-'^Wi [H^^^Y , ui e {5, ■ ■ ■ , 8} 
Three ifW's : 2k^k - 1) =^ rf^-^^ETr ^) ETr {H^^^Y , U2 E {4,5,6} 

Two H^-^'s: {k - if Not Entangled d^-^ |ETr ([H^^^]^^ }^ 

Two i/(-)'s : {2k - 1) Entangled ^ d^--'WIi [(if(') ® l) (l ® ® l) (l 

(2.4.19) 

For A^ even {k odd terms and k — 1 even terms) 



57 



Four //W's -.kik- if {k-2)^ d^-"iETr (i/W)^ , e {5, • • • , 8} 
Three //W's : A; (A; - 1) (2A; - 3) ^ d^-^^ETr ([H^^^ ETr (//(O)^ , G {4, 5, 6} 

Two i/W's : (A; - 1) (A; - 2) Not Entangled ^ d^-^ jlTr ([i^^'^^) }^ 

Two i/(-)'s : 2 (A; - 1) Entangled ^ ci^-^ETr [(i/C) O l) (l ® ® l) (l ® 

(2.4.20) 

Here i^i and U2 indicate the number of sites that the local terms act on (i.e., oc- 
cupy). Therefore, ^E Tr (AQ^-BQg) ^ is obtained by multiplying each type of 
terms, weighted by the counts and summing. For example for ui — h and U2 — 3, 
when N is odd. 



-E 



IV [AQlBQ.Y] = ^ [d^-^k^ (k - 1)^EIV 



2^ {k - 1) rf^-^ETr ( [i/W] 2 j (//«)' + (k - if d^-^ |ETr ( [ijW] ^ j | + 



(2A; - 1) d^-^ETr ® l) (l ® i/('+i)) (ifC) ® l) (l ® i/^+i))] } 



(2.4.21; 



and similarly for N even. 



^E 



"Tt {AQ^BQ.Y] = (A; - 1) (A; - 2) d^-^ETr (//O)' + 

A; (2A; - 3) cI^-^ETr ETr {H^'^f + {k-2) rf^-^ |ETr ( }' + 

2d^-3ETr [(i/C) ® I) (I ® i/C+i)) ® I) (I ® //('+!))] } . 

(2.4.22) 

The expectation values depend on the type of random matrix distribution from 
which the local terms are drawn. We will give explicit examples in the following 
sections. In the following lemma, we use E (i/C)) ^ijld2 and E (/f^)) =m2ld2. 



Lemma 2. In calculating the KTr (^AQ^ BQgf if at least one of the odds (evens) 
commutes with one of the evens ( odds ) then the expectation value is the same as the 
classical expectation value. Further if the local terms have permutation invariance of 
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eigenvalues then the only quantum expectation value that differs from classical is of 
Type II (see the proof and the diamonds in figure 2.4- !)■ 



Proof. This can be shown using the trace property Tr (MP) = Tr (PM). In calcu- 
lating ETr (fjodd j^even j^odd j^even^ . j£ ^^^^ q£ ^^^q odd (cven) temis commutes with any 

of the even (odd) terms to its left or right then they can be swapped. For example 
one gets ETr {Hf'^H^'"''^Hl'"'''Hf'^) = ETr {Hf'^Hf'^H^'""^Hl''^'') which is just the 
classical value. Hence the only types of expectations that we need to worry about 
are 



Type I 



and 



Type II 



now we show that with permutation invariance of the local eigenvalues the first type 



are also classical leaving us with the "diamond terms" alone (Fig. 2.4.1). Consider 
a Type I term, which involves three independent local terms, 



;i;ETr[(l^2®/i'(3'®IdJV-4)(nCE3H(2)®I^^_3)(l^2Cg>H(3)®I^jV-4)(ld3®H(4)®I^iV_5)] 



/i2m2. 



This follows immediately from the independence of H^^^ , which allows us to take its 
expectation value separately giving a fi and leaving us with 



m 



ETr 



(1^2 (g) if O (I (g) if (g) 



Therefore the only relevant terms, shown by diamonds in Fig. 2.4.1 are of Type 
II. As an example of such terms consider (here on repeated indices are summed over) 



59 



-ETr 



(P I V n«2jli2 «lP2jlfc2 / \ i2«3,fc2fc3 «2«3,P2fc3 



(2.4.23) 



where the indices with subscript 2 prevent us from treating the two expectation 
values independently: if*-^-' and H'^'^^ overlap at the second site. The number of such 



terms is 2k — 1, where k = 



□ 



Therefore, we have found a further reduction of the terms from the departure 
theorem, that distinguishes the quantum problem from the other two. Luckily and 
interestingly the kurtosis of the quantum case lies in between the classical and the 
iso. We emphasize that the only inputs to the theory are the geometry of the lattice 
(e.g., the number of summands and the inter- connectivity of the local terms) and the 
moments that characterizes the type of the local terms. 

Comment: The most general treatment would consider Type I terms as well, 
i.e., there is no assumption of permutation invariance of the eigenvalues of the local 
terms. This allows one to treat all types of local terms. Here we are confining to 
random local interactions, where the local eigenvectors are generic or the eigenvalues 
locally are permutation invariant in the expectation value sense. 

The goal is to find p by matching fourth moments 

ETr {AW^BUY - ETr {AQ'^BQgY 
ETr (AUTBU)^ - ETr [AQ^BQf 



for which we calculated the denominator resulting in Eq. 2.4.18, where E [\qij\ ) = 
m(m/3+2) ^'^^ /3— Haar Q (Table 2.4). If the numerator allows a factorization of the 



moments of the local terms as in Eq. 2.4.18 then the value of p will be independent 



of the covariance matrix (i.e., eigenvalues of the local terms). 



Lemma. (Universality) p p {N,d, (3) , namely, it is independent of the distribu- 
tion of the local terms. 
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Proof. We use a similar techniques as we did in the isotropic case. The general form 
for the numerator of Eq. 2.4.12 is (denoting Lemma |2] by L3) 



-ETr 



m 



{AU^BU)'-{AQ^BQ,y 



L3 (2A;-1) 



rf3 



ETr { {QT'AiQi ® (Id 8) Qr+Wig^+i)' 



(2fc - 1) 

- [{QT'AiQi^h) (ld8)Qr+Wigm)]'} (2.4.24) 

where the expectation on the right hand side is taken with respect to the local terms 
and if ('+^) . The right hand side is a homogeneous polynomial of order two in 



the entries of A;, as well as, in the entries of A;+i; consequently Eq. 2.4.24 necessarily 
has the form 



Ci (A^^^^, Qq(1(1, Qeven) m^^^ + C2 (ii^^^^, Qq^^, Qeven) 



odd 



but Eq. 



2.4.24 



must be zero for A; = J, for which m^^^ = m^^*^ = 1. This implies 
that Ci = —C2- By permutation invariance of the local terms we can factor out 
^^odd _ ^odd~j Similarly, the homogeneity and permutation invariance of ii(^+^) 
implies. 



m 



odd ^odd^ (g^^^^ Q^^^^) ^even ^ (^^^^^ q^^^^^ ^evem 



The right hand side should be zero for A/+i = /, whereby we can factor out ( 



m 



J • 



even 



m 



even\ . 
1,1 J' 



hence the right hand side of Eq. 2.4.24 becomes 



f^odd ^odd\ / even 



rf3 



m 



m 



/,(Qodd. Qeven) (2.4.25) 



where fg (Qq^j^, Qeven) is a homogeneous function of order four in the entries of 
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^odd ^^^^ Qeven- To evaluate /g, it suffices to let A/ and A^+i be projectors 

of rank one where A; would have only one nonzero entry on the i position on its 

th 

diagonal and A;+i only one nonzero entry on the j position on its diagonal. Further 

1/n. 



B 
1,1 



and 



B 



take those nonzero entries to be ones, giving m^^ = m 
Using this choice of local terms the right hand side of Eq. |2.4.24| now reads 



(2fc-i: 



ETr 



(2.4.26) 



where here the expectation value is taken with respect to random choices of local 



eigenvectors. Equating this and Eq. 2.4.25 



/.(Qodd'^even) = n^ETr | | ® /,) (/, ® |gf ^)|)(2.4.27) 



To simplify notation let us expand these vectors in the computational basis |gf ^) = 
Wiij2Ki)N2) and \q^-^^'') = Vi^i^\i2)\i2) ■ The ffist term on the right hand side of Eq. 



2.4.26 the classical term, is obtained by assuming commutativity and using the 



projector properties. 
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Tr 



(kf>tei"i sir "xr"!)^ 



TV 



><r"l) 



Tr (w''")] = Tr uv{uv) 



t 



I l|2 
|Mf lip 



E 



(2.4.2J 



where ||.||p denotes the Frobenius norm and cXj are the singular values of uv. The 
second term, the quantum term, is 



Tr 



Tr[Ui^i^Uj^j2Vj2i^Vk2k3Ujik2Umim2Vm2k3Vi2i3\hh^^ = 

^ ' 32K2 \ / m2K2 \ /m2«2 ^ '10 



{u^uvv'^\ . {u^uvv'^\ . — Tr 



m2j2 



UV {uv)^ 
UV (uv)^ 



J2t2 
2 



where we used the symmetry of (^uv {uv^^ = uv (uv)^ uv {uv^ 
Now we can calculate 



(2.4.29) 



i=l 



fi (^odd' <5even) = | Hii-yllF - uv {uv)^ 



giving us the desired result 



(2.4.30) 
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-ETr 



m 



X 



d{2k-l) (m°dd-mOfd) ( 

uv {uv) 



^even _ ^even^ 



E 



I l|2 

\uv\\y - 



(2.4.31; 



from which 



1 — p 



ETr {AH^BUY - ETr [AQ-^BQ^f 
ETr (An^^n)^ - EIV {AQ-^BQf 



(2.4.32) 



d{2k - 1)E 



I l|2 



fcm(n— 1) 
n(m— 1) 



{1-mE (g^.)} 



The dependence on the covariance matrix has cancelled- a covariance matrix is one 
whose element in the i, j position is the covariance between the i^^ and j*^ eigenvalue. 
This shows that p is independent of eigenvalues of the local terms which proves the 
universality lemma. 

Comment: To get the numerator we used permutation invariance of A and B and 
local terms, to get the denominator we used permutation invariance of Q. 

□ 



Comment: It is interesting that the amount of mixture of the two extremes needed 
to capture the quantum spectrum is independent of the actual types of local terms. 
It only depends on the physical parameters of the lattice. 



2.4.5 The Slider Theorem and a Summary 

In this section we make explicit use of /3— Haar properties of Q and local terms. 



To prove that there exists a < p < 1 such that the combination in Eq. 2.4.2 



IS 



convex we need to evaluate the expected Frobenius norms in Eq. 2.4.31 
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Lemma 3. E||'ut'||^ = \/d and K 
terms have [5—Haar eigenvectors. 



uv [uv] 



/32[3d(d-l)+l]+2/3(3d-l)+4 , , , 
d{l3d^+2f ' ^^^^ 



Proof. It is a fact that G = uXjsd?-, when u is uniform on a sphere, G is a. d ^ d 
/3— Gaussian matrix whose expected Frobenius norm has a distribution denoted 
here by xpiP (similarly for v). Recall that E {xV) — ^ ^ ixt) — h{h + 2). 



¥.\\uv\\\ E(x/3d2)'= E||(GiG2)||f 



(2.4.33) 



E \\uv\ 



1 2 (P — ^ '-^2 



^E II Gi Go II T? — o 

(/3rf2)2 II '11^ (/3rf2)2 
2^(^) ~ - 



fc=l 



(/3ci2) 



The quantum case, v)u 



G{Gi 



W\ and are Wishart matrices. 



wn l|2 ; 
ll<-i||p 



similarly v'^v = = T^hr, where 

II Gr2 II / --^ II 



IIG2II 



E 



uv [uv 



ETr {W1W2) 



ETr {WiW2f 
[d^(3 {d^(3 + 2)y 



(2.4.34) 



hence the complexity of the problem is reduced to finding the expectation of the 
trace of a product of Wishart matrices. 



i<ijki<d I yj^^ ^kVi 

(2.4.35) 

where 11 denotes the product of the elements of the matrix. There are three types 



of expectations summarized in Table 2^ 
In Table 
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Notation 


Type 


Count 


X 




{d - If 


Y 


i^k & j^l 

or 

i^k k 3=1 


2d^ {d - 1) 


Z 


i = k & j = I 


d^ 



Table 2.5: Expectation values. 



X = 
Y = 

Z = 



E [n {xiXk) {yiVi)] 



We now evaluate these expectation values. We have 



X = U 



( X0d \^ 

Xl3(d^l) 











V 



/ 9P 9p 

90 90 
DC DC 



V 



by QR decomposition, where and Xh denote an element with a ^—Gaussian and 
Xh distribution respectively; DC means "Don't Care" . Consequently 



X = 



= n 




9^ 

X/3{d-l) 

aXpd 90a + X0{d-i)C 
bX0d 9i3b + Xi0{d-i)d 
Xl,a'b'9} = P'd. 




n 



aX0d 90a 
bX0d 9 (0b 
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where we denoted the four independent Gaussian entries by a, b, c, d to not confuse 
them as one number. From Eq. |2.4.35 we have 



Y 



E 



E 



Z = E{x^yY = E{xU)^{xt)=PdWd + 2)f3{P + 2). 



Eq. 2.4.34 now reads 



E 



UV [UV] 



2 _ /32 [?,d {d-l) + l]+ 2/3 {U - 1) + 4 
F ~ d{pd^ + 2f ■ 



(2.4.36) 

□ 



Theorem. (The Slider Theorem) The quantum kurtosis lies in between the clas- 
sical and the iso kurtoses, '^2° ^ I2 — 12- Therefore there exists a < p < 1 such 
that 7I = + {I — p) 72^^°. Further, hmTv-s-oo p = 1- 



Proof. We have |l - ^J2Tj=i1ij} ^ 0' since < = ^- The last in- 

equality follows from qf- < 1 . Therefore, Eq. 



2.4.8 



72^° - 72 



m 



odd 



km [n — 1) 
n{m — 1) 



m 
2 



even _ ^even^ ^ 



{mE (gfi) - 1} < 0. 



From Eqs. |2.4.29| and |2.4.28| and using the fact that the singular values o"j < 1 
we have 



UV [UV] 



\uv\\-p 



■t=l i=l 

which proves 7I — 72 < 0. In order to establish 72''° < 72 ^ 72? we need to show that 



72 - 7I < 72 - 72"°- Eq. 2.4.32 after substituting mE (qfj) = j^^^ from Table 



2.4 
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and Eqs. |2.4.33[ [2X36] reads 



P 



(1 - rf-^'^-i) 



k-1 

k 



/3{d^ + -2d + l) +4d-2' 
{d - 1) i(3d^ + 2) _ 



d 



-2/fc+l 



l + /3dy2 J \d+l 



d 



(2.4.37) 



We want to show that < 1 — p < 1 for any integer k > 1, d > 2 and (3 > 1. 
All the factors are manifestly > 0, therefore 1 — p > 0. The first two factors are 
clearly < 1 so we need to prove that the term in the braces is too. Further, k = 1 
provides an upper bound as (l — ^^^rij^Ti^r] < (l — Tvi^crTTrV We rewrite the term 
in the braces 



d {(3d^ + 2) [P + -2d+l) + Ad-2] 
iPd^ + 2f {d + if {d - 1) ' 

but we can subtract the denominator from the numerator to get 

iPd + 2) [/3 (d^ - 2d^) +2{d^ -d^ -1)]>0 Wd>2. 



(2.4.38) 



This proves that (2.4.38) is less than one. Therefore, the term in the braces is 



less than one and hence < p < 1. Let us note the following limits of interest (recall 
N~l = 2k) lim(l-p) = '^hjll'^^^i 

d^oo k 



lim (1 

iV-5-oo 



P) 



N 







the first limit tells us that if we consider having two local terms and take the local 
dimension to infinity we have essentially free probability theory as expected. The 
second limit shows that in the thermodynamical limit (i.e., N — )■ oo) the convex 
combination slowly approaches the classical end. In the limit where /3 — )■ oo the (3 
dependence in (1 — p) cancels out. This is a reconfirmation of the fact that in free 
probability theory, for /3 — )■ oo, the result should be independent of (3. We see that 
the bounds are tight. 
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□ 



Comment: Entanglement shows itself starting at the fourth moment; further, in 
the expansion of the fourth moments only the terms that involve a pair of local 
terms sharing a site differ. Note that when the QMBS possesses a translational 
symmetry, there is an additional complication introduced by the dependence of the 
local terms. Though, in this case, the non-iid nature of the local terms complicates 
the matter theoretically, we have not seen a practical limitation of IE in our numerical 
experiments. 

Comment: One could from the beginning use free approximation instead of 
isotropic (m — )■ oo), in which case the proofs are simplified. 

We now summarize the main thesis of this work. We are interested in the eigen- 
value distribution of 



^ = ^odd + ^even = I ® ® I + Yl ^ ® ^/,m ® I, 

1=1,3,5,- l=2,i,6,- 
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which in a basis in that H^^^ is diagonal reads H — A + Q^^BQq. Since this 
problem has little hope in being solved exactly we consider along with it two known 
approximations : 



H = A + Q-'BQ, 
Hiso = A + Q-'BQ. 

We proved that the first three moments of the three foregoing equations are 
equal. We then calculated their fourth moments as encoded by their kurtoses (72's) 
analytically and proved that there exists a < p <l such that 

7|=P72 + (l-p)7r- 

It turned out that the only terms in the expansion of the fourth moments that were 
relevant were 



ETr \^(AU-'BUf - {AQ-^BQ^^j 



(2.4.39) 



ETr {{AYl-^BHf - {AQ-^BQf} ' 

Through direct calculation we found that the numerator ETr | {AU'^BU)^ — (^AQ~^BQqj ^ | 
evaluates to be 



d {2k - 1) (mOdd - mO?d) (mf^^ - mf^^) E (^\\uvfp - 
and the denominator ETr l^{AU-^BUf - {AQ-^BQf^ 



uv (uv) 



2 



m 



odd _ ^odd^ ^^even _ ^eyen) ] {l - mE (4)] . 



odd\ / 
1,1 V 



J 
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Therefore 1 — p does not depend on the local distribution and can generally be 
expressed as 



d{2k - 1)E \\uv\\y 



1 — p 



uv (UV) 



If we further assume that the local eigenvectors are /3 — Haar distributed we get 



p = (1 - [i - (V 

I3{d^ + -2d + l)+Ad-2 
{d - 1) (/3c/2 + 2) 



1 -d 



-2k+l 



1 + /3rfV2 



d 



d+l 



X 



Next we asserted that this p can be used to approximate the distribution 



diy'^ ^ du''' = pdu' + {l-p) du''°. 

We argued that the spectra obtained using Isotropic Entanglement (IE) are accurate 
well beyond four moments. 

For illustration, we apply IE theory in full detail to a chain with Wishart matrices 
as local terms. Other types of local terms (e.g. GOE, random ±1 eigenvalues) can 



be treated similarly; therefore in Section 2.6 we show the plots comparing IE with 
exact diagonalization for these cases. 



2.5 A Detailed Example: Wishart Matrices as Lo- 
cal Terms 

As an example take a chain with odd number of sites and for the local terms in 



Eq. 2.1.1 pick i/C) = W^W, where W is a rank r matrix whose elements are picked 
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randomly from a Gaussian distribution (/3 = 1); these matrices W'^W are known as 
Wishart matrices. Clearly the maximum possible rank is r = (P for each of the local 
terms. 

Any cumulant is equal to the corresponding cumulant of one local term, denoted 



by ft, times the number of summands in Eq. 2.1.1 In particular, the fourth cumulant 
of H is ^4^"^'* = (A^ — 1) K4. Below we drop the superscripts when the quantity 
pertains to the whole chain. Next we recall the definitions in terms of cumulants of 
the mean (/x), the variance (cr^), the skewness (71), and the kurtosis (72) 



^ = or2 = «:2 71 = ^ ^2 = ^ = ^_3 . (2.5.1) 

2.5.1 Evaluation of p 



72-72 

The moments of the local terms are obtained from MOPS 



mi = f3r 
m2 = Pr[l3 {r + n - 1) + 2] 
mg = /3r {(3^ [n^ + (r - 1) (3n + r - 2)] + 6/3 (n + r - 1) + 8} 
7714 = Pr {48 + [n^ + 6n'^{r - 1) + n (6r - 11) (r - 1) - 6 (r^ + 1) + + llr] 
+2/32 [6 + r^) + 17 (n(r - 1) - r) + 11] + 44/3 + r - 1)} 

mi l = P'^f' {f' ~ 1) 



(2.5.2) 



which for real matrices (3 = 1 yields 



mi = r 
m2 = r {r + n + 1) 
ms = r (n^ + 3n + Srri + 3r + + 4) (2.5.3) 
m4 = r (672^ + 21n + Qrn^ + 17rn + 21r + Qnr'^ + 6r'^ + + + 20) 

mi^i = r (r — 1) . 
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The mean, variance, skewness, and kurtosis are obtained from the foregoing rela- 
tions, through the cumulants Eq. |2.5.1[ We drop the superscripts when the quantity 
pertains to the whole chain. Therefore, using Eq. 2.5.1 we have 



7i 



/i = (iV-l)r 

_ n^+3n+4 



(„+l)3/2^r(Ar-l) 



a"^ =r{N -l){n + 1) 

(c) _ n^{n+6)-rn{n+l)+21n+2r+20 
^2 — r{N-l){n+lf ' 



From Eq. 2.4.10 we readily obtain 



(2.5.4) 



m 



ETr {AU^BUf = r'^k^ {rk + n + if . 



(2.5.5) 



By The Matching Three Moments theorem we immediately have the mean, the 
variance and the skewness for the isotropic case 



H={N -l)r 



a'^ = r{N -l){n + 1) 



7i 



n2+3n+4 



{n+lf/^yJr(N-l)' 



Note that the denominator in Eq. 2.4.2 becomes 



72 - ir 



(c) (iso) 
^4 



2E{Tr 



m 



r2 (AT - 1)2 (n + 1)2 



(2.5.6) 



In the case of Wishart matrices, m^^^ = mf^^^ = r, and m^'^'^ = mf^^^ 

and 



r (r + n + 1), m'^^^ = mfj^^ = r (r — 1) given by Eqs. 
tively. Therefore we can substitute these into Eq. 2.4.18 



2.4.10 



2.4.11 



respec- 



-E<^Tr 



m 



{AU^BuY-{AQ^BQy 



m. 



{A) 



m 



(^) 
1,1 



m 



gM{l-mE(gl)} 



/3 (m — 1) / km {n — 1) 



n [m 

2 



(m/3 + 2) 

Pk'^m^ {n — 1) 
(m/3 + 2) (m — 1) n' 



(m2 - mi^iY 



{1712 - mi^iY 



(2.5.7) 
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One can also calculate each of the terms separately and obtain the same results (see 
Appendix for the alternative). 



From Eq. 2.4.31 we have 



-E 



m 



Tr [AU^BU]' - Tr {AQ^BQ,)' 



d {2k — 1) (m2 — mi l] 



X E 



I l|2 

\uv\\y 



UV [UV] 



= {2k — 1) (m2 — nil i)^ 

X n_±Ail±^\, (2.5.8) 

We can divide Eq. |2.5.8| by Eq. |2.5.7| to evaluate the parameter p. 
2.5.2 Summary of Our Findings and Further Numerical Re- 
sults 



We summarize the results along with numerical experiments in Tables 2.6 and 



2.8 to show the equivalence of the first three moments and the departure of the three 
cases in their fourth moment. As said above, 



Qc = Irfiv Qiso = Q Haar c/^ x c?^ Qg = {Q^^Y Q 



{B) 
1 



where, {q^^^^ Q^^^ is given by Eq. 



2.3.10 



In addition, from Eq. 



2.5.4 



(2.5.9) 



we can 



define A to be the part of the kurtosis that is equal among the three cases 

A _ c 2m^m.f _ c _ 1 (rfc+n+1)^ 
^-72 ^4— - l2 2 {n+lf ■ 

Using A we can obtain the full kurtosis for the iso and quantum case, and there- 



fore (see Table 2.6 for a theoretical summary): 



P 



ll-lt 
12 -iT 



(2.5.10) 
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/3 = 1 Wishart 


Iso 


Quantum 


Classical 


Mean fi 


r{N -1) 


Variance 


r{N-l) {(f + 1) 


Skewness 71 


Vr(Af-l)(<i2 + l)3 


Tt{AQIBQ.Y 


m^m^— Eq. 2.5.7 




m^m^ - Eq. 2.5.8 




r^/c^ {rk + n + 1)^ 


Kurtosis 72*^ 




+A 


^E i:r{AQlBQ,Y 


+A 


Eq. 2.5.4 



Table 2.6: Summary of the results when the local terms are Wishart matrices. The 
fourth moment is where the three cases differ. 



The numerical convergence of the kurtoses to the theoretical values were rather 
slow. To make sure the results are consistent we did a large run with 500 million 
trials for A^ = 5, (i = 2, r = 3 and /3 = 1 and obtained four digits of accuracy 



72 - ir 



72 



0.39340 Numerical experiment 
0.39347 Theoretical value. 



Convergence is faster if one calculates p based on the departing terms alone (Eq. 
2.4.39). In this case, for full rank Wishart matrices with N = 5 and d = 2 



(3 = 1, trials: 5 Million 


1 — p 


Numerical Experiment 


0.57189 


Theoretical Value 


0.57183 



/3 = 2, trials: 10 Million 


1 — p 


Numerical Experiment 


0.63912 


Theoretical Value 


0.63938 



Below we compare our theory against exact diagonalization for various number 



of sites A^, local ranks r, and site dimensionality d (Figures 2.5.1 2.5.5). 
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Experiments based on 500000 trials 




Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean n 


8 


8.007 


8.007 


7.999 


Variance cr^ 


40 


40.041 


40.031 


39.976 


Skewness 71 


|v/lO = 1.01192 


1.009 


1.009 


1.011 


Kurtosis 72 


if = 0.590 


i = 0.660 

50 


i = 0.960 


0.575 


0.645 


0.953 


Experiments based on 500000 trials 


N = 5 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean fi 


16 


15.999 


15.999 


16.004 


Variance 


80 


79.993 


80.005 


80.066 


Skewness 71 


Ay5 = 0.716 


0.715 


0.715 


0.717 


Kurtosis 72 


H - 0.087 = 0.255 


i = 0.480 


0.085 


0.255 


0.485 


Experiments based on 300000 trials 


N = 7 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean /j, 


24 


23.000 


23.000 


24.095 


Variance cr^ 


120 


120.008 


120.015 


120.573 


Skewness 71 


0.584 


0.585 


0.585 


0.588 


Kurtosis 72 


16904 n noo 
2n6.S75 - U.082 


rk = 0.153 


i = 0.320 


-0.079 


0.156 


0.331 


Experiments based on 40000 trials 


^ = 9 


TliCMjrotical \'alue 


Nunieiical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean // 


32 


32.027 


32.027 


31.777 


Variance cr^ 


160 


160.074 


160.049 


157.480 


Skewness 71 


4v^ = 0.506 


0.505 


0.506 


0.500 


Kurtosis 72 


539142 ^ Q^g4 
3283175 i^-J-u-i 


— = 109 
800 u.±uy 


^ = 0.240 


-0.165 


0.109 


0.213 


Experiments based on 2000 trials 


= 11 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean 


40 


39.973 


39.973 


39.974 


Variance o"^ 


200 


200.822 


200.876 


197.350 


Skewness 71 


IV2 = 0.452548 


0.4618 


0.4538 


0.407 


Kurtosis 72 


-Miiii = -0.213 ^ = 0.084 


^ = 0.192 


-0.189 


0.093 


0.102 



Table 2.8: The mean, variance and skewness of classical, iso and quantum results 
match. However, the fourth moments (kurtoses) differ. Here we are showing results 
for 0? = 2, r — 4 with an accuracy of three decimal points. 



N=3, d=2, r=4, trials=500000 



0.14 
0.12 
0.1 
£r0.08 



Q 0.06 



0.04 
0.02' 
0. 







0.03 



0.02 
I 0.015 

Q 

0.01 
0.005 
0, 







0.07 
0.06 
0.05 
.&-0.04, 



0) 

Q 0.03 



0.02 
0.01 
0. 
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• Exact Diagonalization 




— I.E. Theory: p = 0.19 
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N=3, d=4, r=16, trials=50000 
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• Exact Diagonalization 






— I.E. Theory: p = 0.07 
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• Exact Diagonalization 




— I.E. Theory: p = 0.07 
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Figure 2.5.1: = 3 examples. Note that the last two plots have the same p despite 
having different ranks r. This is a consequence of the Universality Lemma since they 
have the same and d. 
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N=5, d=2, r=4, trials=500000 



1 1 1 1 

/ m * ^^Sft. * 
/ # ' ^^^^ 
/ ■ ' \^^^ 

IK f ^^\y 


— p = 

---p = 1 

• Exact Diagonalization 

— I.E. Theory: p = 0.43 
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N=5, d=3, r=5, trlals=2000 



1 1 1 1 1 1 


p = 




P=1 




• Exact Diagonalization 




— I.E. Theory: p = 0.35 
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Figure 2.5.2: = 5 
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N=7, d=2, r=4, trials=300000 
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N=9, 6=2, r=4, trials=40000 
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• Exact Diagonalization 


/i \ 


— I.E. Theory: p = 0.68 
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P=1 

• Exact Diagonalization 
— I.E. Theory: p = 0.68 
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Figure 2.5.4: = 9. Note that the two plots have the same p despite having different 
local ranks. 
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N=11,d=2, r=4, trials=2800 



p = 
p = 1 

Exact Diagonalization 
I.E. Theory: p = 0.73 
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Figure 2.5.5: = 11 



2.6 Other Examples of Local Terms 



Because of the Universality lemma, p is independent of the type of local distri- 
bution. Furthermore, as discussed above, the application of the theory for other 
types of local terms is entirely similar to the Wishart case. Therefore, we only show 
the results in this section. As a second example consider GOE's as local terms, 
i.e., Hi^i+i 
random numbers. 



'^^2'^ , where G is a full rank matrix whose elements are real Gaussian 
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N=5, d=2, r=3, trials=70000 




Figure 2.6.1: GC^'s as local terms 
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• Exact Diagonalization 
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— IE Theory: p =0.046 
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Eigenvalue 

Figure 2.6.2: Local terms have a random binomial distribution. 



Lastly take the local terms to have Haar eigenvectors but with random eigenvalues 
±1, i.e., Hi^i^i = QjAiQi, where A/ is a diagonal matrix whose elements are binary 



random variables ±1 (Figure 2.6.2) 



In this case the classical treatment of the local terms leads to a binomial distri- 
bution. As expected p = 1 in Figure [2.6. 2| has three atoms at —2,0,2 corresponding 



to the randomized sum of the eigenvalues from the two local terms. The exact di- 
agonalization, however, shows that the quantum chain has a much richer structure 
closer to iso; i.e, p = 0. This is captured quite well by IE with p = 0.046. 



2.7 Beyond Nearest Neighbors Interaction: L > 2 



If one fixes all the parameters in the problem and compares L > 2 with near- 
est neighbor interactions, then one expects the former to act more isotropic as the 



number of random parameters in Eq. 2.1.1 are more. When the number of random 
parameters introduced by the local terms, i.e., (A^ — L + 1) and are compara- 
ble, we find that we can approximate the spectrum with a high accuracy by taking 
the summands to be all isotropic [55] (See Figures 2.7.1 2.7.3). 
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-x 10 



N=6, d=4 ,L=3, r=64, trials=14 



• Exact Diagonalization 

— I.E. Theory: p = 

- - Standard Fit 
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Figure 2.7.1: IE method approximates the quantum spectrum by H^^ 

ZlLi Ql^i,- ,i+2Qi 
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• Exact Diagonalization 

— I.E. Theory: p = 

- - Standard Fit 
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Figure 2.7.2: IE method approximates the quantum spectrum by H^^ 
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• Exact Diagonalization 

— I.E. Theory: p = 

- - Standard Fit 




2000 3000 
Eigenvalue 



6000 



Figure 2.7.3: IE method approximates the quantum spectrum by H 
Ya=i QlHi,- ,1+aQi 



IE 



Most distributions built solely from the first four moments, would give smooth 
curves. Roughly speaking, the mean indicates the center of the distribution, variance 
its width, skewness its bending away from the center and kurtosis how tall and skinny 
versus how short and fat the distribution is. Therefore, it is hard to imagine that 
the kinks, cusps and local extrema of the quantum problem (as seen in some of our 



examples and in particular Figure in 2.6.2) could be captured by fitting only the first 



four moments of the QMBS Hamiltonian to a known distribution. It is remarkable 
that a one parameter (i.e., p) interpolation between the isotropic and classical suffices 
in capturing the richness of the spectra of QMBS. 

2.8 Conjectures and Open Problems 



In this paper we have offered a method that successfully captures the density 
of states of QMBS with generic local interaction with an accuracy higher than one 
expects solely from the first four moments. We would like to direct the reader's 
attention to open problems that we believe are within reach. 
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1. We conjecture that the higher moments may be analyzed for their significance. 
For example, one can show that the fraction of departing terms in the expansion 
of the higher moments (e.g. analogous to bold faced and underlined terms in 
Eqs. 2.4.4|2.4.6 but for higher moments) is asymptotically upper bounded by 
In Section2.2.3 we conjectured that their expectation values would not 
change the moments significantly. It would be of interest to know if 



ETr{...Q~^B^^QA-^Q-^B-^Q...} < 
ETr {. . . Q^'B^'Q^A^'Q-'B^'Q^ . . .} < 
ETr {...Il~^B^^UA^^U-^B^^U...} . 

For example, one wonders if 

ETr I {AQ-^BQY^ < ETr | {AQ-^BQgY^ < ETr | [AU-^BuY^ 

for k > 2] we have proved that the inequality becomes an equality for k = 1 
(Departure Theorem) and holds for A; = 2 (Slider Theorem). 



2. Though we focus on a decomposition for s, we believe that the main theorems 
may generalize to higher dimensional graphs. Further rigorous and numerical 
work in higher dimensions would be of interest. 



3. At the end of Section 2.4.4 we proposed that more general local terms might 



be treated by explicitly including the extra terms (Type I terms). 

4. Application of this method to slightly disordered systems would be of interest 
in CMP. In this case, the assumption of fully random local terms needs to be 
relaxed. 

5. In our numerical work, we see that the method gives accurate answers in the 
presence of an external field. It would be nice to formally extend the results 



86 



and calculate thermodynamical quantities. 

6. We derived our results for general (3 but numerically tested (3 = 1,2. We 
acknowledge that general /3 remains an abstraction. 

7. Readers may wonder whether it is better to consider "iso" or the infinite limit 
which is "free". We have not fully investigated these choices, and it is indeed 
possible that one or the other is better suited for various purposes. 

8. A grander goal would be to apply the ideas of this paper to very general sums 
of matrices. 



2.9 Appendix 

To help the reader with the random quantities that appear in this paper, we provide 
explanations of the exact nature of the random variables that are being averaged. 
A common assumption is that we either assume a uniformly randomly chosen eigen- 
value from a random matrix or we assume a collection of eigenvalues that may be 
randomly ordered. Random ordering can be imposed or a direct result of the eigen- 
vector matrix having the right property. Calculating each of the terms separately 
and then subtracting gives the same results. 



-E 



m 



Tr {AQ^BQY 



-E 



m 



^ ^ Ctj-^ djj i.QiljlQilj2^i2jl^i2J2/ 



f2.9.11 



where and bj are elements of matrices A and-B respectively. The right hand side of 
Eq. 2.9.1 can have terms with two collisions (i.e., ii = 12 and ji = j2), one collision 



(i.e. ii 7^ ^2 exclusive-or ji 7^ ^2), or no colhsions (i.e., ii 7^ 22 and ji 7^ 72)- Our goal 
now is to group terms based on the number of collisions. The pre-factors for two. 



one and no collisions along with the counts are summarized in Table 2A Using the 
latter we can sum the three types of contributions, to get the expectation 
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TT{AQ^BQy 



(/9+2) 
(m/3+2) 



E(a2)E(62) + 



fc§ [E (62) E (aia^) + E {a') E (6162)] - gg^E (a^a^) E (6162) • 



(2.9.2) 



If we take the local terms to be from the same distribution we can further express 
the foregoing equation 



-E 



m 



TT{AQ^BQy 



[(/3 + 2) + /3 (m - 1) E (aiOs) {2m2 - E (aiOa)}] 



(m/3 + 2) 

(2.9.3) 

The quantity of interest is the difference of the classical and the iso (see Eq. 



2.4.8) 



^ETr {AU^BUY - ^ETr {AQ^BQf 



{E (a2) E (62) - E (62) E (aiag) - E {a^) E (6162) + E (0102) E (6162)} 



(m/3+2) 



/3(m-l) 
(m/3+2) 



1712 "^2 "~ ""^2 ''^i 1 ~ ^2 "^i 1 + "^1 1 "^1 1 



(m/3+2) 



(^) 



(A) 
1,1 



{B) 



m 



1,1 



(2.9.4) 



If we assume that the local terms have the same distribution: m2 = m 



(A) 



m. 



(B) 

2 ' 



mil = 1 = ^i 1 as m bq 
1 



2.9.3 



ETr {Ali^BIi) ETr [AQi^ BQ) 

m m 



the foregoing equation simplifies to 

/3(m- 1) 



(m/3 + 2) 

In the example of Wishart matrices as local terms we have 



(m2-mi_i) . 



m2 = E (a2) = rk {rk + n + 1) 
mi^i = E (aia2) = k{k-l)r^ + ^{ {tn^-^ - l) (n + r + 1) + tn''~^ {n - 1) (r - 1)} 

= A; (A; - 1) r2 + ^ {tn^-^{nr + 2)-n-r-l}. 

(2.9.5) 
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Chapter 3 



Calculating the Density of States 
in Disordered Systems Using Free 
Probability 

In the previous chapter we saw that when a decomposition into two commuting sub- 
sets is possible, one can with a high accuracy obtian the density of states of quantum 
spin chains. Spin chains that we discussed do not possess transport properties such 
as hopping of electrons. In this chapter we like to extent the ideas of decomposing 
the Hamiltonian to 'easier' pieces and treat different systems. In particular we will 
focus on one particle hopping random Schrodinger operator. What follow in the rest 
of this chapter also apprears in [56] . 
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3.1 Introduction 



Disordered materials have long been of interest for their unique physics such as 
localization [SZIEH], anomalous diffusion [321110] and ergodicity breaking [HI]. Their 
properties have been exploited for applications as diverse as quantum dots |62l [63] , 
magnetic nanostructures [6l], disordered metals [651 ES], and bulk hetero junction 
photovoltaics [671 EH EH] • Despite this, theoretical studies have been complicated by 
the need to calculate the electronic structure of the respective systems in the presence 
of random external potentials. Conventional electronic structure theories can only be 
used in conjunction with explicit sampling of thermodynamically accessible regions 
of phase space, which make such calculations enormously more expensive than usual 
single-point calculations [70] . 

Alternatively, we aim to characterize the ensemble of electronic Hamiltonians that 
arise from statistical sampling directly using random matrix theory; this would in 
principle allow us to sidestep the cost of explicit statistical sampling. This naturally 
raises the question of whether accurate approximations can be made for various char- 
acteristics of random Hamiltonians such as their densities of state (DOSs). We use 
techniques from free probability theory, which allow the computation of eigenvalues 
of sums of certain matrices [7T]. While this has been proposed as a tool applicable 
to general random matrices [72] and has been used for similar purposes in quantum 
chromo dynamics [73] , we are not aware of any quantification of the accuracy of this 
approximation in practice. We provide herein a general framework for quantitatively 
estimating the error in such situations. 
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3.2 Quantifying the error in approximating a PDF 
using free probability 



We propose to quantify the deviation between two PDFs using moment expansions. 
Such expansions are widely used to describe corrections to the central limit theorem 
and deviations from normality, and are often applied in the form of Gram-Charlier 
and Edgeworth series [711 ES]- Similarly, deviations from non-Gaussian reference 
PDFs can be quantified using generalized moment expansions. For two PDFs w (^) 
and w (^) with finite cumulants Hi, K2, ■ ■ ■ and Ki, K2, . . . , and moments /ii, /i2, . . . and 
/ii,/i2, . . . respectively, we can define a formal differential operator which transforms 
w into w and is given by [76], [71] 



w (0 = exp 



E- 



71=1 



(3.2.1) 



This operator is parameterized completely by the cumulants of both distributions. 

The first k for which the cumulants Kk and kk differ then allows us to define 
a degree to which the approximation w ^ w is valid. Expanding the exponential 
and using the well-known relationship between cumulants and moments allows us to 
state that if the first k — 1 cumulants agree, but the kth cumulants differ, that this 
is equivalent to specifying that 



«; (0 = * (0 + 



k\ 



(3.2.2) 



At this point we make no claim on the convergence of the series defined by the 



expansion of (3.2.1 ), but use it as a justification for calculating the error term defined 



in (3.2.2). We will examine this claim later. 
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3.3 The free convolution 



We now take the PDFs to be DOSs of random matrices. For a random matrix Z, the 
DOS is defined in terms of the eigenvalues of the M samples Zi, . . . , Z^^, . . . , Zm 

according to 

^<^'K)=to„^E^E^(«-e'). (3.3.1) 

m=l n=l 

The central idea to our approximation scheme is to split the Hamiltonian H — A + B 
into two matrices A and B whose DOSs, p^"^^ and p^^^ respectively, can be determined 
easily. The eigenvalues of the sum is in general not the sum of the eigenvalues. 
Instead, we propose to approximate the exact DOS with the free convolution ASB, 
i.e. p^^^ ~ a particular kind of "sum" which can be calculated without exact 

diagonalization of H. The moment expansion presented above quantifies the error of 
this approximation in terms of the onset of discrepancies between the kth moment 
of the exact DOS, iJ.\f\ and that for the free approximant P'^^^^^ ■ By definition, the 
exact moments are 

l^f'-f^r'^'-iiA + B)'), (3.3.2) 

where (Z) = E Tr (Z) /N denotes the normalized expected trace (NET) of the NxN 
matrix Z. The kth moment can be expanded using the (noncommutative) binomial 
expansion of {A + B)'^; each resulting term will have the form of a joint moment 
{A^'^B'^'^ ■ • ■ A^^'B'"'') with each exponent ng^nis being a positive integer such that 
Yll=i (^s + ^s) = k. The free convolution pk is defined similarly, except that A and 
B are assumed to be freely independent, and therefore that each term must obey, by 
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definition [77] , relations of tlie form 



= (n;=i - (5™= - (5™=))) (3.3.3a) 

= (n;^iA"^5'"^) + lower order terms, (3.3.3b) 

where the degree k is the sum of exponents Ug, ms and the second equality is formed 
by expanding the first line using linearity of the NET. Note for A; < 3 that this 
is identical to the statement of (classical) independence [77j. Testing for /i^"^^^^ 7^ 
A^fc^^^^ then reduces to testing whether each centered joint moment of the form in 



(3.3.3a) is statistically nonzero. The cyclic permutation invariance of the NET means 
that the enumeration of all the centered joint moments of degree k is equivalent to 
the combinatorial problem of generating all binary necklaces of length fc, for which 
efficient algorithms exist [80] . 



The procedure we have described above allows us to ascribe a degree k to the 
approximation p'^^'^ ^ p(-4fflB) given the splitting H = A + B. For each positive 
integer n, we generate all unique centered joint moments of degree n, and test if they 
are statistically nonzero. The lowest such n for which there exists at least one such 
term is the degree of approximation k. We expect that > 4 in most situations, 
as the first three moments of the exact and free PDFs match under very general 
conditions |HI]. However, we have found examples, as described in the next section, 
where it is possible to do considerably better than degree 4. 
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3.4 Decomposition of the Anderson Hamiltonian 



As an illustration of the general method, we focus on Hamiltonians of the form 



hi J 
J h2 



■•• J 

J Hn 



(3.4.1) 



where J is constant and the diagonal elements hi arc identically and independently 
distributed (iid) random variables with probability density function (PDF) Ph{i)- 
This is a real, symmetric tridiagonal matrix with circulant (periodic) boundary con- 
ditions on a one-dimensional chain. Unless otherwise stated, we assume herein that 
hi are normally distributed with mean and variance cr^. We note that cr/J gives 
us a dimensionless order parameter to quantify the strength of disorder. 



So far, we have made no restrictions on the decomposition scheme H — A + B 
other than p^^^ and p^^^ being easily computable. A natural question to pose is 
whether certain choices of decompositions are intrinsically superior to others. For 
the Anderson Hamiltonian, we consider two reasonable partitioning schemes: 



H = Ai + Bi 



hi 



\ 



h. 



h. 



( 



+ 
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J 
J J 
J 



(3.4.2a) 



H = A2 + B2 



hi J 
J 



h3 J 
J 



V 



+ 



h2 J 
J 



V 



(3.4.2b) 



We refer to these as Scheme I and II respectively. For both schemes, each fragment 
matrix on the right hand side has a DOS that is easy to determine. In Scheme I, 
we have pAi = Ph- Bi is simply J multiplied by the adjacency matrix of a one- 
dimensional chain, and therefore has eigenvalues A„ = 2J cos{2n7c/N) [82j. Then 
the DOS of Bi is p^i (0 = J2n=i ^ ~ -^n) which converges as — )• oo to the arcsine 
distribution with PDF pAs (0 = V {t^V^J^~^-^) on the interval [-2 | J| , 2 | J|]. In 

hi J 
J 

Since X has eigenvalues e± {^) = hi {^) /2 ± \/hu{ (,^) /4 + J^, their distribution can 
be calculated to be 



Scheme II, we have that pA^ = PB2 = Px where px is the DOS of X 



Px (0 



(3.4.3) 



3.5 Numerical free convolution 

We now calculate the free convolution A ffl i? numerically by sampling the distribu- 
tions of A and B. We define as simply the average DOS of the free approxi- 
mant Z = A + Q^^BQ, where Q is a x random matrix of Haar measure. For 
real symmetric Hamiltonians it is sufficient to consider orthogonal matrices Q, which 
can be generated from the QR decomposition of a Gaussian orthogonal matrix [78] . 
(This can be generalized readily to unitary and symplectic matrices for complex and 
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21 



-2.5 2.5 



-2.5 2.5 



Figure 3.5.1: Calculation of the DOS, p{^), of the Hamiltonian H of (3.4.1) with 
M = 5000 samples of 2000 x 2000 matrices for (a) low, (b) moderate and (c) high 
noise {a/J=0.1, 1 and 10 respectively with a = 1). For each figure we show the 
results of free convolution defined in Scheme I (p(^i^^i); black solid line). Scheme II 
^p(^2ffl-B2). gYQQY]^ dashed line) and exact diagonalization (p^^^; red dotted line). 



quaternionic Hamiltonians respectively.) In the iV — )■ oo limit, this converges to the 
free convolution ASB jZUCS]. 



The exact DOS p(^+^) and free approximant p(^^^) are plotted in Figure 3.5.1 'a)- 



(c) for both schemes for low, moderate and high noise regimes (cr/J =0.1, 1, 10 
respectively). 

We observe that for Scheme I we have excellent agreement between p^^^ and 
p(Aiffl_Bi) g^Qj-Qgg all values of cr/J, which is evident from visual inspection; in contrast. 
Scheme II shows variable quality of fit. 

We can understand the starkly different behaviors of the two partitioning schemes 
using the procedure outlined above to analyze the accuracy of the approximations 
^(AifflBi) ^j^^ p{H) ^ p(A2ffl-B2)^ YoT Scheme I, we observe that the ap- 



P 



proximation (3.2.2) is of degree k = 8; the discrepancy lies solely in the term 



'(Aii?i)^) [H3]- Free probability expects this term to vanish, since both Ai and 



Bi are centered (i.e. {Ai) = (Bi) = 0) and hence must satisfy (3.3.3b) with 
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--CO- 




+ n lT3 V 




Figure 3.5.2: Diagrammatic expansion of the term (AiBiAiBiAiBiAiBi) in terms 



of allowed paths dictated by the matrix elements of Ai and Bi of Scheme I in (3.4.2a) 



rii = mi = ■ ■ ■ = 77-4 = 7714 = 1. In contrast, we can calculate its true value 
from the definitions of Ai and Bi. By definition of the NET (■), only closed paths 
contribute to the term. Hence, only two types of terms can contribute to (^{AiBi)^y, 



these are expressed diagrammatically in Figure 3.5.2 The matrix Ai weights each 
path by a factor of h, while Bi weights each path by J, and in addition forces the 
path to hop to an adjacent site. 

Consequently, we can write explicitly 

i 

+ ^ ^ E (hiJhi+iJhiJhi^iJ) 

i 

+ {hiJK.iJhiJhi^iJ) 

i 

+ ^ ^ E (hiJhi+iJhiJhi+iJ) 

i 

=2J%(/i,)'E(/i2) + 2J^E(/i2)' = + 2JV, (3.5.1) 
where the second equality follows from the independence of the /ij's. As this is 
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the only source of discrepancy at the eighth moment, this explains why the agree- 
ment between the free and exact PDFs is so good, as the leading order correction 
is in the eighth derivative of p(^iffl^i) with coefficient 2aV^/8! = (cxJ)^ /20160. In 
contrast, we observe for Scheme II that the leading order correction is at = 4, 
where the discrepancy lies in {A^B"^)- Free probability expects this to be equal to 
{AlBl) = {AD {El) = {X^^ = (J2 + a^/2f, whereas the exact value of this term 
is (J^ + cr^). Therefore the discrepancy is in the fourth derivative of p^^^^^ with 
coefficient (-aV4) /4! = -(7^96. 

3.6 Analytic free convolution 

Free probability allows us also to calculate the limiting distributions of p(^^^) in the 
macroscopic limit of infinite matrix sizes N ^ oo and infinite samples M — )• oo. In 
this limit, the DOS p^^^^^ is given as a particular type of integral convolution of 
p^^^ and p^^\ We now calculate the free convolution analytically in the macroscopic 
limit for the two partitioning schemes discussed above, thus sidestepping the cost of 
sampling and matrix diagonalization altogether. 

For our first example, we take A and B as in Scheme I, but with each iid hi 
following a Wigner semicircle distribution with PDF pw (^) = 1^4 — ^'^/4:Tt on the 
interval [—2,2]. (Using semicircular noise instead of Gaussian noise simplifies the 
analytic calculation considerably.) Then, p*^^^ = pw and p^^^ = pas- The key tool 
to performing the free convolution analytically is the /^-transform r {w) = (w) — 
[84j, where is defined implicitly via the Cauchy transform 




(3.6.1) 
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Figure 3.6.1: DOS, p(^), of the Hamiltonian (3.4.1) with M = 5000 samples of 
2000x2000 matrices with (a) low, (b) moderate and (c) high semicircular on-site noise 
(cr/ J=0.1, 1 and 10 respectively with a = 1), as calculated with exact diagonalization 
(red dotted line), free convolution (black solid line), and perturbation theory with 
Ai as reference (blue dashed line) and Bi as reference (gray dash-dotted line). The 
partitioning scheme is Scheme I of (3.4.2a). 



For freely independent A and B, the /^-transforms linearize the free convolution, 
i.e. r^"^^^) (^w) = r^^'^ (w) + r^^^ (w), and that the PDF can be recovered from the 
Plemelj-Sokhotsky inversion formula by 

p(^=^)(e) = ^Im((^(^=^))-^(e)) (3.6.2a) 
g(^^B) (w) = r(^=^) (w) + w-\ (3.6.2b) 

Applying this to Scheme I, we have r*^"^) [w) = w and r*^^^ (w) = (— 1 + y/AJ^ + w^^ / w, 
so that g^'^^^') (w) = w + (v^4J2 + w^) /w. The need to calculate the functional 
inverse ^ in this procedure unfortunately precludes our ability to write 

piAWB) (^^^ in a compact closed form; nevertheless, the inversion can be calculated 
numerically. We present calculations of the DOS as a function of noise strength a/ J 
in Figure |3.6.1 
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3.7 Comparison with other approximations 



For comparative purposes, we also performed calculations using standard second- 
order matrix perturbation theory [85] for both partitioning schemes. The results 



are also shown in Figure 3.6.1 Unsurprisingly, perturbation theory produces results 
that vary strongly with cr/J, and that the different series, based on whether A is 
considered a perturbation of B or vice versa, have different regimes of applicability. 
Furthermore it is clear even from visual inspection that the second moment of the 
DOS calculated using second-order perturbation theory is in general incorrect. In 
contrast, the free convolution produces results with a more uniform level of accuracy 
across the entire range of cr/J, and that we have at least the first three moments 
being correct [8T| . 

It is also natural to ask what mean-field theory, another standard tool, would 
predict. Interestingly, the limiting behavior of Scheme I as — oo is equiva- 
lent to a form of mean-field theory known as the coherent potential approximation 
(CPA) [HHl EH EH] in condensed matter physics, and is equivalent to the Blue's func- 
tion formalism in quantum chromodynamics for calculating one-particle irreducible 
self-energies [73]. The breakdown in the CPA in the term is known [891157]; 

however, to our knowledge, the magnitude of the deviation was not explained. In 
contrast, our error analysis framework affords us such a quantitative explanation. 

Finally, we discuss the predictions of isotropic entanglement theory, which pro- 
poses a linear interpolation between the classical convolution 

POO 



and the free convolution p(^^^) (^) in the fourth cumulant [SD ED]- The classi- 
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cal convolution can be calculated directly from the random matrices A and B; by 
diagonalizing the matrices as A = Q^^AaQa and B = Q]^^AbQb, the classical con- 
volution p(^*^) (0 can be computed from the eigenvalues of random matrices of the 
form Zci = Aa + n~^A5n where 11 is a x random permutation matrix. It is 
instructive to compare this with the free convolution, which can be sampled from 
matrices of the form Z' = Aa + Q^^AbQ, which can be shown by orthogonal invari- 
ance of the Haar measure random matrices Q to be equivalent to sampling matrices 
of the form Z = A + Q^^BQ described previously. 

As discussed previously, the lowest three moments of Z and H are identical; this 
turns out to be true also for Z^i [81]. Therefore IE proposes to interpolate via the 
fourth cumulant, with interpolation parameter p defined as 



We observe that IE appears to always favor the free convolution limit {p = 0) 
as opposed to the classical limit {p = 1); this is not surprising as we know from 
our previous analysis that 4 , and that the agreement with the exact 

diagonalization result is excellent regardless of a/ J. In Scheme II, we observe the 
unexpected result that p can sometimes be negative and that the agreement varies 
with a/ J. From the moment expansion we understand why; we have that the first 
three moments match while k^^^'^^''^^ — k^^^^^^'^ = —a^ /A. The discrepancy lies in the 
term {A^B^) which is expected to have the value {A2) (-B|) = {-P + in free 

probability but instead has the exact value (J^ + a^). Furthermore, we have that 
^{A2*B2) _^ /t^"^^ ^-^2 ^ where the only discrepancy lies is in the so-called departing term 
l^A^B^A^B^) [m SO]. This term contributes to 4^^^^ but has value (A^) = 



p = 




(3.7.1) 
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(J2 + a^/2f in 4 



(^2*62) 



, since for the classical convolution we have that 



(3.7.2) 



This therefore explains why we observe a negative p, as this calculation shows that 



which is manifestly negative. 

In conclusion, we have demonstrated that the free probability of random matrices 
can provide unexpectedly accurate approximations for the DOS of disordered Hamil- 
tonians, both for finite dimensional systems and in the macroscopic limit N ^ 00. 
Our results illustrate variable accuracies of the approximations predicated on our 
particular choices of partitioning schemes, which can be quantitatively estimated 
using moment expansions. These results represent an optimistic beginning to ad- 
dressing the electronic structure of disordered condensed matter systems using the 
tools of random matrix theory. We are currently investigating the predictions of 
free probability on the localization properties of eigenvectors, as well as studying the 
general applicability of free probability to lattices in higher dimensions, as well as 
systems with off-diagonal and correlated disorder. 




(3.7.3) 
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Part II 
Eigenvectors 
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Chapter 4 

Generic Quantum Spin Chains 



In this chapter I discuss frustration free condition, evolution in time and imaginary 
time within MPS representation along with numerical methods. I then discuss the 

degeneracy and frustration free condition for generic quantum spin chains with local 
interactions. I leave the discussion of the entanglement of the ground states for the 
next chapter. 

4.1 Frustration Free Condition 

An L— local Hamiltonian where each interaction acts nontrivially on L particles can 
be written as 



1=1 

where i indexes M groups of L spins. By spectral decomposition any such Hamilto- 
nian on N spins each with (i— states can be written as 



M 
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H = rrigXg \g){g\ + ^ Aa \a){a\ , withamax + rrig = d^^ 

a=l 

where g refers to the ground state quantities, rrig the possible multiphcity due to 
degeneracy and a's label the excited states. Similarly for each of the local terms 

r 

H,^mlXl\g,){gi\ + J2K\'^i){^\ withr + m^ = d^. (4.1.1) 

p=i 

Energy is meaningful in a relative sense, hence one can shift the local Hamiltonian 
such that the ground state has energy zero. This defines an effective Hamiltonian H'- 

p=l 

The effective local Hamiltonian has zero energy ground states and eigenstates 
corresponding to the excited states with positive eigenvalues. By definition, A^ — A* > 
for all p, where we have not made as of yet any assumptions on the possible 
numerical values they can take. The span and Kernel of H'^ ^ is therefore the same 
as 

^".=i:i<>«i, (4.1.2) 

p=l 

where we replaced all Ap — > by ones to get projectors of rank r as our local 
terms. As discussed at the end of this section, there are advantages in doing so. 
Locally there are (f — r zero energy states; however, the projectors are not mutually 
exclusive. This makes the problem of counting the kernel of H non-trivial, namely 
we cannot assert that the number of ground states is {(P — r)^ . In particular the 
ground state of the global Hamiltonian may not be a local ground state of some local 
term(s): frustrated. 
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Lemma 4. Suppose E^^. and E^^ are the smallest eigenvalues of Hamiltonians 
and respectively and E^i^ that of + H^. Further, let E^^^^ and E^^^ be the 
largest eigenvalues of Hamiltonians and H"^ respectively andEl^^,j. that ofH^ + H"^ 
then 



E' +E^ < E''^ 

gr ' gr — gr 



/?! _|_ rr^ \ pi 



nl,2 

"max ' "-"max — "-"max" 



In the first (second) case equality holds if and have the same eigenstate for 
their smallest (largest) eigenvalues. 

Proof. The proof follows from the fact that the average is greater than the least 
eigenvalue and smaller than the largest eigenvalue. Suppose we are in the states in 
which H^+H'^ takes the eigenvalue E^^^, then {H^ + H"^) = E];^ , but {H^ + H"^) = 
(H^) + {H"^) > Eg^ + Egj.. The equality holds if and if^ have the same eigenstate 
for their smallest eigenvalues. Similarly if we are in the state in which + takes 
the eigenvalue E^il^, then {H^ + H^) = E^il^ , but {H^ + H^) = {H') + (H^) < 
^max + ^max- The equality holds if and have the same eigenstate for their 
largest eigenvalues. □ 

One can prove this using Rayleigh quotient or min-max theorems as well. 



Therefore, the summation in the first Eq. 4.1.2 can result in "lifting" the global 



ground state relative to that of the local ones. We seek conditions under which global 



ground state has zero energy, i.e., remains "unlifted" (see Section 4.2 for motivation). 

It is necessary and sufficient for a state to be the ground state if it is orthogonal 
to all the local terms H'^ because all the terms in the Hamiltonian will be zero. If the 
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local terms commute by the foregoing lemma the global and local ground states will 
be the same, i.e., E^:^ = Eg^ + E^^. = + = 0. This is too strong a requirement. 
It could be that the lowest eigenvectors of each summand are aligned, in which case 
the system is Frustration Free (FF) or unfrustrated. A classical analogue would be 
the ferromagnet. 

Definition 1. The ground state is unfrustrated if it is also a common ground state 
of all of the local terms H^. 



Let us go back to Eq. |4.1.3| and investigate chains of rf-dimensional quantum 
spins (qudits) on a line with nearest-neighbor interactions. The Hamiltonian of the 

system, 

H=Y,Hk,k+i (4.1.3) 
fc=i 

is 2-local (each Hk^k+i acts non-trivially only on two neighboring qudits). Our goal 
is to find the necessary and sufficient conditions for the quantum system to be un- 
frustrated 

The local terms can be written as 

Hk,k.. = Ef^P^l, + Y: E^'^^pI-U (4-1-4) 



where Eq''^ is the ground state energy of Hk^k+i and each Pjfl^i is a projector onto 
the subspace spanned by the eigenstates of Hk,k+i with energy E^t"^ . The question 
of existence of a common ground state of all the local terms is equivalent to asking 
the same question for a Hamiltonian whose interaction terms are 



— ® Pk,k+1 ® Ifc+2,...,Ar, 
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(4.1.5) 



with Pk,k+i = J2p=i Pk'k+i projecting onto the excited states of each original interac- 
tion term H^ k+i- When this modified system is unfrustrated, its ground state energy 
is zero (all the terms are positive semi-definite). The unfrustrated ground state be- 
longs to the intersection of the ground state subspaces of each original Hk^k+i and is 
annihilated by all the projector terms. 

As far as the question of (un)frustration and count of the ground states are 



concerned, Eq. 4.1.5 yields the same result as Eq. 4.1.4 



4.2 Why Care About Frustration Free Systems? 

There are many models such as the Heisenberg ferromagnetic chain, AKLT, Parent 
Hamiltonians that are frustration free (FF) |116[ 11171 SU [33] • Besides such models 
and the mathematical convenience of working with projectors as local terms, what is 
the significance of FF systems? In particular, do FF systems describe systems that 
can be realized in nature? Some answers can be given: 

1. It has been proved in |1U6] all gapped Hamiltonians can be approximated by 
frustration free Hamiltonians if one allows for the range of interaction to be 
0(\ogAf). It is believed that any type of gapped ground state is adequately 
described by a frustration free model [H]. 

2. The ground state is stable against variation of the Hamiltonian against pertur- 
bations [95j 

Hi9) = ^ 9kHk,k+i , 9k >0 

k 

as the kernels of the local terms remain invariant. 

3. In quantum complexity theory the classical SAT problem is generalized to the 
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so called qSAT |103] . The qSAT problem is: Given a collection of /c— local pro- 
jectors on n qubits, is there a state ip) that is annihilated by all the projectors? 
Namely, is the system frustration free? 



4. Ground states of frustration free Hamiltonians, namely MPS can be prepared 
by dissipation [94]. 

4.3 Why Does Imaginary Time Evolution Work? 

By imaginary time evolution we mean it — )■ r. The Hamiltonian evolution becomes 

Intuitively one might like to see imaginary time evolution as dissipation of energy 
such that for sufficiently long time the system relaxes to its lowest energy state. 
Mathematically, one can do a spectral decomposition H = where 
are energies associated with states \a). The imaginary time evolution 

which implies an exponential suppression of the overlap of {ip) with states that have 
energies higher than that of the ground state. My numerical implementation of the 
imaginary time evolution is the same as described in 
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4.4 Numerical Study of Quantum Spin Chains Us- 
ing MPS 



Suppose we want to evolve the MPS representation of the quantum spin system, (to 



be exphcit we include the A's see Eq. 1.3.9) starting at time t 



■N-l 

i^t)='P{(^T,{t,)A'f^\z,)\ OBC. 
p=i 



The quantum mechanical time evolution in At is given by 

i^t+At) = |0exp(-«Atifp,p+i)| \iJt). 

As in Chapter 2, we once again decompose the Hamiltonian into two pieces H = 

Hi + H2, where Hi = Z]p=i,3,... -^p,p+i -^2 = Z]p=2,4,- -^p,p+i made up of 
terms that all commute with one another. One can show [391 Exercise 4.47] 



Q-itHi _ g-it-f/i,2g-jt-ff3,4 . . . ^-itHff-2,N-l ^ 
Q-itH2 _ g-jt-f^2,3g-jt-f^4,5 . . . g-«t-ffjV-l,JV 

In order to evolve the system in time we need to make use of Trotter's formula 
[391 Thm 4.3] 

Theorem. (Trotter formula) Let A and B he Hermitian operators. Then for any 
real t, 
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(4.4.2) 



Among other things, one can use the proof of this theorem to show [39| see Eq. 
4.103], 



(4.4.3) 



where we used A = Hi and B = H2. 
Using Eqs. |4.4.1| and |4.4.3 



i^t+At) = < (g) exp i-iAtHp^p+i) (g) exp {-lAtHp^p+i) \ ^Pt) + O (At^) 
Uodd peven J 

which imphes that Trotterization and evolution in time can be implemented by evolv- 
ing the even terms first and then the odd terms for small time intervals. 

An advantage of MPS is that one can locally update the state ^97, Lemma 2], 
i.e., apply the operator on the local terms n and n + 1, we first define 

Let the imaginary time evolution operator be 



'ni'n + l V 'n''7i + l / 



which after the imaginary time evolution update reads 
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To get the updated Matrix Products we use SVD 



Tt+At {^n) A^Atr^+At (^n+i) = SVD (qIXa^) (4.4.5) 

The state after imaginary time evolution will in general not be normalized. As 
in [96], I normalize the state after every update step. The work in [96] was confined 
to translationally invariant chains with d = 2 where F's and A's could be taken to 
be the same at every site and bond respectively. I generalized the numerical code to 
implement the updates for any d and regardless of whether there was translational 



invariance. A major hinderance was implementation of G update (Eq. 4.4.4) in an 
efficient manner. I overcame this by first reshaping the matrices and applying the 
update rules followed by undoing of the reshaping. I then took the SVD to update 
the chain and normalized the state. This was done for every pair of nearest neighbors 
qudits in Hi followed by updates on every pair of nearest neighbor qudits in H2 (see 
the code in Algorithm [T]). See the section 4.7 for further discussion and numerical 
results. 
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Algorithm 1 This code updates Eqs 



4.4.4 



an 



% Kami 5 Hovasgagli 

% ApiiL 2012 

% AIL riglita reserved 

% writing tne portion of tafi code for calculating Tneta and its update 
funct ioa [at ate 

Loa annrin ] -recQmpiiteGL2_'re atTr ivial ( Jian, state, EtrengtJi,bitl,bit2, boiid_a ) 
global imaginaryt iiiie 
gleJbal uaeswd 

global florjrtbyjmiLtipLiGatioa 
global Cutoff Lajribda 

% interact a tUe bonda a connecting bitl to bit2 

% note I bond g is not interacted 

i 

i 

i b a c 

% bitl — bits — 

1-atate .Knaij 

dit - atate.ditj 

I b.ardcDde the Line Here. . . 

bQnd_b - bQQd_a - If 

bQnd_c - bQad_a + If 

% bond indices in -[bitl, bit 2 & bit3> 
al-f ind ( at ate . opt ions . bonds (bitl,i) — bond_a)j 
bl- f ind ( St ate . opt ions . bonds (bitl,i) — bond_b)|[ 
a2- f ind ( St ate . opt ions . bonds (bit2,i)~bond_a3| 
c2-find(Bt ate . opt iona . bonda ( bit 2 , ■ ) --bQnd_c ) f 

i lanlbdaa of interacting bonda 
La-diag( atate.L(bQnd_a, ■ ] ) ; 
Lb-atate.L(bond_b, i ) j 
liC-Btate .Ei(boad_c, ■ ) j 

% tnis is for tne interaction of nearest aitea. Tlieae are jas.t tne 
ganunaa i 

Hba - aeroa (Iflfdit ) ( 
Hac - aeroa (l,l,dit ) p 
r - (l.dat), 

Hba( I , I ,r3 - atate.H( 1 , 1 ,r,bitl J p * mis is Ganmia_ab 
Haic(ififr) — atate.H( 1 , 1 ,r ,bit2 J ( % Tliia ia G-aiimia_ac 
% we need inverae vectora for La and Li 
nonaeroLa - f ind( atate.L(bond_a, i ) xjutofElanbda) j 
liail-zeros (X, 1 ] f 

Lai 1 ( nonseroLa ) - 1 . / at ate . L ( boad_a , nonseniLa) f 

Lai - diag(Lail ) p 

% b a c 

% bitl bit2 

bb-expni(— atrengtti. *liaiii( 1 , 1 ,bitl ) ) p * It ia called V in ttue ttieaia 

■met aVec - =eroa (dit*dit*l*l,l)p 
tnet a_be f ore - zeros (1,1, dit , dit ) p 
for a - lidit 

for t- lidit 

%at evey entry s,t "nieta naa a matriic tHat Ib X k X. 



% Ton rion't iified to store tnis: FOE CHECKING 




imfita - E; 




%met a - Bba ( : , : , a ] ^Lai ^Bac ( : , : , t ) ; 




fttnet a_bie f ore ( : , : , s ,t ) - ■cneta; 




ttiet a_bef ore (:,:,B,tJ — Hba(:,:,B)*Lai*Hacfi,i,t)) 




% s r SE_rHErA- s i se ( met a j 




%Here Lai is tne inverse of Lambda in between tHe sites. 


Eba 


and 




%Hac are just tHe L Gajimia L eacn. 




tHetamter - reanape ( tHet a_be£Qre ( : , : , s , t ) , [ ] , 1 ) ; 




THet aVec ( ( p- 1 ) 1 : p*l*l ) -tHet al nter ; 




P^p4-l; 




end 




end 




met a2 - re snape ( ■niet avec , i*i , [ ] ) j 




TTtifitaSt - TtietaS' ; 




THet a_iipdate - bb*Tlieta2t ; 4 returns a dit'2, by 




% Here every coLuJim contains a Loop over apina (Eiiced pair of 1) 




tUGt a_iipdat62t - ■Elieta_updat6' ; * by dit*2 




riTD - sise (tlieta_update2t ) ; 




tHetaB - aeros(riirD(2) ,1,1) ; 




for i - l:rin)(2) 




tli6taE(i, : , : ) - resnape(tlieta_iipdate2t( i , i) , K, K) j 




end 




tHet a_iipdat6_back - aenoB (l,l,dit, dit) j 




% Heed to make it a sq-uaire with diiaeiiaioiia }:*dit by K*dit 




'nieta_flat-=eros(dit*l,dit*l) ; * Will be tHe updated THeta 




for a - l:dit 




for t - l:dit 




tlieta_update_back( : , : , s ,t ) - tlietaE( ' i ']} % Can be omitted 


™fita_f Lat( (s-l)*l+l:s*l, (t- 




1 ) *I-4-l : t *1 ) -tHet a_update_back ( : , : , a , t ) j 




end 




end 




%% use tne SVD for decomposition of Tneta 




[VV,EE,mf] - Bvd('niifita_fLat, 'econ' ) ; % THeta - W * EE * HW' 




LL - diag(EE ) ; 




%Note on Size: W and EE and HH ME ALL (dj.t*K) x (dj.t*X)i aaaie 


as 


met a 




%LL rS (dit*!) I 1 




[LLb , How] -sort (LL, 'descend' ) ; 




i.e LLs - LL(nov) ; 




LOHsnorm - hujii(LLb (l+l :end) . *2 ) ; 




flt ate . L ( bond a , : ) — ; 




i renormalize tUe lanlbdas 




Etate.L(T-ond a, : )-LL(liow( 1 :1) ) ; 




atat6.L(bQnd_a, : )-atate.L(bQnd_a, i ) /1101M15 state. L( boiid_a, ■ ] ) ; 




EEcutoff - diag( state. L(bQad_a, :)) ; 




for r - l:dit 




state.H( : , : ,r,bitl) - W(l* (r-1 :r*l,lioif ( 1 :1) ) * EEcutoff; 




state.H( : , : ,r,bit2) - EEcutoff * HW{K* [r~l :r*l,liow( 1 :l) ) ' 


? 


end 




% % just to cneck tnat gajiunas do wliat tHey are auppoaed to. 




* rr - (l:dit) ; 




% Hba_aft6r( : , : ,rr) - state.H( : , : ,rr,bitl ) j * THIS IS Gaimiia_ab 




* HaE_aft6r( : , : ,rr) - state. H( :,:, rr, bit2 ) j * THIS IS Gaiiima_ac 





4.5 Degeneracy and Non- Frustration Condition for 
Generic Local Terms 



We choose to investigate chains of ci-dimensional quantum spins (qudits) with 2- 
local nearest-neighbor interactions Our first resuh is an analytic derivation of 
the necessary and sufficient conditions for such quantum systems to be unfrustrated. 
Second, we look at their ground state properties and find a range of parameters 
where we conjecture that these states are highly entangled and thus may be difficult 
to find computationally. We then corroborate this by a numerical investigation using 
a Matrix Product State (MPS) method. 

While the MPS formulation has been shown to work very well numerically for 
most one-dimensional particle systems, complexity theory issues seem to show there 
must be exceptions to this rule. Finding the ground-state energy of a one-dimensional 
qudit chain with d = 11 has been shown to be as hard as any problem in QMA 
[1011 1102 J. It is not believed that classical computers can efficiently solve problems 
in QMA. However, to our knowledge until now there have not been any concrete 
examples (except at phase transitions) for which MPS methods do not appear to work 
reasonably well. This research was undertaken to try to discover natural examples 
of Hamiltonians for which MPS cannot efficiently find or approximate the ground 
states. 



In Section 4.6 we show that the question of non-frustration for qudit chain Hamil- 
tonians with general nearest-neighbor interactions can be simplified to only Hamilto- 
nians that are sums of projector terms [103] . We then analytically show under what 

we 



conditions zero energy ground states for this system exist. Second, in Section 4.7 



-'^The rest of this chapter is based on [133) : however, I have corrected the erroneous assertion we 
made regarding choosing solutions. 
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show how to search for and approximate the ground states numerically and analyze 
the efficiency of finding the required MPS. We identify an interesting class of unfrus- 
trated qudit chain Hamiltonians, on which our MPS methods do not work well. Led 
by our numerical work, we conjecture that these ground states are highly entangled. 
Finally, we summarize our results and conclude with an outlook to further work in 
Section 

4.6 Generic Interactions 

We now choose to focus on a class of Hamiltonians whose local terms have generic 
eigenvectors. This implies 



is a random rank r projector acting on a d^-dimensional Hilbert space of two qudits, 
chosen by picking an orthonormal set of r random vectors (a different set for every 
qudit pair - we are not assuming translational invariance). 

We now find conditions governing the existence of zero energy ground states (from 
now on, called solutions in short). We do so by counting the number of solutions 
possible for a subset of the chain, and then adding another site and imposing the 
constraints given by the Hamiltonian. 

Suppose we have a set of D„ linearly independent solutions for the first n sites 
of the chain in the form 



r 
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similar to MPS, with = 1, . . . ,d and = 1, . . . , -D^; here and below all the re- 
peated indices are summed over. The F's satisfy the linear independence condition^ 

We now add one more site to the chain, impose the constraint Pn,n+i and look for 
the zero-energy ground states for n + 1 sites in the form 

j h,...,i„+i _ I ii,...,i„-i-ni„,ln\ -pin+ijn+l] / j o\ 

The unknown matrix ra„tal+t^^ must satisfy 



for all values of a„_i, a;„,+i and p, with \v^ vectors defined in (4.6.1 ). This results 
in a system of linear equations 

C pi„+l,[n+l] = Q (A a 5^ 

with Cp„„_i,i„+,„„ = (t;^^„+i|z„z„+i)ra,:!"i,a„ a matrix with dimensions rD^^i x 



If dDn > rDn-1 and the matrix C has rank rDn-i, the conditions given by Eq. 4.6.5 



are independent and we can construct dDn — rDn-i linearly independent T]^^^a^^\ 
corresponding to solutions for the n + 1 qudit chain (see the appendix for a proof 



^Note that this is not the standard MPS form, which also requires hnear independence in the 
other direction, i.e. 

In this case Sk would be the Schmidt rank for the partition of the qudits into (l,...,fc) and 
{k + l,...,n). 
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Figure 4.6.1: The existence of zero energy ground states for a qudit chain with 
(i-dimensional qudits and r projectors per pair. We highhght two notable cases: 
d = 2,r = 1 and d = 4,r = 4. 

of rank (C) = rD„_i). The freedom we have now is to use only a subset of them 
for constructing solutions (see the next chapter). Previously in our work [133] we 
asserted that the choice made implies 

Sn+i < dsn - rs„_i, (4.6.6) 

valid for all n. In order for the forgoing inequality to hold, one needs to prove that 
the possible dependences resulting from making choices do not affect the inequality 
(see next chapter). For example, excluding a subspace at a given step could break 
the full rankness of C at a later step. 

Let Dq = 1 and Di = d as the only constraint on Faj'*^ is linear independence. 
The recursion relation above at each gives linearly independent zero energy states, 
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where 



D„ = dDn-i-rD.n-2, (4.6.7) 



for all n with. The solution of this recursion relation is 



with f + g = d and fg = r. Hence, 



f-9 



^ d d /rf2 

There are three interesting regimes for r and d which yield different behaviors of Dn 



(Figure 4.6.1 ): 



1. r > ^ gives D„ = r 2 with cos0 = ttt- D„ becomes negative when 

n + 1 > I and thus no zero energy states can be constructed for a long chain 



if r > f . 



2. r = ^ results in D„ = (^)" (ra + 1), an exponential growth in n (except when 
d = 2, which gives linear growth). 



3. r < ^ implies / > | and f > g so for large n, Dn ~ /" ^1 — and the 
number of zero energy states grows exponentially. 
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Figure 4.7.1: (Color online) Ground state energy from imaginary time evolution vs. 
X for different ranks of the Hamiltonian. This is a plot for d = 4, and projector ranks 
of 2, 4, 6. Exact description would require x = = 2^°- 

4.7 Numerical investigation using Matrix Product 
States 



In this Section we use the methods described above to numerically search for the 



ground states of our class of random projector Hamiltonians (4.1.5). We probe 
the relations obtained in the previous Section, and see how well the energy coming 
from our small-x MPS imaginary time evolution converges to zero. The numerical 
technique we use is similar to Vidal's [971 1104] . We use imaginary time evolution to 
bring the system from a known state to its ground state: l^&grd) = hm^_>o< 



-ff-|*o>| 



Our experimentation with the parameters for a linear chain of length = 20 is 
shown in Figures 4.7.1[4.7.2 and 4.7.3 all the plots are on semi- log scale and the 
quantities being plotted are dimensionless. 

We see that for r < d the final energy converges to the zero energy ground state 
relatively fast with x ^ d^^"^- This can be seen in all the figures by the lowest curves 
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Figure 4.7.2: This is a plot for d = 5, and projector ranks of 4, 6, 8. Exact description 
would require x = 5^°- 



(marked by triangles). As can be seen the final energy obtained from imaginary time 
evolution tends toward zero with a steep slope, indicating that the ground state can 
be approximated efficiently with a small x iii MPS ansatz. 

The r > case, marked by squares, is shown by the top curves in all the 
figures. One sees that the final energy plateaus relatively fast in all three cases. This 
shows that the numerics have converged to a nonzero value and that increasing x 
will not yield a lower value of energy. Therefore, the numerical results suggest that 
there are no ground states with zero energy. 

In the previous section we analytically showed that when d < r < d'^/A there are 
many zero energy ground states. However, when we try to numerically find these 
states we see that the final energy converges to zero slowly. This is shown in all the 
Figures by the curves marked by circles. Out of these there are the critical cases, 
where r = ^. These correspond to the curves marked by closed circles in Figures 3 
and 5. The numerical investigation of the case d < r < d'^/4 is interesting because 
it suggests that for large number of spins finding the ground state with small x^ 
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Figure 4.7.3: The case of d = 6, and projector ranks: 5, 7, 9, 11. Exact description 
in general would require x = 6^°. 



tractable on a normal computer, is very hard. We interpret this as high amount of 
entanglement among the zero energy ground states and leave the analytical proof of 
this statement for a follow up paper. 



4.8 Summary 

We have investigated the no-frustration conditions for a system of qudits on a line 
with d states per site and random rank r local projector Hamiltonians acting between 
the nearest neighbor sites. We proved that there are no ground states with zero 
energy for r > ^ and sufficiently large A^. The system is not frustrated for r < ^. 
This second parameter region further splits into two. For d < r < many entangled 
zero energy ground states exist. On the other hand, for r < d we can also construct 



separable zero-energy ground states (see the next chapter and also Figure 4.6.1). 



We have verified the above numerically, in particular we have seen that when 
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d < r < d'^/A approximating the ground state energy (finding the ground states) is 
hard as the states seem to be highly entangled. 



4.9 Appendix 



We would like to say that for random \v) the rank of C is generically the maximum 



rank allowed, min(rD„_i, dDn)- The full rankness of C in Eq. 4.6.5 is not obviously 
true. In particular in the regime r < ^, to which we restrict ourselves from now 
on, Dn grows exponentially in n, while the number of parameters in the on 
which C depends only grows linearly. Thus C is far from a generic matrix of its size, 
but we now prove that its rank is indeed rD^^i. 

The argument used by Laumann et al [108J to prove their "geometrization theo- 
rem" also applies to our problem. It shows that for a chain of N qudits with random 
I'^fcfc+i)' ^-^^ f*-*^ ^ Hamiltonian H as in equations (2), (4) and (5) the number of 
zero-energy states, i.e. dim(ker(if)), is with probability one (which is what we mean 
by generic) equal to its minimum value. The calculation leading to the recursion 



relation Eq. 4.6.7 and its solution, shows that in the regime r < d'^/A this minimum 

is > D]^, since if the rank of the rD^^i x dDj. matrix C is ever less than rDk_i we 

can choose D^^i of them. Hence it is sufficient to find a single set of Iff^+i) for 

which dim(ker(iJ)) = Dn to prove that Dn is the generic value, i.e. that greater 

values occur with probability zero. This implies that the rank of each C is generically 

rD. _i, since otherwise at the first k where C had smaller rank we could construct 

more than Dk+i solutions for a chain of length k + 1. 

We construct with the property {v'l. ^_^i\ikik+i) = unless ^fc < | and 

ik+i > |. This can be done for r linearly independent if r < ^. We now assume 
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d is even; the modifications for d odd are obvious. We proceed by induction on n. 
Assume that in eachT^^l^ak with k < n, at runs from 1 to D^. From the definition 



of C (following Ec 
and so from Eq. 



|4.6.5 ) and the special choice of \v), Cpari^^,i„+ian = for in+i < 



4.6.5 



choose, for 1 < < 



Ta'jtal"^^^^ is unconstrained for in+i < |. This allows us to 

d 



Ea+ii+t^' = 1 when On+i = - 1) + in+i 



with 1 < a„ < 1 < in+i < I 

Ei"+ii"t^' = otherwise. 
As part of our induction, we assume that for 1 < < |-Dn_i, 

Ea„'L"\a„ = 1 when a„ = - 1) + in 

with 1 < < Dn-i, I <in< 



2 

EL"„'!"i,a„ = otherwise. 
Now we can show that the rows of Cpa„_,-|^,^„_^.J^a„ are linearly independent. Eor 
if' Zlp,a„_i 2/p,«n-iCpan-i,in+ian = ^^i all , this is true in particular for all 
in+i > d/2, an < |l^„-i, when it becomes Ep2/pan-i(^n,n.+iKri^n+i) = for all 

< |, in+i > |, and < -Dn-i- Since the are linearly independent, this is 
only true if |/pa„_i = for all p, Hence the rank of C is rDn-i and a^+i can 

take altogether dDn — rDn-i = Dn+i values, which is what we wanted to prove. 
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Chapter 5 

Entanglement of The Ground 
States 

The next natural question to ask is: how entangled are the ground states of qudit 
chains with generic local interaction? In this chapter I summarize my efforts in 
proving various results, which have not appeared elsewhere. As of now we have 
not succeeded in proving the main theorem (see the conjecture below); I include the 
partial results with the hope that they inspire future progress. In the previous chapter 
we showed that there are exponentially many ground states when d < r < (P/A (see 
Table). 



Paraiiietor range 


Nuiiil)er of grouiid states 


Frustrated? 


d<r< (P/A 


~ exp (n) 


No 


r > d'^/A 




Yes 



We will see that when r < d the ground states can be product states. How 
entangled are the ground states in the regime d < r < d^ /A when d > Al Simple 
theorems of algebraic geometry tell us that among the many ground states there are 
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highly entangled ground states. Can there be product states? If not, are all the 
ground states highly entangled (see the conjecture below)? By highly entangled we 
mean their Schmidt rank is exponentially lower bounded. 

There are two methods that can yield the bounds needed: 1) Choosing solutions 
as we build the ground states marching from one end of the chain and proving lower 
bounds on the number of solutions that need to be kept to build any state. 2) Relating 
the construction of the previous chapter to the Schmidt rank by working from both 
ends and matching solutions. 

5.1 Set Up 

5.1.1 Genericity of C^"] 

Due to generic local interactions, the complete set of eigenstates has parameter count 
equal to a polynomial in the number of spins, though the dimensionality is expo- 
nential. This implies that C^*^! is non-generic despite its entries being functions of 
polynomial random variables. However, as long as the variables are continuous, state- 
ments about the rank of C^"' can be made in a 'generic' sense. That is if the rank 
is full for some choice of random variables, the full rankness holds with probability 
one for random choice of those variables. 

5.1.2 Basis for Solutions and C^'^^ 

We want to understand the Kernel of C'"] as give by 

^pj„-i,i„+ia„ = ('^n,n+lNn^n+l)rLi'_i,a„ (5.1.1) 

Further, 
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Kn+l) = /3ll,i„+iKnWl) (5.1.2) 

with (3' s drawn randomly, say from a Gaussian distribution. C^"! solves: 



/^N ■p[ri+l],i„+i ^ Q 



(5.1.3) 



Inserting Eq,5.1.2 in Eq S.l.l (/3's are real) we get 



(5.1.4) 



Let -B be the r x d? matrix of /3's which we can write as d blocks of matrices [Bk] 
of size r X d put next to one another (not multiplied): 



B = [B1B2 ■ ■ ■ Bd 



(5.1.5) 



with, Bk 



d,k 



k = 1, . . . ,d. The matrix B has rank r; whereas 



Plk ' ' ' Pd,k 

rank (5^) = min (r, d) (both with probability one). 
In matrix notation. 



(5.1.6) 



Remark: In what follows, by the word generic we mean with probability one. 

Lemma 5. For r < d there are product states, i.e., there exists x = 1 solutions. 
Proof. We show that the states can be satisfied by product states. In this case. 
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taking the Schmidt rank to be one, we propose that F'"' = 7^ is a vector of size d 
corresponding to the physical index z„. The constraint matrix C'"' takes the form 
^pjn+i = -^lOn -^27^ ■ ■ ■ Bd'yt ' which is a r x d matrix. It clearly has full 
rank equal to r. The question becomes, can we build 7„+i to be a vector of size 
d7 The system of r < d equations with random coefficients can always be satisfied. 
Therefore indeed we can take r["+-^l = ■jn+i- n 

Corollary 1. When r > d, generically there are no product states. 

Proof. The only solutions is the trivial solution 7^4. 1 = 0. The problem is over 
specified; r constraints and d variables generically cannot be satisfied for r > d. □ 

Comment: When r > c? , (I ® is a full rank injective map. Pt"] has Dn inde- 
pendent columns and is a tall rectangular matrix and is therefore an injective map. 
Their composition (I (g) -Bi) F'"' is yet another full rank injective map, whose rank is 
the number of columns. It is however not clear that when we concatenate this map 



by (I (g) Bk) r'"! with 1 < k < d to obtain Eq. 5.1.6, the resulting matrix has rank 
equal to the number of rows. 

Comment: In our previous work we proved that there are many zero energy ground 
states when d < r < d'^/A and none when r > d'^/A. 

5.1.3 Constructing Solutions for d <r < d? /A 



In Eq. |5.1.9| taking Dq = 1 and Di = d, we build the solutions recursively by 
marching along the chain. The very first set of solutions r^j^'*^^, is a d x d diagonal 



matrix. From Eq. 5.1.6 we have C'^^ = (fiiFt^l ■ ■ ■ i^^^F'^l) = B. To find a basis for 
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the kernel we row reduce C'^' to put it in the row echelon form: 



C 



[1] 

echel 



1 A 



<P-r 



/3, 



(i2 



t t 

h ■■■ h^- 



Recall that F'^l is the set of solutions in the Kernel of C^^^. The null space of C'^J is 



span 



/ t ^ 



\ J 

the horizontal line depicts the partition of the ^'s from the unit vectors and does not 
signify any mathematical operation. 

Consequently, F^^ can be expressed as 



/ 



r[2] ^ 



t 

-h 



i 



This will in turn define the rd x d'^ matrix C^^ 
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whose Kernel defines Fl^l Continuing this way we arrive at 



t 

4 



, giving 



where 



t 


t 


t 


t 






■ ■ 7i 


7^: 


; 


i 


i 








Bd 





B. 



( t \ 

(k) 
Th 



l®Bu 



\ ^ I 

Therefore these constraints give 



/ t \ 

-k 
V ^ ) 



Dn+1 > dDn - rDn-i, 



132 



solutions with equality for C^"! being full rank. In our previous work we proved that 
CN is full rank: 

Dn+i = dDn-rDn-i. (5.1.9) 
The main theorem of this chapter remains to be proved: 



Conjecture 1. Generically all the ground states of the qudit chain with generic local 
interactions are all highly entangled in the regime d < r < d'^/4. 



In order to obtain entanglement bounds one can take at least two different ap- 
proaches. 



5.2 First Method: Choosing Solutions 

f h 

Let the n step be the first step in which we make a choice by throwing away w„ 
solutions and keeping s„ = Dn — Un solutions. For example if we throw out the 

second and the last column of Pl"] we obtain 
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t 

-h 
1 





t 



1 



t 

i 



\ 



(5.2.1) 



v 



1 





This in turn defines the rDn~i x dsn matrix C^^' 



(i,„®Ei)rN (i,„®52)r[^ 



whose dimension of Kernel dictates the number of independent solutions we can 
build, i.e., r["+il. We can keep making choices by excluding solutions to build some 
arbitrary state, i.e., Sk = Dk — Uk- 



Lemma 6. A lower bound on Sk is a lower bound on Xk- 



Proof. Xk can be obtained from Sk by applying the canonicality condition con- 
dition 2, in Theorem 1] to the r''^'; equivalently a further constraint of linear inde- 
pendence from right to left (see the previous chapter) needs to be imposed on the 
solutions. Since solutions with Xn are contained in solutions with s„ one can take 
the Xk = Sfc^, (lb denotes lower bound) to allow construction of any state. □ 
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5.2.1 Rank of 



Suppose at some step n we throw away m„ and keep s„, = Dn — Un of the solutions. 
In our previous work we proved that for Un = for all n, C^"''^ is full rank. The full 
rankness however does not generally hold if we throw away solutions. When a choice 



is made, as in the example shown in Eq. (5.2.1), the matrix of constraints becomes 



\ 



t 


t 


t 


t 






■ ■ /I 


(d) 


; 


i 


i 


i 


Bl 




B', 


B/ 



(5.2.2) 



Each solution that is excluded will reduce the columns of C'"' (i— fold and the super- 
scripts on -Bfc's remind us that there may be some columns missing as a consequence 
of having made a choice. It is easy to see how this comes about. Suppose we exclude 
the first column of r'"', then each of the top foremost -B^'s in Eq. (5.1.7) loses its 



first column. In particular if we throw away the last d solutions in F'"', one can see 
that C'"' above would have r rows all zeros as its last rows. 



5.2.2 Towards Entanglement Bounds 

Suppose we march along a semi-infinite line and choose solutions randomly then it 
is plausible to assume the recursion s„+i < dSn — ts„_i (i.e., no rank deficiency) to 
hold. In this case, the entanglement bounds can nicely be obtained 
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Lemma 7. Suppose < dsn — ts„_i holds everywhere, then Sn+i > qSn where 
g = ^ (l + -^jT—^), which for d > 4 and d < r < d'^/A implies exponentially large lower 
bound on Xn+i- 

Proof. We apply the set of inequahties, s„+2 < dsn+i — rSn and existence of solutions 
Sk > for all k. Positivity of s„+2 =^ Sn+i > ^s„. Let s„+i = ^s„ + u„ which implies 

Sn+2 < dUn- Wc nOW boUnd U„ using S„+3 < dSn+2-rSn+l Sn+3 < d'^Un-rSn+i = 



d^Un - - rUn = (rf^ - r) - Now Sn+3 > implies Un > My2-p)S 



Combining this with Sn+i = + we get 



-^Sn r = d 



which proves Sn = are lower bounds on solutions that grow exponentially with n 
for d > A. Using Lemma [6] we conclude Xn is exponentially lower bounded. □ 

Comment: The formulation above does not take into account the finiteness of 
the chain. Further it ignores the affect of possible rank deficiency due to choosing 
solutions by assuming s„+i < dsn — rsn-i holds at every step. In particular, for a 
finite chain of length A^, xa^-i ^ d whereas, q^~^ ~ exp (A^ — 1). 



5.3 Second Method: Matching Solutions 

One can approach the problem by marching along both ways and match solutions 
in between the two ends and ask what the lower bound on the Schmidt rank must 
be. It is notationally convenient to index the marching along from left and right 
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differently; let = m + n. Suppose we march along from left n sites building the 
r matrices up to and including T]^^^^},a„ and suppose we we march along from right 
m sites with Fq'^I^I^i • Recall that marching along from left (right) only imposed 
linear independence of the solutions from left (right). To build any state we need to 
match the solutions and apply the last constraint between sites n and m, 



^ C'pa,i_i,am-i;an,Om^an,am ~ 0) 

where = /3f„,,^radla„rL'::;[ai_i defines rDn-iDm-i constraints on 

DnD,n variables-entries of X. It can be checked that Dn+m = DnD^ — rDn-iD„i-i 
as expected. Let us make a crisp problem definition. Given, 



1. Eq. |5.3.1| has Dn+m = Djy solutions. Let us call the space of solutions S 

2. r'"] has Dn independent columns: If T]^,'}^l,anba„ = 0, Vi„,a„_i =^ ba„ = 

3. has Dm independent columns: If Ca^T'aJ^i-i = 0, V?™, am-i ^ c^^ = 

4. f3f ■ are generic; 
prove that 

{5} n {determinental variety of D„ x D^matrices with rank x ^ Xo} = 

with probability one for some large Xo- Using matrix notation, one wants to bound 
the rank of matrix X that satisfies 
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J2 (W,®5«„)r(z„)xr(z„) = o 



where each Bi^ is ar x d generic matrix as discussed above and F is dDn-i x Dn 
with Dn independent columns and each V (im) for a fixed im is Dm x -Dm-i- 

How do the entries of X depend on the focal constraints that were imposed at 
previous steps? The entries are polynomials of very high degre^- much larger than 
the size of the matrix. One can see this by using the following basis for building 
solutions 



X 



Y 



(adjX)-^ 
-/ 



n — 1 n m m — 1 

X 



The elements of X have entries that are homogeneous polynomial functions of local 
terms at the previous steps 



X., 



poly 

orD„-iU,n-i 



in \r"DoDi — D„-i (n \r^Dn-lD,^-Y 
\P\2) ■ ■ ■ \Pn-\,n) 



(5.3.2) 



(Wm,m-l) 



W 



21 



where we denote the matrix of the local terms between sites k and + 1 by Pk,k+i if 
we are marching from the left and by Wk,k+i if we are marching from the right. For 
example {f^^Y ^o^i '^"-i means that the polynomial dependence of entries of X on 

^The degrees of polynomials, shown below, were obtained by Jeffrey Goldstone. 
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the first set of local terms from left are homogeneous of degree r'^D^Di ■ ■ ■ -D„-i with 
respect to the elements of the random local terms between the first and the second 
sites. 

The conjecture can be answered if the lower bound on rank (X) is found. In 
particular, can one use the fact that the entries are very high powers of random 
parameters to argue in favor of an effective genericity for the matrix X. 

5.4 What Does Algebraic Geometry Buy You? 

Though the Hamiltonian does have generic local interactions, it is highly non-generic. 
In particular, it has to obey the restrictive tensor product structure. Moreover, 
the number of random parameters scale linearly with the system's size and are no 
more than r(PN, whereas the size of the Hamiltonian is d'^ . This prevents one from 
utilizing the techniques of Algebraic Geometry in any direct way to prove statements 
about the entanglement of all the ground states. However, one can use basic ideas 
of Algebraic Geometry [TOl [11] to show that among the many ground states there is 
at least a highly entangled state [TT]. We can consider the space of all the ground 
states and use the following proposition given in [11, prop. 10]. 

Proposition. Every bipartite system C^"^ (g) C^"' has a subspace S of Schmidt rank 
> X, O'^d of dimension dim{S) = {Dn—X + l)(-Dm^X + 1) 

It is straighforward to calculate x needed for the space of solutions with dim (S) = 
D]\f to have an intersection when a cut is made in the middle Dm = Dn = Dn. In 

2 

the previous chapter we obtained the functional dependence of on n 
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with 



For large N one gets x ~ /^/^ when r < ^ and x ~ A?" (I)'^^^ when r = ^. It is 
not surprising to see that there is at least one solutions with a high Schmidt rank; 
the conjecture requires a stronger result. 
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Chapter 6 



Examples of Quantum 2- SAT and 
Combinatorial Techniques 



Here I describe two examples of quantum 2-SAT (both of which are due to Sergey 
Bravyi [107J) on a chain of length 2n with three-dimensional {d = 3) and four- 
dimensional {d = 4) qudits, both of which have unique highly entangled ground 
states. In the d = 3 case, the ground state has a Schmidt rank that grows linearly 
with the number of sites and in d = 4 case, the ground state has Schmidt rank 
^ _ 2"+i _ 1, Below I show that the entanglement entropies for d = 3 and d = 4 cases 



are H = ^ log n + 0.645 and H = (-\/2 — l) n + | log2 n + ^ log2 ( j respectively. 



I provide numerical simulations to verify these formulas. In the next chapter we 
prove that the gap closes polynomially in the d = 3 case. As far as we know the 
technique we use for proving the gap in this case is new. The d = 3 example below 
gives the combinatorial background for the next chapter and d = 4 example does not 
appear elsewhere. 
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Figure 6.1.1: The states of the qudit. 

6.1 Quantum 2-SAT for balanced parentheses {d = 
3) 

We describe an example of a frustration-free 2-local Hamiltonian on a chain of n 
qutrits which has a unique highly entangled ground state ipQ. More precisely, if one 
cuts the chain in the middle, the Schmidt rank of ipo X ^ '^/2, while the entan- 
glement entropy is S ~ (1/2) log2n. The Hamiltonian is likely to have a polynomial 
spectral gap. 



Define a 3-letter alphabet (see Fig. 6.1.1) 

S = {/,r,0}. 



We shall identify / and r with the left and right brackets respectively, that is, / = [ 
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and r =] . Let us say that a string s G S" is balanced iff after removing all zeroes 
from s one gets a balanced sequence of brackets. More formally, for any s G E"- let 
Lj(s) and Ri{s) be the number of left and right brackets in s located in the interval 

Definition 2. A string s G is called balanced iff Li{s) > Ri{s) for all i — 
1, . . . ,n — 1 and Ln{s) — Rn{s). 

For example, for n = 2 there are only two balanced strings: 00 and Ir. For n — 3 
there are four balanced strings: 000, OZr, IrO, and lOr. For n — A there are nine 
balanced strings: 

0000 /OOr 
OOZr ZOrO 
OZOr ZrOO 
OZrO Z/rr 

We would like to construct a Hamiltonian whose unique ground state is the uniform 
superposition of all balanced strings. 

First we need to find a more local description of balanced strings. We shall say 
that a pair of strings s,t G E" are equivalent, s ^ t, if one can obtain s from i by a 
sequence of local moves 

00 i — > Ir, 01 < — > 10, Or < — > rO. (6.1.1) 
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applied to pairs of consecutive letters. For any integers p,q > let Up^q G be the 
string that has p leading r's and q tailing Ts, that is, 



p n—p—q q 



In particular, uqa = 0". 

Proposition. A string s G zs balanced iff it is equivalent to the all-zeros string, 

s ~ 0". Any string s G S" zs equivalent to one and only one string Up^q for some 
integers p,q > 0. 



Proof. Indeed, applying the local moves Eq. (6.1.1) we can make sure that s does 
not contain substrings Ir and 10 . . . Or. It means that if s contains at least one /, 
then all letters on the right of I are I or 0. Similarly, if s contains at least one r, 
then all letters on the left of r are r or 0. Since we can swap with any other letter 
by the local moves, s is equivalent to Up^q for some p, q. It remains to show that 
different strings not equivalent to each other. Indeed, suppose Up^q ~ Mp',q' 

such that p > p' . Then -Rp(s) — Lp{s^ < p' for any string s equivalent to Up>^qi. This 
is a contradiction unless p = p'. Similarly one shows that q = q'- See Figure □ 

It follows that the set of all strings S" is a disjoint union of the equivalence classes 
[up^q]. We shall now introduce projectors Q that "implement" the local moves Eq. 



(6.1.1) and a frustration-free Hamiltonian 



n-l 



HP'oP = J2Q.H^ (6.1.2) 



such that ground states of i/P^^P are 
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Entanglement of Motzkin States is due to 
the mutual information between halves 





















A 


B 






( ) ( ( 


( ) ) ) 

1 1 



Figure 6.1.2: A state in the ground state for d = 3 example. The high amount of 
entanglement is due to the high mutual information between the two halves of the 
chain. 
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Define quantum states \A),\B),\C)eC^<^C^ets 

\A) = |00) - |Zr), 

\B) = |0r) - |rO), 

and 

\c) = \oi)-m. 

Define a projector 

Q^li\A){A\ + \B){B\ + \C){C\) 

If a state I?/;) obeys Qjj+i = for all 1 < j < n — 1 then {s\ip) = (s'lip) for any 
pair of equivalent strings s,s'. Hence i^P^^P is indeed frustration free and its ground 
subspace is spanned by the states ipp^g. 

How can we exclude the unwanted ground states ipp^g with p and/or g 7^ ? 
The key observation is that the equivalence class i3„ = [mo,o] is the only class in which 
every string s satisfies si ^ r and Sn ^ I- Hence we can modify our Hamiltonian as 

H^HP'''P+\r){rU + \l){lU. 

Now if is a frustration- free Hamiltonian with the unique ground state \ipo)- 

Let us now show that the Schmidt rank of \ipn) grows linearly with n. Consider 
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a bipartition {1, . . . , n} = AB where A and B is the left and the right halves of the 
chain (we assume for simplicity that n is even). For any string s G let sa and 
sb be the restrictions of s onto A and B. If s is a balanced string, one must have 
Sa ~ Wo,p and ~ Up^Q for some < p < n/2, since each unbalanced left bracket in 
A must have a matching unbalanced right bracket in B. It follows that the Schmidt 
decomposition of ipo can be written (ignoring the normalization) as 

n/2 

IV^o) = ^ \i^O,p)A ® \i^p,o)B, 
p=0 



where the states ijjp^q are defined in Eq. (6.1.3). Therefore, the reduced density matrix 
of A has rank x = 1 + n/2. 

Numerical simulation shows that the entanglement entropy of A grows logarith- 
mically, S{A) ^ (1/2) log2(^), while the spectral gap of H decays polynomially, 
A ~ 1/n^ (see next chapter). 



6.1.1 Entanglement entropy 

To have a balanced string, we need to have an even number of slots available for r's 
and /'s after O's have been removed. First, recall that the Catalan numbers are 



Ck = ^— , A; = 0,1,- 




which among many other things, count the number of ways one can deposite and 
withdraw a dollar a day such that after k days one starts and ends with zero dollars 
without ever going negative. It is immediate to see that a string that has 2n — 2k 

147 



zeros has number of configurations 




, with (6.1.4) 



k = {!,■■■ ,n} 

This count (with zeros taken into account) is also known as Motzkin Numbers. Let 
^271 = Ylk=i '^ki ground state reads 




where i is any of the walks that start and end with height zero without ever going 
negative on 2n qutrits allowing only \, — > moves. 

In order to calculate the entanglement entropy, wc need to consider the correlation 
between the left and the right when a cut is made at some arbitrary bond. For 
simplicity we put the cut in the middle, i.e., at n. Let us define Cm to be the space 
of all states with m excess left parentheses on the first n qutrits; moreover TZm would 
be the space of all states with m excess right parentheses on the remaining n qutrits: 

Cm — {\Sm)l--n\Sm G {0, T, /}" with m eXCCSS Vs} 

Tim = {|s,„)„+i...2„|sm e {0,r,/}" with m excess r's} 
We want the states of the form 
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V^^ m=0 

Tracing over the right n qutrits we have 



m=o \|x)e£™ / \|^/>e7^„ 



Recalhng that we can have or / and r, the number of states in Cm denoted by 



Mm become 



n—m 



Mm,n = ( ) (# nn walks of height m on remaining (n — A;) qutri^1&.)l.5) 



k=0 

n 



k 



N = E<,- (6-1-6) 

m=0 

where nn means none-negative. The number of walks such that we start with zero 
and end with zero without going negative is given by the Catalan numbers. The 
same problem but ending with a positive height is given by a theorem, originally due 
to D. Andre (1887), the so called Ballot problem [1251 P- §]: 

Theorem. (D. Andre 1887) Let a,b be integers satisfying 1 < b < a. The number 
of lattice paths M (p) joining the origina O to the point (a, b) and not touching the 
diagonal x = y except at O is given by 

a + b \ ^ 

In other words, given a ballot at the end of which candidates P, Q obtain a, b votes 
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a—b 



respectively, the probability that P leads Q throughout the counting of votes is . 
First note that a = b + 1 gives the Catalan numbers. For us a + b = m — k + 1 



and a — b = m + 1. Using this in Eq. 6.1.5 and taking care of the parity 



M = + 1 [ ^ W ri-k + 1 

- i:(.:™)-{m-e:-r)} 

n!(m + l)V- — (6.1.9) 



.^^ {i + m + ly.il {n — 2i — m)\ 



(6.1.10) 



To check Eq. 6.1.7, we see that M„ = 1, corresponding to all left parentheses and 



fc=o ^ \ k J y lin- k) I k=o\k 

as expected. This count is also known as Motzkin triangles [126\ p. 4]. Consequently 
the Schmidt numbers become Pm = and the entanglement entropy is H (n) = 
— Ylm=oP"i^og2Pm- Before making approximations, using Maple one can express 

^ _n!(m + l) 2Fi{[-l{n-m) ,-l{n~m-l)],[m + 2],4:) 

r (m + 1) r (n — m + 1) 
where 2F1 denotes Hypergeometric function. 

For simplicity we have made the cut in the middle and will confine to this restric- 
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tion below. However, more generally, one can place the cut at site 1 < h < 2n — 1 
then 



mm(h,2n-h} / \ / 

h E EWE 



mm{h,2n—h) 

N 

m=0 \\x)£Cm I \|^/>67^„ 

where 



Cm = I |sm)i--h|'5m £ {0, r, with m excess /'s| 

Tim = ^\sm)h+i-2n\sm & {0, r, Z}^"^'' with m excess r'sj 

and Schmidt numbers would become pm = , where Mh,m is defined as 

above but of height m on h sites, similarly for M2n-h,m- The normalization being 



2n—h,m- 



Now we analyze the sum given by Eq. |6.1.9| carefully, 

Mn,m = (m + 1) ^ 



{i + m + l)\i\{n — 2i — m)\ 



First, though, let's do a little calculation. We will analyze a trinomial coefficient, 
where x + y + z = 0. 

We first use Stirling's formula to get 



+x '^+y ^ + zj \87T^{n + 3x){n + 3y){n + 3z) 



n/3+x / \ n/3+v / \ n/3+z 



n + 3xj \n + 3y J \n + 3z 



3 



n 
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Let's expand this by saying 



, n/3+a; 



n + 3x 



exp [ — [ — + X 



n f X 1 x'^\ 
exp -- 3- ---9- -3- 
6 \ n I ) n 

3x2 

exp I —X 



Thus, since a; + y + ^ = 0, we get 



n 



3v^ I n / 3^2^^^x 

exp — 3 



l+x f + y l + V 27r Y (n + 3x)(n + 3y)(n + 3z) V 2 n 



Now, we take the formula for M„, let m — a^/n and i = | + /3i/n, and use this 
approximation. The term inside the square root is approximately l/n^, so we make 
this substitution to get 



(m + 1) 



n+ 1 



n + 1 V^ + "T. + 1 i n — 2i — m 



27TO n 

3^3 
27rn3/2 

3v^ 
27rn3/2 



exp 



-3 



((a + + + (a + 2/3)') 



gn+l 



a exp 



(2a^ + 6a/3 + 6/3') 



3''+Vexp (-3a' - 9a/3 - 9/3') . 



We need to evaluate the sum of M„^^^j from i — to i — n. We approximate 
this by integrating over i. Since we have i = | + /3\/n, we get di = ^/ndji, Since the 
maximum is near i = |, we can turn this sum into an integral from — oo to oo. The 
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integral we need to evaluate is thus 



^^S'^+^a / exp (-3a^ - 9aP - 9/3^) ^/^d^ 
T+^a j exp {^9 (/3 - a/2f - -^a^^ y/^dp 
3"+^aexp • 



27rn 



2y/7rn 



This is maximized when 



^ / 3 2 1 rv 

— aexp — a ) = U, 
ao; V 4 



or a 



7273. 



Now, we need to figure out the entropy of the probability distribution proportional 
to M^^. Recalling that m = a^/n, and noticing that the normalization factor 
cancels, the entropy is 



/f({MU)«4E^exp(-^!^)log 

m=0 ^ / 



1 m 



7^ — exp , 
T n ^ \ 2 n 



3 



where 



m / 3 m 

— exp -- — 
n \ 2 n 

m=Q ^ 



We can approximate the sum with an integral. Since the integrand goes to rapidly, 
we can extend the upper limit of integration to 00, getting 



^({MU)--^/ - exp (---j log 



1 m 



3 



dm. 



where 



'- r ^ ^ /_3m2\ 



dm. 
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Now, we can replace m/ ^Jn by a again. This gives 



1 /"°° 3 



1 2 / 3 2 

—a exp \--a 



da. 



OO o 



T" = j exp(— -a^)(iQ;. 



where 



Here the extra log ^/n term comes from the fact that T' 
just a constant, so we can evaluate it to get 



^JnT" . The integral is 



^({^n,„^}) ~ ^logn-^ + 7+^(log2 + log7r-log3) nats (6.1.11) 
pa -logs^ + 0.64466547 bits. 



where 7 is Euler's constant. In Figure 6.1.3 we compare Eq. 6.1.11 with numerical 



evaluation of Eq. 6.1.7 
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Entanglement entropy S(n) vs. n for qutrits 




n 

Entanglement entropy S(n) vs. n for qutrits 




Figure 6.1.3: Entanglement entropy for the case of qutrits 



6.2 Quantum 2-SAT for Mirror Symmetric States 

(d = 4) 

We describe an example of a quantum 2-SAT on a chain of 2n four-dimensional 
{d = 4) qudits which has a unique satisfying state with the Schmidt rank growing 
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exponentially with n |1U7] . 



Let A and B be the left and the right halves of the chain containing n qudits 
each, that is, 

A = {1,2, . . . ,n} and 5 = {n + 1, n + 2, . . . , 2n}. 
Basis states of each qudit will be labeled using the alphabet 

S = {0,a,/3,7}. 

Let us first informally describe the idea behind the construction. The letters a, /3,7 
represent three 'particle types' while the state represents the 'vacuum'. The par- 
ticles can propagate freely through the vacuum, although they cannot pass through 
each other. Furthermore, the boundary between A and B is impenetrable for a and 
[3 particles, while 7 particles can propagate freely across the boundary. Let us first 
describe the role of a and /3 particles. The only place where a,/? can be created or 
annihilated is the boundary between A and B. Specifically, one can create/annihilate 
pairs aa or from the vacuum at qudits (n, n + 1). This will create a 'gas' of a 
and /3 particles such that the gas contained in A is the 'mirror image' of the gas 
contained in B if one ignores all zeroes. For example, aaOOO/3 : 0/3aOaO represents 
an admissible gas of particles for n = Q (here : represents the boundary between A 
and B). All possible admissible configurations of the gas will appear in superposition 
in the ground state. The mirror symmetry between A and B will be responsible for 
the exponentially large Schmidt rank. To maintain the mirror symmetry will will 
forbid pairs a/3 and fia on the boundary between A and B. This however is not 
sufficient by itself because it does not guarantee that A and B contain the same 
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number of particles. For example, a string 000 : a /3a does not have forbidden pairs 
on the boundary but it cannot propagate to any other string. Such strings could 
give rise to unwanted ground states with no entanglement between A and B. This 
is where 7 particles come to play. The rules for creation/annihilation of 7 particles 
are as follows: 

• Every a or (3 particle located in A can emit / absorb 7 particle on its right, that 
is, aO 0:7 and f30 -H- f3'j. 

• Every a or /5 particle located in B can emit/absorb 7 particle on its left, that 
is, Oa -ir^ 'ja and 0(3 ^ 7/3. 

• 7 particles are forbidden at qudits 1 and 2n 

Note that 7-particles cannot tell the difference between a and /3 particles, so they 
cannot maintain the mirror symmetry between A and B by themselves. The purpose 
of 7-particles is to ensure that the total number of a and /3 particles is the same in A 
and B. In the above example a string 000 : a^a can now propagate to a forbidden 
string: the leftmost a-particle in B emits 7-particle obtaining OO7 : a^a which now 
can propagate to a forbidden string 7OO : a^a since 7-particles can move freely 
through the vacuum. On the other hand, a balanced string like OOaO : OOaO cannot 
propagate to a forbidden string since 7-particles are confined to the interval between 
the two a particles. We will prove below (see Lemma |8| that the only strings that 
cannot propagate to a forbidden string are those representing a gas of a, (3 particles 
where the intervals between adjacent particles may be filled by O's and 7's and the 
gas contained in A is the mirror image of the gas contained in B (if one ignores all 
O's and all 7's). 

Let us now describe this construction more formally. We shall say that a pair of 
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strings s,t G E^" is equivalent, s ~ t, iff one can obtain s from t by a sequence of 
local moves listed below. These moves can be applied to some pairs of consecutive 
qudits {j,j + 1). We say that the pair is inside Aiffl<j<n — 1. We say that the 
pair is inside B iS n + 1 < j < 2n — 1. We say that the pair is on the boundary iff 
j = n. 

Move 1: Oa -H- aO, 0/3 -H- /30, O7 -H- 7O (inside A or inside B) 
Move 2: aO ■<->■ a^, /30 ■<->■ ^7 (inside A or on the boundary) 
Move 3: Oa 70;, 0/3 -H- 7/9 (inside B or on the boundary) 
Move 4: 00 <H- aa, 00 -H- O7 <H- 7O (on the boundary) 

In addition to these moves we shall impose several constraints: 

Constraint 1: Pairs a/3, /3a are forbidden on the boundary. 
Constraint 2: The first qudit of A is not 7. 
Constraint 3: The last qudit of B is not 7. 

Definition 3. A string s e E^" is called good iff all strings in the equivalence class 
of s obey Constraints 1,2,3. Otherwise a string s is called bad. 

Given a pair of strings s, s' e {a, /3}"^, we shall say that s' is the mirror image of s 
iff s'^ — Sm-i+i for alH = 1, . . . , m. For any string s e E^" let us denote sa and sb 
the restrictions of s onto A and B. 

Definition 4. A string s — {sa,sb) G E^" has mirror symmetry iff after removing 
all zeroes and all 7's the strings sa and sb become mirror images of each other. 

Lemma 8. A string is good iff it has mirror symmetry, the leftmost particle in A (if 
any) is a or /3, and the rightmost particle in B (if any) is a or /3. Any good string 
is equivalent to the all-zeroes string. 
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Proof. Let s = {sa, sb) G S^" be any good string. It is clear that the leftmost 
particle in A cannot be 7 since otherwise Move 1 would propagate 7 to the first 
qudit of A violating Constraint 2. By the same reason the rightmost particle in B 
cannot be 7. Let us show that s has mirror symmetry. If both sa,sb are all-zeroes 
strings we are done, so let us assume that sa contains at least one non-zero. Suppose 
sa contains at least one 7. Consider the left-most 7 and push it to the left until it 
gets absorbed by a or /3 (Move 2). Applying this to every 7-particle in sa we can 
assume that sa contains only 0, a, and f3. Applying Move 1 we can transform sa to 
the following canonical form 

Sa = xi,..., Xm) where x = {xi, Xm) e {a, P}"^ 

n—m 

for some m > 0. If Sb is all-zeroes string, we can apply Move 2 to the rightmost 
particle in A (which is Xm G {ct, f3}) to emit 7-particle, x^O XmJ- Propagating 
this 7-particle to the last qudit of B we violate Constraint 3. It shows that sb must 
contain at least one a or /3. Using the same arguments as above, we can apply 
Moves 1,3 to transform Sb into canonical form 

Sb = {Vk: ■ ■ • ,yi,0^^) where y = {yk, . . . G {o;,/?}^ 

n—k 

for some A; > 0. Using Move 4 and keeping in mind that s satisfies Constraint 1, 
we can consecutively annihilate all pairs Xiyi until we arrive at sa = or ss = 0. 
However the same arguments as above show that if sa = then the sting sb cannot 
contain a or /3, that is, = (and vice verse). Thus we proved that any good 
string is equivalent to the all-zeroes string. Since all moves used above preserve the 
mirror symmetry, we also proved that any good string has mirror symmetry. 
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Conversely, suppose a string s G S^" has mirror symmetry, the leftmost particle 
in A (if any) is a or /3, and the rightmost particle in B (if any) is a or /3. Let t be 
any string equivalent to s. We have to show that t obeys Constraints 1,2,3. Since all 
Moves 1,2,3 preserve the order of a, P particles in A and B, it clear that t satisfies 
Constraint 1. Since a or /3 particles located in A can emit 7 particles only on their 
right (Move 2), no 7 particle emitted in A can violate Constraint 2. A 7-particle 
emitted in B or on the boundary (Move 3) can violate Constraint 2 only if B contains 
at least one a or /3 particle while A does not. However this contradicts to the mirror 
symmetry. Thus t satisfies Constraint 2. The same arguments show that t satisfies 
Constraint 3. □ 

Consider a state \ip2n) G (C^)*^^" defined as the uniform superposition of all good 
strings, 

s is good 

Lemma 9. The state \1p2n) considered as a bipartite state shared by A and B has 
Schmidt rank 

Xn = 2"+^ - 1. (6.2.1) 

Proof. Indeed, the Schmidt basis of \1p2n) can be easily constructed using Lemma [sj 
Choose any integer m G [0,n] and any string x G {a,P}"^. Let \A{m,x)) be the 
uniform superposition of all strings s G S" of the form 

S = {Zq, Xi, Zi, X2, Z2, . . . , Xm, Zm), 

where Zq is a string of zeroes, and Zi, . . . , are arbitrary strings of zeroes and 7's. 
Any of the strings Zq, . . . , can be empty. Similarly, let \B{m,x)) be the uniform 
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superposition of all strings s G S" of the form 

where Zq is a string of zeroes, and Zi, . . . , Z^ are arbitrary strings of zeroes and 
7's. Any of the strings Zq, . . . , Zm can be empty. Using the characterization of good 
strings given by Lemma |8] we conclude that 

n 

\i^2n) = Y^ Yl I^Ka;)) ® |5(m,a;)) 

"1=0 x£{a,l3}'" 

is the Schmidt decomposition of \ip2n) (up to normalization of the Schmidt basis 



vectors). It immediately implies Eq. (6.2.1). □ 



Since the set of good strings is specified by 2-local moves and constraints, we 
can specify the state \ip2n) by 2-local projectors acting on nearest-neighbor qudits. 
Define auxiliary states 





|0a) 


- |aO) 


M/s) ~ 


m 


-1/30) 


M^) ~ 


IO7) 


-I7O) 


1-)- 




-|7) 


\Ca) ~ 


|00)- 


- \aa) 


IC/3) ~ 


|00)- 





We assume that all above states are normalized. Define a propagation Hamilto- 
nian responsible for 'implementing' Moves 1,2 for consecutive pairs of qudits 

inside A, namely 
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^prop,A _ \M^){M^\ + \Mp){Mp\ + \M^){M^\ 



Define a propagation Hamiltonian Hp^"^'^ responsible for 'implementing' Moves 1,3 
for consecutive pairs of qudits inside B, namely 



^prop,B ^ |M,)(M,| + |M^)(M^| + |M^)(M^| 



Define a propagation Hamiltonian H^^^' responsible for 'implementing' Moves 2,3,4 
on the boundary, namely 



^prop,AB ^ + + 

+ |-)(-| ® + |-)(-| ® + |Ca)(Ca| + m{Cp\ 

This Hamiltonian acts on the pair of qudits (n, n+1). Finally, define Hamiltonians 
imposing Constraints 1,2,3, namely. 



Here ifcon,^^ ^^^^ ^^^^ qudits (n, n+1). 
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Lemma 10. The state \il'2n) is the unique state annihilated by all the Hamiltonians 

j^prop,A jjprop,B jjprop,AB ^con,A jjcon,B Qfid JJcon,AB 

Proof. Indeed let H be Hamiltonian defined as the sum of all above Hamiltonians. It 
is clear that \ip2n) is annihilated by H. Since if is a stoquastic Hamiltonian, it suffices 
to consider ground states with real non-negative amplitudes. If has a positive 
amplitude on some string s, the propagation Hamiltonians ensure that has the 
same amplitude on any string equivalent to s. The Hamiltonians implementing the 
constraints then ensure that only good strings can appear in \-ip). Lemma [S] implies 
that there is only one equivalence class of good strings. Hence H has unique ground 
state \ip2n)- □ 

6.2.1 Entanglement entropy 

Recall that we impose three constraints on the states of the 2n qudits: 

1. Pairs a/? and Pa are forbidden at the boundary. 

2. The first qudit of A is not 7. 

3. The last qudit of B is not 7. 

The particles 7 can propagate freely through the the vacuum state given by states. 

We want to count the number of mirror symmetric states that obey the constraints 
2 and 3 (the first constraint is implied by mirror symmetry). First let us ask: how 
many strings can there be in A alone? Well out of the 4" possible strings the ones 
that violate constraint 2 need to be excluded (we are not worrying about B yet). 
The complete list of the excluded states is (each row represents a forbidden string in 
Sa) 
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7 


# 


# 




# 


4"-i 







7 


# 


# 


# 












7 


# 


# 


4"-3 


(6.2.2) 














7 


1 





where 7^ G {a, /3, 7, 0} denotes any state and the counts are written to the right. 
Therefore, the total number of possible strings in sa is, 



Number of allowed strings in = 4" - T,'^Zl^''. 

Definition, (m— dense string) A string of size m is m— dense if it has no O's or 7's. 

We wish to find the number of symmetric states where every string on A is 
m— dense (i.e., there are m qudits on A that are not 7 or 0). The number of mirror 
symmetric states becomes 



(Count of n — m particles of type 0, 7 in ^)^ 



(Count of n — m particles of type 0, 7 in ^4)^ 



m 



( \ ( \ 

m 



V 1 ; 




m 



m 



m 



where as before <| + IM^ number of ways that a 

and P particles can be positioned in m slots. To find the number of allowed m— dense 
strings in A we first count all possible (unconstrained) ways of putting n — m of 
or 7 particles and m of a and /3 particles on the n qudits. We then subtract from it 
the forbidden states. The total number of ways one can have an m— dense chain is 
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(without imposing the constraints) 



(Number of ways to choose n — m slots for and 7) x 
(Number of ways to place 0, 7 on the n — m qudits) x 
(Number of ways to place a, /3 on the remaining m) 



Mathematically 



Number of unconstraint m — dense chains = | | 2"^ "^2™ 

m 



The number of states that we need to exclude in A are 



71 # # # # ••• #: 1. I ^ I {2-™-i}2™ 

n — m — 1 

72 # # # ••• #: 1. I """^ I {2— -2} 2- 

n — m — 2 



7fc # # ... #: 1. I ^ I {2-™-'^} 2^ 

n — m — /c 



7„_„ ... #: 1. I ^ 1 I I 2 



m 
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where, l<k<n — mis the first k — 1 zeros followed by a 7; # is means it can 
be any state as long as we have a total of n — m of O's and 7's and m of a, p. The 
number of states that need to be excluded are therefore 



n—m I 7 

El ^ I L' 

I 2"" . (6.2.6) 

fc=i \ m 

Comment: Mathematica erroneously expresses the foregoing equation in terms of a 
Hypergeometric function, that has poles for integer n. 

In summary the number of allowed m— dense state on the n qudits are 



)n—m I J 
2" - ^ I 2"" • (6-2.7) 

k=i \ rn 



The ground states are 



m=0xe{a,/3}'" 

The number of ways that 0, 7 can be put on n qudits are 



M„,„ = 2"— <J I - 5^ 1 ~ I 2-^= ^ . (6.2.8) 

m j k=i\rn 

Consequently, the Schmidt numbers are prn,n = ]^'" where Mm,n is given by Eq. 



6.2.8 



and = X]m=o^™^mn is normalization constant. The entanglement 



entropy becomes 
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H {{Pm,n}) = - ^ 2™p„,„ \og2Pm,n- 

Let us rewrite Eq. |6.2.8| as 



(6.2.9) 



m=0 




where r = \ {^~~)- Next we use Stirling's approximation n\ ~ (n/e)" a/27™ , 
similarly for m! and (n — m)!, to obtain (below all the logrithms are in base 2 unless 
stated otherwise) 



f{n,m) 



2nm exp [/ (m, n)] 
71 {n -m) (n + m)^ 
2 [{n — m) log (2) + n log n — m log m 



m) log 



(6.2.12) 

m)] 



Let m = an, giving 



/(«) 



2a; exp [n/ (a)] 
7rn(l-a) (1 + a)^ 

2[(1 - a) log (2) + logn-alogan- (1 - a) log n (1 - a)] (6.2. 13) 



2[(a-l)log(l 



a] 



a log a + (1 — a) log (2)] 



(6.2.14) 
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N 

9 (n, m) 



E2" 

m=0 



m,n 



E 



2nmexp [g {n, m)] 
^ TT (n — m) (n + m)^ 



/ {n, m) + mlog 2. 



We can approximate this sum with an integral over a 



N ~ 



^ 2a cxp [ng («)] 
go; 



/o TT (1 — a) (1 + a 
5((q;) = / (a) + a log 2 



2 ' 



(6.2.15) 



Oi "I 
(a - 1) log (1 - a) - a log q; - - log 2 + log (2)J , 



note that the factor of n cancelled because of change of variables from m to a. In 
anticipation of the steepest descent approximation to the entanglement entropy, we 
evaluate 



9 = 



_ dg {a) 



da 
d'gja) 



log(l-a) - logo; - ^^og2 



= 2 



a{a — 1) 



g' — Q ^ aQ — — "n/S — 1 and g" (cto) = — (3 + 2-\/2), which implies ckq 



is a maximum. Let us proceed in calculating the entanglement entropy given by 
H (m, n) = — Xlm=o log2 ^ by first approximating the sum with an integral 

over a and then performing the steepest descent approximation (in nats) 
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/ r , ^ / MX 1 , 2aexp\nq(a)] , M^(a,n) 

H{{M{a,n)}) ~ -- da——''^ V , ) 



(1 -«)(! + af 



_ M^(Qio,n) / 1 r^^^ 2Qiexp[ng(Qi)] \ 

- TV \7VX ''7r(l-a)(l + a)'J 

{ao, n) 



^ -log ^ , 
by the definition of N . It remains to calculate — log ^ 



log 



{ao,n) 



N 



= -log-f 



h exp [nf (ao)] 



log 



Jg do; exp [ngi (a)] 

^ exp [n/ (gp)] 

f^^ da exp {n [g {ao) + \g" (ao) (a - "o)^] } 



■ log ^ ^ exp [nf (ao) - (ao)] y 



= — log < exp [— nao log 2] 



(3 + 2^^2) 



n 



> . 



Expanding this we obtain 



H{{M{a,n)]) ~ (^2 - 1^ nlog2 + ^logn+ ^ log ^ g^^ j ^^^^ 
= (v^-l)n+ilog,n+ilog,(^^^j bits. 
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Checking normalization: N(n) vs. n 

1 0^^° I : : 




50 100 150 200 



Figure 6.2.1: Normalization as a function of n. 



Entanglement Entropy vs. n 




n 



Figure 6.2.2: Entropy H ({M^}) vs. n for d = 4 case. We include the figure on left 
to demonstrate the closeness of the approximation. 
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Chapter 7 

Criticality Without Frustration for 
Quantum Spin-1 Chains 

In the previous chapter we showed two examples of FF qudit chains with high en- 
tanglement and introduced the mathematical techniques needed for calculating their 
entanglement entropies. Here we elaborate on the d = 3 model- balanced parenthesis 
model. While FF spin-1/2 chains are known to have unentangled ground states, the 
case s=l remains less explored. We propose the first example of a FF translation- 
invariant spin-1 chain that has a unique highly entangled ground state and exhibits 
some signatures of a critical behavior. The rest of this chapter also appears in |1U9] . 

7.1 Motivation 

The presence of long-range entanglement in the ground states of critical spin chains 
with only short-range interactions is one of the most fascinating discoveries in the 
theory of quantum phase transitions [121 IllOl 1111] . It can be quantified by the 
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scaling law S{L) ~ logL, where S{L) is the entanglement entropy of a block of L 
spins. In contrast, non-critical spin chains characterized by a non- vanishing energy 
gap obey an area law |112^lll3^lll4j asserting that S{L) has a constant upper bound 
independent of L. 

One can ask how stable is the long-range ground state entanglement against 
small variations of Hamiltonian parameters? The scaling theory predicts |110l 1115] 
that a chain whose Hamiltonian is controlled by some parameter g follows the law 
S{L) ~ logL only if L does not exceed the correlation length C, ^ \g — gd'", where 
z/ > is the critical exponent and gc is the critical point. For larger L the entropy 
S{L) saturates at a constant value. Hence achieving the scaling S{L) ~ log L requires 
fine-tuning of the parameter g with precision scaling polynomially with 1/L posing 
a serious experimental challenge. 

The stringent precision requirement described above can be partially avoided for 
spin chains described by frustration-free Hamiltonians. Well-known (non-critical) 
examples of such Hamiltonians are the Heisenberg ferromagnetic chain |116j . the 
AKLT model |117] . and parent Hamiltonians of matrix product states [HI 1118] . 
More generally, we consider Hamiltonians of a form H = gjUjj+i, where ^j.j+i is 
a projector acting on spins j,j + 1 and gj > are some coefficients. The Hamiltonian 
is called frustration-free (FF) if the projectors ^jj+i have a common zero eigenvector 
ip- Such zero eigenvectors ip span the ground subspace of H. Clearly, the ground 
subspace does not depend on the coefficients gj as long as they remain positive. 
This inherent stability against variations of the Hamiltonian parameters motivates a 
question of whether FF Hamiltonians can describe critical spin chains. 

In this Letter we propose a toy model describing a FF translation- invariant spin-1 
chain with open boundary conditions that has a unique ground state with a loga- 
rithmic scaling of entanglement entropy and a polynomial energy gap. Thus our FF 
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model reproduces some of the main signatures of critical spin chains. In contrast, 
it was recently shown by Chen et al |119j that any FF spin-1/2 chain has an un- 
entangled ground state. Our work may also offer valuable insights for the problem 
of realizing long-range entanglement in open quantum systems with an engineered 
dissipation. Indeed, it was shown by Kraus et al |12U] and Verstraete et al |121] that 
the ground state of a FF Hamiltonian can be represented as a unique steady state of 
a dissipative process described by the Lindblad equation with local quantum jump 
operators. A proposal for realizing such dissipative processes in cold atom systems 
has been made by Diehl et al |122j . 

Main results. We begin by describing the ground state of our model. The three 
basis states of a single spin will be identified with a left bracket I = [, right bracket 
r =], and an empty space represented by 0. Hence a state of a single spin can be 
written as a|0) + + 7|r) for some complex coefficients a,/3,7. For a chain of 
n spins, basis states |s) correspond to strings s G {0,/,r}". A string s is called a 
Motzkin path \123\ ?] iff (i) any initial segment of s contains at least as many /'s as 
r's, and (ii) the total number of /'s is equal to the total number of r's. For example, 
a string lllrOrlOrr is a Motzkin path while lOlrrrllr is not since its initial segment 
lOlrrr has more r's than Vs. By ignoring all O's one can view Motzkin paths as 
balanced strings of left and right brackets. We shall be interested in the Motzkin 
state \Ain) which is the uniform superposition of all Motzkin paths of length n. For 
example, jA^a) ~ |00) + \lr), {Ms) ~ |000) + \lrO) + \lOr) + \Olr), and 

IMa) ~ |0000) + |00/r) + \OlOr) + |/00r) 

+ |0/rO) + |/OrO) + |/rOO) + \llrr) + \lrlr). 

Let us first ask how entangled is the Motzkin state. For a contiguous block of spins 
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A, let Pa = T^t^ji^a \-M.n){.Mn\ be the reduced density matrix of A. Two important 
measures of entanglement are the Schmidt rank x{A) equal to the number of non- 
zero eigenvalues of pA, and the entanglement entropy S{A) = — Tr logs p^i- We 
will choose A as the left half of the chain, A = {1, . . . , n/2}. We show that 

xiA) = l + n/2 and S{A) = ^log^n + (7.1.1) 

where lim„_>oo Cn = 0.14(5). The linear scaling of the Schmidt rank stems from the 
presence of locally unmatched left brackets in A whose matching right brackets belong 
to the complementary region B = [l,n]\A. The number of the locally unmatched 
brackets m can vary from to n/2 and must be the same in A and B leading to 
long-range entanglement between the two halves of the chain. 

Although the definition of Motzkin paths may seem very non-local, we will show 
that the state \A4n) can be specified by imposing local constraints on nearest- 
neighbor spins. Let 11 be a projector onto the three-dimensional subspace of ® 
spanned by states |0/) — |/0), |0r) — |rO), and |00) — \lr). Our main result is the 
following. 

Theorem 1. The Motzkin state \A4n) is a unique ground state with zero energy of 
a frustration-free Hamiltonian 

i/=|r)(r|i + |0(/|n + $^n,,,+i, (7.1.2) 

where subscripts indicate spins acted upon by a projector. The spectral ga'^ of H 
scales polynomially with 1/n. 

^Here and below the spectral gap of a Hamiltonian means the difference between the smallest 
and the second smallest eigenvalue. 
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The theorem remains true if H is modified by introducing arbitrary weights gj > 1 



for every projector in Eq. (7.1.2). A polynomial lower bound on the spectral gap 
of H is, by far, the most difficult part of Theorem [Tj Our proof consists of several 
steps. First, we use a perturbation theory to relate the spectrum of H to the one of 
an effective Hamiltonian H^s acting on Dyck paths — balanced strings of left and 
right brackets!^ This step involves successive applications of the Projection Lemma 
due to Kempe et al |129j . Secondly, we map H^s to a stochastic matrix P describing 
a random walk on Dyck paths in which transitions correspond to insertions/removals 
of consecutive Ir pairs. The key step of the proof is to show that the random walk 
on Dyck paths is rapidly mixing. Our method of proving the desired rapid mixing 
property employs the polyhedral description of matchings in bipartite graphs |13U] . 
This method appears to be new and might be interesting on its own right. Exact 
diagonalization performed for short chains suggests that the spectral gap of H scales 



as A ~ see Fig. 7.1.1 Our proof gives an upper bound A = 0{n ^/^) and a 



lower bound A = Q{n for some c ^ 1. 



Previous work. Examples of spin chain Hamiltonians with highly entan- 
gled ground states have been constructed by Gottesman and Hastings |131] . and 
Irani |132] for local dimension d = 9 and d = 21 respectively (here and below 
d = 2s + 1) . These models exhibit a linear scaling of the entropy S{L) for some 
blocks of spins while the spectral gap is polynomial in 1/n. The model found in |132] 
is FF and translation-invariant. Ref. |133] focused on 'generic' spin chains with a 
Hamiltonian H = Hjj+i where the projectors are chosen randomly with 

^One can regard Dyck paths as a special case of Motzkin paths in which no '0' symbols are 
allowed. 
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Figure 7.1.1: The spectral gap A of the Hamiltonian H defined in Eq. (7.1.2) for 
3 < < 13 obtained by the exact diagonalization. The dashed line shows a linear 
fit log A = 0.68 — 2.91 log n. Our numerics suggests that the first excited state of H 
belongs to the subspace spanned by strings with exactly one unmatched bracket. 



a fixed rank r |^ The authors of |133j identified three important regimes: (i) frus- 
trated chains, r > (ii) FF chains, d < r < d'^/A, and (iii) FF chains with 
product ground states, r < d. It was conjectured in |133j that generic FF chains in 
the regime d < r < d'^/4 have only highly entangled ground states with probability 
one. This regime however requires local dimension d > A. The new model based on 
Motzkin paths corresponds to the case d = r = 3 (ignoring the boundary terms) and 
thus it can be frustrated by arbitrarily small deformations of the projectors making 
them generic. In addition, results of |133] imply that examples of FF spin-1 chains 
with highly entangled ground states have measure zero in the parameter space. The 
question of whether matrix product states specified by FF parent Hamiltonians can 
exhibit quantum phase transitions has been studied by Wolf et al |135J . However, 
the models studied in jl35j have bounded entanglement entropy, S{L) = 0(1). 

Hamiltonian. Let us now construct a FF Hamiltonian H whose unique ground 

■^Though the results of Ref. |133l are applicable to more general Hamiltonians, the convenient 
restriction to random projectors is sufficient for addressing the degeneracy and frustration condition. 
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state is \A4n)- First we need to find a more local description of Motzkin paths. Let 
S = {0, /, r}. We will say that a pair of strings s, t G S" is equivalent, s ~ t, if s can 
be obtained from t by a sequence of local moves 



00 ^ Ir, 01 O 10, Or ^ rO. (7.1.3) 

These moves can be applied to any consecutive pair of letters. For any integers 
p,q > such that p + q < n define a string 

p n—p—q q 

Lemma 1. Any string s &Y7^ is equivalent to one and only one string Cp^q. A string 
s eT/^ is a Motzkin path iff it is equivalent to the all-zeros string, s ~ co,o- 



Proof. Indeed, applying the local moves Eq. (7.1.3) one can make sure that s does 
not contain substrings Zr or ZO . . . Or. If this is the case and s contains at least one /, 
then all letters to the right of / are I or 0. Similarly, if s contains at least one r, then 
all letters to the left of r are r or 0. Since we can swap with any other letter by 
the local moves, s is equivalent to Cp^q for some p, q. It remains to show that different 
strings Cp^q are not equivalent to each other. Let Lj{s) and Rj{s) be the number of 
/'s and r's among the first j letters of s. Suppose Cp^q ~ Cpi^q' and p > p'. Then 
Rp{s) — Lp{s) < p' for any string s equivalent to Cp'^qi. This is a contradiction unless 
p = p'. Similarly one shows that q = q'. □ 

The lemma shows that the set of all strings can be partitioned into equivalence 
classes Cp^q, such that Cp^q includes all strings equivalent to Cp^q. In other words, s G 
Cp^q iff s has p unmacthed right brackets and q unmatched left brackets. Accordingly, 
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the set of Motzkin paths coincides with the equivalence class Co,o- 

Let us now define projectors 'implementing' the local moves in Eq. (7.1.3). Define 
normalized states 



10) ~ |00) - |/r), IV') ~ |0/) - I/O), It/;") ~ |0r) - |rO) 

and a projector 11 = |0)(0| + \ip^){ip^ \ + \il)'^){ip'^\. Application of 11 to a pair of spins 
J, J + 1 will be denoted Hj- j+i. If some state ip is annihilated by every projector 
Iljj+i, it must have the same amplitude on any pair of equivalent strings, that is, 
{s\ip) = (^IV) whenever s ^ t. It follows that a Hamiltonian = "^jZl^jj+i 
is FF and the ground subspace of is spanned by pairwise orthogonal states 
\Cp^q), where \Cp^g) is the uniform superposition of all strings in Cp^g. The desired 
Motzkin state |A^„) = |Co,o) is thus a ground state of iJ^. (It is worth mentioning 
that not all states \Cp^q) are highly entangled. For example, |C„,o) = |r)®" is a 
product state.) How can we exclude the unwanted ground states \Cp^g) with p 
or g 7^ 0? We note that Cq.o is the only class in which strings never start from r 
and never end with /. Hence a modified Hamiltonian H = |r)(r|i + |/)(/|n + that 
penalizes strings starting from r or ending with / has a unique ground state |Co,o)- 
This proves the first part of Theorem [TJ We can also consider weighted Hamiltonians 
H^ig) = Y.'^Zl 9j'^j,j+i a.nAH{g) = go\r){r\i+gn\l){l\n+H^{g), where go, . . . , gn > I 
are arbitrary coefficients. One can easily check that the ground state of H{g) does 
not depend on g and H{g) > H. It implies that the spectral gap of H{g) is lower 
bounded by the one of H. 

Entanglement entropy. We can now construct the Schmidt decomposition of 
the Motzkin state. Let A = {1, . . . ,n/2} and B = {n/2 + 1, . . . ,n} be the two halves 
of the chain (we assume that n is even). For any string s G let sa_ and sb be 
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the restrictions of s onto A and B. We claim that s is a Motzkin path iff sa ~ Co,m 
and ~ Cm,o for some < m < n/2. Indeed, sa (sb) cannot have unmatched right 
(left) brackets, while each unmatched left bracket in sa must be matched with some 
unmatched right bracket in sb- It follows that the Schmidt decomposition of \Mn) 
can be written as 

n/2 

\Mn) = ^y/p^\Co 
m=0 

where \Cp,q) is the normalized uniform superposition of all strings in Cp^q and pm are 
the Schmidt coefficients defined by 

l<-'0,o(^)| 

Here we added an explicit dependence of the classes Cp^q on n. For large n and m one 
can use an approximation pm ~ m^exp (— 3m^/n), see the Supplementary Material 
for the proof. Note that pm achieves its maximum at m* ^ a/?t./3. Approximating 
the sum '^^Pm = 1 by an integral over a = m/y/n one gets Pm ~ n~^/'^qa{m), where 
Qa is a normalized pdf defined as 







4 ■ 33/2 ■ 



It gives 

POO 

S{A) = - y^Pm log2 Pm ^ log a/^ - / da logg Qo 

^ Jo 



Evaluating the integral over a yields Eq. (17.1. 1|). The approximation p^ ^ n ^^'^qa{m) 



also implies that max^Pm = 0{n ^Z^). This bound will be used below in our spectral 
gap analysis. We conjecture that one can achieve a power law scaling of S{A) in 



Eq. (7.1.1) by introducing two types of brackets, say / = [, r =], /' = {, and r' =}, 
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such that bracket pairs Ir and I'r' are created from the 'vacuum' 00 in a maximally 



entangled state (|Zr) + |Z'r'))/V2. The local moves Eq. (7.1.3) must be modified as 
Oa; ^ xO, where x can be either of /, r, r', and 00 ^ (/r + I'r') / ^J2. We expect the 
modified model with two types of brackets to obey a scaling Si^A) ~ ^Jn^ while its 
gap will remain lower bounded by an inverse polynomial. 

Spectral gap: upper bound. Let A2 > be the smallest non-zero eigenvalue 



of the Hamiltonian defined in Eq. (7.1.2). We shall use the fact that the ground 



state \Ain) is highly entangled to prove an upper bound A2 < 0{n ^/^). Fix any 
k G [0,n/2] and define a 'twisted' version of the ground state: 

n/2 
Tn=0 

where 6m = for < m < k and 6m = tt otherwise. Note that |0) and |A^„) have 
the same reduced density matrices on A and B. Hence {(j)\H\(f)) = {(j)\Ilcut\4') , where 
Ucut = n„/2,n/2+i- Since max„j9„ = 0(n~^/^) and Y^mPm = 1, there must exist 
k G [0, n/2] such that X]o<m<fcPm = 1/2 ± e for some e = 0(n~^/^). This choice of k 
ensures that (0|A^„) = XlmPme*^™ < 2e, that is (p is almost orthogonal to the ground 
state. Define a normalized state |0) ~ |0) — {Ain\(p) ■ |A^n)- Then (0|A^„) = and 
{4>\H\4>) = {4>\Il,ut\4>) < (0|n,„t|0)+O(e). The difference {<f)\U,^t\<l>) - {Mn\Ilcut\Mn) 
gets contributions only from the terms m = k,k ±1 in the Schmidt decomposition, 
since Hcut can change the number of unmatched brackets in A and B at most by one. 
Since {Mn\Tlcut\Mn) = 0, we get 

imcuS) < 0{pk+Pk-i+Pk+i) = 0{n-'/'). 
We arrive at {(p\H\(f)) = 0{n-'^/'^). Therefore A2 is at most 0(n"^/2^. 
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Spectral gap: lower bound. It remains to prove a lower bound A2 > n^'^^^\ Let 
Tip^q be the subspace spanned by strings s G Cp^q and Hm = Tiofi be the Motzkin 
space spanned by Motzkin paths. Note that H preserves any subspace V-p^q and the 
unique ground state of H belongs to Hm- Therefore it suffices to derive a lower 
bound n'^^^^'^ for two quantities: (i) the gap of H inside the Motzkin space "Hm, and 
(ii) the ground state energy of H inside any 'unbalanced' subspace Tip^g with p 7^ 
or g 7^ 0. Below we shall focus on part (i) since it allows us to introduce all essential 
ideas. The proof of part (ii) can be found in the Supplementary Material. 

Recall that a string s G {/jr}^™ is called a Dyck path iff any initial segment of s 
contains at least as many /'s as r's, and the total number of /'s is equal to the total 
number of r's. For example, Dyck paths of length 6 are lUrrr, Urlrr, llrrlr, Irlrlr, 
and Irllrr. The proof of part (i) consists of the following steps: 
Step 1. Map the original Hamiltonian H acting on Motzkin paths to an effective 
Hamiltonian ifeff acting on Dyck paths using perturbation theory. 
Step 2. Map ifeff to a stochastic matrix P describing a random walk on Dyck paths 
in which transitions correspond to insertions or removals of consecutive Ir pairs. 
Step 3. Bound the spectral gap of P using the canonical paths method |136[ ?]. 

To construct a good family of canonical paths in Step 3 we will organize Dyck 
paths into a rooted tree in which level-m nodes represent Dyck paths of length 2m, 
edges correspond to insertion of Ir pairs, and each node has at most four children. 
Existence of such a tree will be proved using the fractional matching method jl30] . 

Let T)m be the set of Dyck paths of length 2m, T) be the union of all "Dm with 
2m < n, and M.n be the set of Motzkin paths of length n. Define a Dyck space 
T-Ld whose basis vectors are Dyck paths s E T). Given a Motzkin path u with 2m 
brackets, let Dyck('u) G T>m be the Dyck path obtained from u by removing zeros. 
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We shall use an embedding V : T-Lo ^ Hm defined as 

v\s) = ^= Yl 1^)' sevnVm. 

V ^^"-^ Dyck(M)=s 

One can easily check that V'^V = I, that is, V is an isometry. For any Hamiltonian 
if, let \2{H) be the second smallest eigenvalue of H. 



Step 1. The restriction of the Hamiltonian Eq. (7.1.2) onto the Motzkin space Hm 
can be written as H = Hmove + -f^mt, where Hmove describes freely moving left and 
right brackets, while Hint is an 'interaction term' responsible for pairs creation. More 
formally, Hmove = Z]j=i n^/+i and Hm = Z]j=i n*^/+i5 where n™°'"' projects onto 
the subspace spanned by |0/) — |/0) and |0r) — |rO), while 11™* projects onto the 
state 1 00) — \lr). Note that the boundary terms in Eq. (7.1.2) vanish on T-Lm- We 
shall treat Hint as a small perturbation of Hmove- To this end define a modified FF 
Hamiltonian H^ = Hmove + (.Hint-, where < e < 1 will be chosen later. One can 
easily check that \M.n) is the unique ground state of H^ and \2{H) > \2{H^) (use 
the operator inequality H > H^). Note that Hmovei' = iff is symmetric under 
the moves 0/ 10 and Or -H- rO. It follows that the ground subspace of Hmove is 
spanned by states V \ s) with s E V. To compute the spectrum of Hmove, we can 
ignore the difference between Vs and r's since Hmove is only capable of swapping 
zeros with non-zero letters. It follows that the spectrum of Hmove must coincide with 
the spectrum of the Heisenberg ferromagnetic spin-1/2 chain, that is, H'"°''^ can be 
replaced by the projector onto the singlet state |01) — |10), where 1 represents either 
/ or r. Using the exact formula for the spectral gap of the Heisenberg chain found by 
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Koma and Nachtergaele we arrive at X2{Hmove) = 1 — cos (^) = Q{n ^). Let 

be the first-order effective Hamiltonian acting on the Dyck space I-Ld- Applying the 
Projection Lemma of |129j to the orthogonal complement of \A4.n) in we infer 
that 

> eX,{H^,) - ■ .i^^^'f^-f^ „ ■ (7.1.6) 

Choosing e <^ guarantees that e||iJmt|| is small compared with \2{Hmove)- For 
this choice of e one gets 

X^iH) > X^m > eX2iH,s) - 0(eV). (7.1.7) 

Hence it suffices to prove that X2{Hcs) > n^'^^^\ 

Step 2. Recall that H^q = V'^HintV acts on the Dyck space "Hd- Its unique ground 
state \V) eI-Ld can be found by solving \M.n) = V It yields 



Let 7r(s) = (slr*)^ be the induced probability distribution on V. Given a pair of 
Dyck paths s,t eV, define 



n V vr s 



We claim that P describes a random walk on the set of Dyck paths V such that 
P{s,t) is a transition probability from s to t, and tt{s) is the unique steady state 
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of p. Indeed, since A/7r(t) is a zero eigenvector of Hcs, one has ^fP{s,t) = 1 and 
J2s = ^(^)- Off-diagonal matrix elements {s\Hes\t) get contributions only 

from terms — |00)(Zr| and — |/r)(00| in Hint- It implies that {s\Hcs\t) < for s ^ t 
and hence P{s,t) > 0. Furthermore, P{s,s) > 1/2 since {s\Hes\s) < n/2. In the 
Supplementary Material we shall prove the following. 

Lemma 2. Let s,t E V be any Dyck paths such that t can he obtained from s by 
adding or removing a single Ir pair. Then P{s,t) = Q{l/n^). Otherwise P{s,t) = 0. 



Let A2(-P) be the second largest eigenvalue of P. From Eq. (7.1.9) one gets 
^2{Hcs) = n{l — X2{P)). Hence it suffices to prove that the random walk P has a 
polynomial spectral gap, that is, 1 — \2{P) > n^'^^^\ 

Step 3. Lemma |2] tells us that P describes a random walk on a graph G = {T>, E) 
where two Dyck paths are connected by an edge, (s, t) e -E, iff s and t are related 
by insertion/removal of a single Ir pair. To bound the spectral gap of P we shall 
connect any pair of Dyck paths s, t G P by a canonical path 7(5, t), that is, a sequence 
So, Si, . . . , G such that Sq = s, Si = t, and {si, Sj+i) G E for all i. The canonical 
paths theorem |136l I137j shows that 1 — A2(-P) > l/(p/), where / is the maximum 
length of a canonical path and p is the maximum edge load defined as 

p = max , , ^, — -r n{s)n(t). (7.1.10) 

{a,b)GE 7T(a)P(a,b) ^ \ I \ l \ I 

The key new result that allows us to choose a good family of canonical paths is the 
following. 



Lemma 3. Let T)^ be the set of Dyck paths of length 2k. For any k > 1 there exists 
a map f : — ?■ T>k^i such that (i) the image of any path s G T>k can be obtained 
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from s by removing a single Ir pair, (ii) any path t G T>k_i has at least one pre-image 
in Vk, and (Hi) any path t G Pfc-i has at most four pre-images in Vk- 



The lemma allows one to organize the set of all Dyck paths T> into a supertree 
T such that the root of T represents the empty path and such that children of any 
node s are elements of f~^{s). The properties of / imply that Dyck paths of length 
2m coincide with level-m nodes of T, any step away from the root on T corresponds 
to insertion of a single Ir pair, and any node of T has at most four children. Hence 
the lemma provides a recipe for growing long Dyck paths from short ones without 
overusing any intermediate Dyck paths. It should be noted that restricting the 
maximum number of children to four is optimal since \Vk\ = Ck^ 4Vv/^A;3/2, where 
Cfc is the k-th Catalan number. Our proof of Lemma [3] based on the fractional 
matching method can be found in the Supplementary Material. Five lowest levels of 



the supertree T are shown on Fig. 7.1.2 



We can now define the canonical path 7(5, t) from s G T>m to t G V^.. Any 
intermediate state in 7(5, t) will be represented as uv where u G Vii is an ancestor 
of s in the supertree and v G P;" is an ancestor of t. The canonical path starts 
from u = s, V = ^ and alternates between shrinking u and growing v by making 
steps towards the root (shrink) and away from the root (grow) on the supertree. 
The path terminates as soon as u = and v = t. The shrinking steps are skipped 
whenever m = 0, while the growing steps are skipped whenever v = t. Note that any 
intermediate state uv obeys 

min(|s|, |t|) < \u\ + \v\ < max(|s|, |t|). (7.1.11) 
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Figure 7.1.2: (Color Online) Five lowest levels of the supertree T. Nodes of T are 
Dyck paths — balanced strings of left and right brackets. Depth-A; nodes are in 
one-to-one correspondence with Dyck paths of length 2k (the set P^). Any step 
towards the root requires removal of a consecutive [] pair. Any node has at most 
four children. The supertree can be described a family of maps / : — )■ 
such that f{s) is the parent of s. The maps / are defined inductively such that 
/([X]) = [/(X)], f{\]Y) = []fiY) for any node, fi[X]Y) = [fiX)]Y for black 
nodes, and /([X]F) = [X]/(y) for red (shaded) nodes. See the proof of Lemma |3] 
in the Supplementary Material for more details. 



Since any path 7(5, t) has length at most 2n, it suffices to bound the maximum edge 
load p. Fix the edge (a, h) E E with the maximum load. Let p(m, fc, /', /") be the 
contribution to p that comes from canonical paths '^{s,t) such that a = uv E 
where 

s e Vm, t eVk, u e Vi>, v e v^, 

and such that b is obtained from a by growing v (the case when b is obtained from a by 
shrinking u is analogous) . The number of possible source strings s G Vm contributing 
to p(m, k, I") is at most 4™~'' since s must be a descendant of u on the supertree. 
The number of possible target strings t G "Dfc contributing to p{m, k, I', I") is at most 
^k-i" gjm^jg f must be a descendant of v on the supertree. Taking into account that 
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7r(s) and nit) are the same for all s G and t eVj. we arrive at 



p{m, k, r, I") < 4" 



n{a)P{a,b) ni,+inP{a,h) 



with TT; = A'-{2^/\Mn\- Here we used the identity 7r(w) = and Eq. (7.1.8). 

Lemma [2] implies that 1/P(a,b) < rP^^\ Furthermore, the fraction of Motzkin 
paths of length n that have exactly 21 brackets is ai = Ci[^j) /\Ain\- However Ci ~ 
^ I yfixt^l'^ coincides with 4^ modulo factors polynomial in 1/n. Hence 

By definition, oi < \ for all I. Also, one can easily check that oi as a function 
of / has a unique maximum at / ^ n/3 and decays monotonically away from the 
maximum. Consider two cases. Case (1): I' + 1" is on the left from the maximum of 



(J;. From Eq. (7.1.11) one gets min (m, k) < I'+l" and thus o"mO"fc < cTmin (m,fc) < crr+«"- 



Case f^Sj.' /' + /" is on the right from the maximum of o"/. From Eq. (7.1.11) one 

gets max (m, k) > I' + I" and thus am^k < o"max(m,fc) < o-v+i"- both cases we get 
a bound p(m, /c, Z") < ra'^'^^^ Since the number of choices for m, fc, Z', /" is at most 
rf^^\ we conclude that p < rf'^^'^ and thus 1 — \2{P) ^ n^*^^"'^^ 

Open problems. Our work raises several questions. First, one can ask what is 
the upper bound on the ground state entanglement of FF spin-1 chains and whether 
the Motzkin state achieves this bound. For example, if the Schmidt rank x{L) for 
a block of L spins can only grow polynomially with L, as it is the case for the 
Motzkin state, ground states of FF spin-1 chains could be efficiently represented 
by Matrix Product States |138] (although finding such representation might be a 
computationally hard problem |139j ). One drawback of the model based on Motzkin 
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paths is the need for boundary conditions and the lack of the thermodynamic hmit. 
It would be interesting to find examples of FF spin-1 chains with highly entangled 
ground states that are free from this drawback. We also leave open the question 
of whether our model can indeed be regarded as a critical spin chain in the sense 
that its continuous limit can be described by a conformal field theory. Finally, 
an intriguing open question is whether long-range ground state entanglement (or 
steady-state entanglement in the case of dissipative processes) in ID spin chains can 
be stable against more general local perturbations, such as external magnetic fields. 



7.2 Supplementary Material 



7.2.1 Schmidt coefficients of the Motzkin state 



In this section we compute the Schmidt coefficients Pm defined in Eq. (7.1.5) and 
show that for large n and m one can use an approximation 

~ ''7^^ exp (— 3m^/n). (7-2.1) 

Let P„ A: ^ r}^"^'^ be the set of balanced strings of left and right brackets of length 
2n + k with k extra left brackets. More formally, s G T>n,k iff any initial segment of 
s contains at least as many /'s as r's, and the total number of /'s is equal to n + k. 

Lemma 4 (Andre's reflection method). The total number of strings in 'Dn,k is 

k + l f2n + k 



n + k -\- 1\ n 

Proof. For any bracket string s let L{s) and R{s) be the number of left and right 
brackets in s. Any s G {/,r}^"+'^ such that s ^ Vn^^ can be uniquely represented as 



s = urv, where r corresponds to the first unmatched right bracket in s, while m is a 
balanced string (Dyck path). Let v' be a string obtained from f by a reflection r -H- I 
and s' = urv'. Then 

L{s') = L{u) + L{v') = R{u) + R{y) = R{s) - 1 = n - 1 

and 

R{s') = R{u) + 1 + R{y') = L{u) + 1 + L{y) = L{s) + 1 
= n + k + 1. 



Furthermore, any string s' with n — 1 left brackets and n + k + 1 right bracket can 
be uniquely represented as s' = urv'. Hence the number of strings in Vn,k is 



n,k 



2n + k\ f2n + k\ k + l (2n + k 



n / \ n — 1 J n + k + 1 \ n 



□ 



One can easily check that \Cp^q{n)\ = |Co,p+q(ra)| since the identity of unmatched 
brackets does not matter for the counting. Let 



Mn,m = \Co,m{,n)\. 

Lemma |4] implies that 

\ - m + 1 ( 

^-^ t + m + 1 \ 



i>0 
2i+m<n 



m + 1 f n \ f2i + m 
i + m + l\2i + m J \ i 
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It can be rewritten as 



(7.2.2) 



where 



i>0 
2i+m<n 



fm + 1) ■ n\ 



{i + m + l)\i\{n — 2i — m)\ 
Let a = m/y/n and [3 = {i — n/?))ly/n. Using Stirling's formula one can get 



(7.2.3) 



3^3 
27rn3/2 



3"+^aexp (-3a2 - 9a/3 - 9/3^) . 



(7.2.4) 



We approximate the sum in Eq. (7.2.2) by integrating over i. Since z = + /S-y/n, 



we get di = y/nd[5, Since the maximum is near i = |, we can turn this sum into an 
integral from — oo to oo. The integral we need to evaluate is thus 



3^ ,n 



27rn3/2 



^3"+^a r e 

„+l 3^2 



2y^n 



3"+^ae-3° ~ mexp(-3m74n) 



Recalling that 



|Co,^(n/2)p 



(7.2.5) 



we arrive at Eq. (7.2.1) 
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7.2.2 Proof of Lemma d 



Suppose s EVk and t G T>k±i. Using the definition of P{s,t) one can easily get 



Here Mn[s] = {u e Mn ■ Dyck(u) = s}. Note tliat {u\Hint\v) = -1/2 if u 
and V differ exactly at two consecutive positions where u and v contain 00 and Ir 
respectively or vice versa. In all other cases one has {u\Hint\v) = 0. 

Suppose t G T>k+i and P{s, t) > 0. Let us fix some j G [0, 2k] such that t can be 
obtained from s by inserting a pair Ir between Sj and Sj+i. For any string u G A^ni-s] 
in which Sj and Sj+i are separated by at least two zeros one can find at least one 
V G A^n[i^] such that {u\Hint\v) = —1/2. The fraction of strings u G A^n[s] in which 
Si and separated by two or more zeros is at least 1/n^ which implies 



Suppose now that t G V^-i. Let us fix some j G [1,2A; — 1] such that t can be 
obtained from s by removing the pair SjSj^i = Ir. For any string u G A^„,[s] in 
which Sj and sj+i are not separated by zeros one can find at least one v G A^n[^] 
such that {u\Hint\v) = —1/2. The fraction of strings u G A^nl-s] in which sj and Sj+i 
are not separated by zeros is at least 1/n which implies 




P{s,t) > 



1 

2^' 



P{s,t) > 



1 

2^' 



191 



7.2.3 Proof of Lemma M 

Let us first prove a simple result concerning fractional matchings. Consider a bipar- 
tite graph G = {AU B,E). Let x = {xejees be a vector of real variables associated 
with edges of the graph. For any vertex u let 5{u) be the set of edges incident to u. 
Define a matching polytope |130] 

V = {x : Xe>0 for all e e E, 

1< ^ Xe<4:, ^ Xe = l 
ee5(a) ee<5(fe) 

for all a G v4 and b G B}. 

Lemma 5. Suppose V is non-empty. Then there exists a map f : B ^ A such that 
(i) f{h) = a implies {a,b) G E, (ii) any vertex a G A has at least one pre-image in 
B, and (Hi) any vertex a E A has at most four pre-images in B. 

Proof. Since V is non-empty, it must have at least one extremal point x* G V. Let 
E* <Z E he the set of edges such that x* > 0. We claim that E* is a forest (disjoint 
union of trees). Indeed, suppose E* contains a cycle Z (a closed path). Then x* ^ < 1 
for all (a, h) & Z since otherwise the cycle would terminate at h. Hence < x* < 1 
for all e E Z. Since the graph is bipartite, one can label edges of Z as even and odd 
in alternating order. There exists e 7^ such that x* can be shifted by ±e on even 
and odd edges of Z respectively without leaving V. Hence x* is a convex mixture 
of two distinct vectors from V. This is a contradiction since x* is an extreme point. 
Hence E* contains no cycles, that is, E* is a forest. We claim that x* = 1 for all 
e e E*. Indeed, let T C E* he the subset of edges with < x* < 1. Obviously, 
T itself is a forest. Degree-1 nodes of T must be in A and there must exist a path 
'J T starting and ending at degree-1 nodes u, u' G A. Since < x* < 1 for all 
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e G 7) there exists e 7^ such that x* can be shifted by ±e on even and odd edges 
of 7 respectively without leaving V. This is a contradiction since x* is an extreme 
point. Hence = 1 for all e G -E*. We conclude that x* G {0, 1} for all edges of G. 
The desired map can now be defined as f{h) = a iff x*^ = 1. □ 

We are interested in the case where 

A = Vn-i and B = Vn 

are Dyck paths of semilength n — 1 and n respectively. Paths a G "Pn-i and b G Vn 
are connected by an edge iff a can be obtained from b by removing a single Ir pair. 
Our goal is to construct a map / : Vn — )■ I^n-i with the properties (i),(ii),(iii) stated 
in Lemma |5] According to the lemma, it suffices to choose / as a stochastic map. 
Namely, for any b G Vn we shall define a random variable G Vn^i with some 
normalized probability distribution. It suffices to satisfy two conditions: 

Pr[/(6) = a] > only if a can be obtained (7.2.6) 
from b by removing a single 
Ir pair, 

and 

Pr[/(6) = a] = X„ for all a G P„-i. (7.2.7) 

Here 1 < X„ < 4 is some function of n that we shall choose later. We shall define / 
using induction in n. 

Lemma 6. For every n > 1 there exists 1 < Xn < 4 and a stochastic map f : Vn — > 
Vn-i satisfying Eqs. fTJ^djjJ^. 
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Proof. Any Dyck path b G P„ can be uniquely represented as 6 = Isrt for some 
s G Vi, t G Vn-i-i, and i G [0,n — 1]. We shall specify the map / : Vn — )■ I?n-i by 
the following rules: 



h G Vn 


fib) G Vn-i 


probability 


Isrt, s E Vi, 1 < i < n — 2 


lf{s)rt 


Pi 


Isrt, s (zVi, 1 < i < n — 2 


Isrfit) 


i-Pi 


Irt, t G Vn-i 


Irfit) 


1 


Isr, s G Vn-i 


lf{s)r 


1 



Here we assumed that / has been already defined for strings of semilength up to 
n — 1 such that Eqs. (7.2.6|7.2.7) are satisfied. By abuse of notation, we ignore 



the index n in /, so we regard / as a family of maps defined for all n. It is clear 



that our inductive definition of / on Vn satisfies Eq. (7.2.6). The probabilities 
Pi, ■ ■ . ,Pn~2 € [0, 1] are free parameters that must be chosen to satisfy Eq. (7.2.7). 
Note that these probabilities also implicitly depend on n. The choices of f{b) in 
the first two lines of the above table are represented by black and red nodes in the 



example shown on Fig. 7.1.2 Consider three cases: 

Case 1: a = Irt' for some t' G I?n-2- Then f{b) = a iff 6 = llrrt' or 6 = Irt for 
some t G Vn-i such that f{t) = t' . These possibilities are mutually exclusive. Hence 



Y,^v[f{h) = irt'] = p,+ mit) = t'] 

= Pi+Xn-l. 



Substituting it into Eq. (7.2.7) gives a constraint 



Pi 



Xn — X, 



n-1- 



(7.2. 
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Case 2: a = Is'r for some s' G 2^„-2- Then f{b) = a iS b = Is'rlr oi b = Isr 
for some s G I?n-i such that f{s) = s'. These possibihties are mutually exclusive. 
Hence 

5^ Pr[/(6) = ZsV] = l-pn-2+ Yl Pr[/(s) = s'] 

= l-p„_2+X„_i. 



Substituting it into Eq. (7.2.7) gives a constraint 



Pn-2 — 1 — (yXn — 



n-1 



(7.2.9) 



It says that X„ must be a non-decreasing sequence. 



Case 3: a = Is'rt' for some s' G Vi, t' G Vn-i-2, and i = 1, . . . ,n — 3. In other 
words, both s' and t' must be non-empty. Then /(6) = a iff 6 = Isrt' for some 
s G "Pj+i such that /(s) = s', or 6 = ZsVt for some t G "Dn-j-i such that f{t) = t' . 
These possibilities are mutually exclusive. Hence 



Pr[/(6) = Is'rt'] = p.+i = ^1 



Pi+i Xi+i + (1 - Pi)X„_i„i. 



Substituting it into Eq. (7.2.7) gives a constraint 



71— i— 1 



(7.2.10) 
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for each i = 1, . . . , n — 3. Let us choose 



4(2 - 1/2) 



i + 1 



(7.2.11) 



Combining Eqs. ( 7.2.8[[7^.9[7.2.10 ) we obtain a hnear system with unknown variables 
Pi, . . . ,Pn-2 G [0, 1]. We shall look for a solution {pi} having an extra symmetry 



Pi + pn-i-i = 1 for i = 1, . . . , n - 2. 



(7.2.12) 



One can check that the system defined by Eqs. ( |7.2.8|7.2.9[7.2.10[7.2.12D has a solu- 
tion 

i{i + l)(3n -2i-l) 



Pi 



,n-2. 



(7.2.13) 



n{n + l){n — 1) 

Hence we have defined the desired stochastic map / : Vn — )■ T>n~i- This proves the 
induction hypothesis. 

It remains to note that for n = 1,2 the map / is uniquely specified by Eqs. ( 7.2.6[7.2.7 ) 
and our choice of Xn- Indeed, one has V2 = {llrr, Irlr}, Vi = {Ir}, and Vq = 0. To 



satisfy Eq. (7.2.6), we have to choose f{llrr) = f(lrlr) = Ir for n = 2 and f{lr) 



for n = 1. It also satisfies Eq. (7.2.7) since X2 = 2 and Xi = 1. This proves the base 
of induction. □ 



7.2.4 Ground state energy for unbalanced subspaces 

Recall that the unbalanced subspace Hp^q is spanned by strings s G Cp^g that have 
p unmatched right and q unmatched left brackets. Our goal is to prove that the 
restriction of H onto any subspace Hp^g with p > or g > has ground state energy 
at least n^'^^^\ By the symmetry, it suffices to consider the case p > 0. To simplify 
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the analysis we shall omit the boundary term | / ) (/ 1 „ . Note that such omission can only 
decrease the ground state energy. Accordingly, our simplified Hamiltonian becomes 

n-1 

H^\r){r\, + J2^,,j+i- (7-2.14) 

Recall that 11 is a projector onto the subspace spanned by states |00) — \lr), \Ql) — 1/0), 
and I Or) — |rO). Let Ai(i7) be the ground state energy of H. 

Any string s G Cp^g can be uniquely represented as 

s — Uoruiru2 ■ ■ ■ rUplvilv2 . . .Ivg 

where Ui and Vj are Motzkin paths (balanced strings of brackets) . The remaining p 
right and q left brackets arc unmatched and never participate in the move 00 ^ Ir. 
It follows that the unmatched brackets can be regarded as "solid walls" that can be 
swapped with O's but otherwise do not participate in any interactions. In particular, 
the spectrum of H restricted to Hp^q depends only on p + q as long as p > 0. This 
allows us to focus on the case q — 0, i.e. assume that all unmatched brackets are 
right. 

Given a string s e Cp^, let s e {0, 1, r, x, y}" be the string obtained from s by the 
following operations: (i) replace the first unmatched right bracket in s by 'x\ and 
(ii) replace all other unmatched brackets in s (if any) by 'y'. Define a new Hilbert 
space Hp whose basis vectors are |s), s e Cp^. Consider a Hamiltonian 

n-1 

H^\x){xU + '£^jj^i + ^j,j+i + ^lj+i, (7-2.15) 

where 0^ and are projectors onto the states \0x) — \xO) and \0y) — \yO) respectively 
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(with a proper normalization). One can easily check that {s\H\t) = {s\H\t) for any 
s,t E Cpfl. Hence the spectrum of H on l-Lpfl coincides with the spectrum of H on 
1-Lp. Furthermore, if we omit all the terms ©jj+i in H, the ground state energy can 
only decrease. Hence it suffices to consider a simplified Hamiltonian 

n—1 

= \x){x\, + J]H,„+i + Ql^^, (7.2.16) 

which acts on Hp. Note that positions of y-particles are integrals of motion for H^. 
Moreover, for fixed positions of y-particles, any term in touching a y-particle 
vanishes. Hence can be analyzed separately on each interval between consecutive 
y-particles. Since our goal is to get a lower bound on the ground state energy, we 
can only analyze the interval between 1 and the first y-particle. Equivalently, we can 
redefine n and focus on the case p — 1, q — that is, assume that there is only one 
unmatched right bracket. The relevant Hilbert space l-Li is now spanned by states 

|s) ® \x) (g) where s e Mj-i, t e Mn-j- 

Recall that A^fc is the set of Motzkin paths (balanced strings of left and right brack- 
ets) of length k. 

We would like to treat the terms responsible for the motion and detection of the 
x-particle as a small perturbation. To this end, choose any < e < 1 and define the 
Hamiltonian 

n—1 n—1 

Clearly, Hf < H^, so it suffices to get a lower bound on the ground state energy of 
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Let us first find the ground subspace and the spectral gap of the unperturbed 
Hamiltonian Hq = Yl^Zi^j,j+i- Note that the position of the x-particle j is an 
invariant of motion for Hq. Moreover, any projector Ilj^j+i touching the x-particle 
vanishes. Hence we can analyze separately on the two disjoint intervals A = 
— 1] and B = [j + l,n]. It follows that the ground subspace of Hq is spanned 
by normalized states 

= \Mj-i)®\x)®\Mn-j), j = l,...,n. (7.2.17) 



The spectral gap of Hq can also be computed separately in A and B. Since we have 



already shown that the original Hamiltonian Eq. (7.1.2) has a polynomial gap inside 
the Motzkin subspace, we conclude that X2{Hq) > n^'^^^K 



Let us now turn on the perturbation. The first-order effective Hamiltonian acting 
on the ground subspace spanned hj ipi, . . . ,ipn describes a hopping of the x-particle 
on a chain of length n with a delta-like repulsive potential applied at site j = 1. 
Parameters of the hopping Hamiltonian can be found by calculating the matrix 
elements 



and 



where = \Aik\ is the k-th Motzkin number. We arrive at the effective hopping 
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Hamiltonian acting on C", namely, 



H,s = \l){l\ + Y.^,j+i, (7.2.18) 



where 



+ + b + (7.2.19) 
is a rank-1 projector. Applying the Projection Lemma of |129j we infer that 



where V = |a;)(a;|i + Yl^=i perturbation operator. Since X2{Hq) > 

n~^^^\ we can choose e polynomial in 1/n such that 2e||y|| is small compared with 
\2{Hq). For this choice of e one gets 

Ai(i/f)>eAi(i/eff)-0(e^)nO(i). 

Hence it suffices to show that Xi{Hcs) > n^'^''^\ where Hes is now the single x- 



particle hopping Hamiltonian Eq. (7.2.18). 



Let us first focus on the hopping Hamiltonian without the repulsive potential: 

n-l 
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This Hamiltonian is FF and its unique ground state is 

n 

\g)^J2VMj-iM^-,\j). (7.2.20) 

Our strategy will be to bound the spectral gap of Hmove and apply the Projection 
Lemma to i^eff by treating the repulsive potential |1)(1| as a perturbation of Hmove- 
First let us map Hmove to a stochastic matrix describing a random walk on the 
interval [l,n] with the steady state 7r(j) = For any a, 6 e [l,n] define 



P{j,k) = 5,-, - {j\Hmove\k)^'^. (7.2.21) 

Since a/ vr(j) is a zero eigenvector of H^ove, "we infer that ^f.P{j,k) = 1 and 
^(j)-P(j) ^) = 7r(/c). A simple algebra shows that 

P{jj + l)^Mllzi^ and P{j + l,j)^^ 

are the only non-zero off-diagonal matrix elements of P. We shall use the following 
property of the Motzkin numbers. 

Lemma 7. For any n > 1 one has 1/3 < Mn/Mn+i < 1. Furthermore, for large n 
one can use an approximation 

Mn (7.2.22) 

where c ~ 1.46. 

The lemma implies that 

l<PUJii)<l 
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for all j. Hence the diagonal matrix elements P{j,j) are non-negative, that is, we 
indeed can regard P{j, k) as a transition probability from j to k. Furthermore, using 



Eq. (|7.2.20|) and the above lemma we infer that the steady state vr is 'almost uniform', 

<![(^<^o(i) foralll<7,A;< 



that is. 



n 



-o(i) 



(j) 



n. 



(7.2.23) 



In particular, minj7r(j) > n We can now easily bound the spectral gap of P. 

For example, applying the canonical paths theorem stated above we get 1 — X2{P) > 



l/{pl) where p is defined in Eq. (7.1.10) and the canonical path 7(s,t) simply moves 



the x-particle from s to t. Since the denominator in Eq. (7.1.10) is lower bounded 
by ra-^(^), we conclude that 1 - A2(P) > n'^'^^'l It shows that \2{Hmove) > n'^^^^l 
To conclude the proof, it remains to apply the Projection Lemma to Hes defined 



in Eq. (7.2.18) by treating the repulsive potential |1)(1| as a perturbation. Now the 
effective first-order Hamiltonian will be simply a c-number {l\gy = 7r(l) > n~'~'^^^ 
which proves the bound Ai(ifcfr) > n^'^^^\ 
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